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INTRODUCTION 

Since  the  first  operational  airtanker  drops  in  1956, 
many  different  aircraft,  ranging  from  agricultural  spray 
planes  to  modem  military  transports,  have  been  used  as 
airtankers.  The  tank  and  gating  systems  of  these  aircraft 
are  of  many  different  designs,  with  differing  capabilities. 
In  1969,  the  Forest  Service  initiated  research  studies  to 
quantify  the  capabilities  of  different  tank  and  gating 
systems.  The  studies  entailed  the  dropping  of  retardant 
and/or  water  over  a  sampling  grid  and  determining  the 
groimd  pattern  under  a  variety  of  conditions:  tank  con- 
figuration, door  sequencing  speed,  drop  height,  retardant 
type,  relative  humidity,  temperature,  windspeed,  and 
direction  (George  1975;  George  and  Blakely  1973).  Using 
these  data,  Honeywell  Corporation,  under  contract  to  the 
Forest  Service,  developed  a  pattern  simulation  model 
whose  primary  variable  was  the  flow  history  from  the 
individual  systems  (Swanson  and  others  1975).  The 
rough  model  used  flow  history  derived  from  motion 
pictures  of  airborne  airtanker  drops.  The  model  was  re- 
fined through  the  use  of  accurately  measured  flow  data. 
(Swanson  and  others  1977).  A  BASIC  listing  of  the  model 
is  included  as  appendix  A. 

Information  from  the  refined  model  (PATSIM)  sug- 
gested that  proper  tank  design  could  be  used  to  produce 
drop  patterns  that  would  be  effective  for  nearly  any  fuel- 
fire  condition.  In  order  to  define  the  tank  parameters  and 
to  fine-tune  and  validate  PATSIM,  an  experimental  tank 
and  gating  system  (ETAGS)  was  designed,  installed  in  a 
Forest  Service  P2V-5  airtanker,  and  drop-tested  over  a 
sampling  grid.  The  ETAGS  was  designed  so  that  tank 
and  gating  parameters,  including  load  size,  fluid  head 
height,  door  opening  length-to-width  ratio,  compartment 
separation,  door-opening  area,  and  fluid  flow  rates  could 
be  examined  along  with  flight  parameters  (drop  height 
and  speed).  The  knowledge  gaiined  from  the  ETAGS  test 
was  used  to  refine  PATSIM  to  its  final  form.  The  accu- 
racy of  PATSIM  was  proven  for  tanks  with  "average" 
length-to-width  ratios;  however,  for  tanks  with  an  ex- 
treme length-to- width  ratio  PATSIM  showed  a  strong 
tendency  to  underpredict  the  length  of  the  ground  pattern 
and  overpredict  the  retardant  coverage  concentrations 
(Swanson  and  others  1978). 

About  90  percent  of  the  current  airtankers  have  been 
static-tested  and  user  guidelines  produced.  New  airtank- 
ers are  tested  as  part  of  Interagency  Airtanker  Board 
(lATB)  approval  process.  PATSIM  is  also  used  as  a  tool  to 
aid  in  the  design  of  new  tank  and  gating  systems  and  for 
evaluating  flow-rate  modifications  to  enable  airtankers  to 
meet  improved  and  more  flexible  standards  set  by  the 
IATB(USDAFS  1987). 


GUIDELINE  DEVELOPMENT 
SEQUENCE 

The  first  step  in  the  generation  of  airtanker  perform- 
ance guidelines  is  to  quantify  the  flow  of  water/retardant 
from  the  airtanker  delivery  system  for  all  possible  drop 
configurations.  This  is  accomplished  by  static  testing 
airtankers  to  determine  flow  history,  door-opening  rate, 
and  internal  tank  pressure  according  to  the  methods  set 
forth  by  Blakely,  George,  and  Johnson  (1982).  The  repre- 
sentative flow  rates  for  each  compartment  and  drop  type 
are  selected  and  used  as  input  to  PATSIM  (the  airtanker 
ground  pattern  response  model)  and  PATADD  (adjunct  to 
PATSIM  for  reducing  data  from  dissimilar  compartments 
in  a  tanking  system).  The  ground  pattern  data  output  by 
PATSIM/PATADD  is  smoothed,  tabulated,  plotted  (using 
computer  programs  developed  at  the  Intermountain  Fire 
Sciences  Laboratory  and  written  in  BASIC),  and  as- 
sembled into  a  standard  guideline  format.  A  schematic  of 
the  steps  followed  in  producing  airtanker  performance 
guidelines  is  shown  in  figure  1. 


Measure  and  static  test  aircraft 
using  computer  to  collect  data. 


♦ 


Develop  algorithm  using  tank  measurements  and 
calculate  flow,  door  opening,  and  pressure  histories. 


Resolve  any  discrepancies, 
i.e.,  negative  pressure. 


Input  representative  flows  to  PATSIM. 

♦ 

Smooth  and  tabulate  PATSIM  generated  data. 


Using  inhouse  software  produce  camera  ready  copy  of: 
Best  strategy       Characterisics,       Footprints  Text 
charts  flow,  technical        and  width 

data,  tables  charts 


Assemble  into  guideline  format. 


Figure  1 — Guideline  production  sequence. 
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OPERATION  OF  PATSIM/PATADD 
MODELS 

The  simulation  program  (PATSIM)  uses  an  algorithm 
consisting  of  four  distinct  parts: 

1.  Generation  of  a  marginal  downrange  distribution 
(percent  of  retardant  delivered  in  each  foot  of  pattern 
length)  from  flow  rate  and  aircraft  velocity. 

2.  Development  of  the  marginal  crossrange  distribution 
as  a  function  of  altitude  and  retardant  type. 

3.  Calculation  of  the  ground  response  pattern  from  the 
distributions. 

4.  Summarization  of  the  ground  response  pattern  as  a 
scaled  pattern  schematic  of  the  coverage  level  distribution 
on  the  ground  and  the  gallons,  area,  and  lengths  associ- 
ated with  various  coverage  levels. 


The  generation  of  the  marginal  range  distribution, 
though  empirically  derived,  has  its  analogy  in  the  physi- 
cal process  of  the  breakup  of  the  mass  of  retardant  re- 
leased from  an  eiirtanker  (fig.  2).  The  door  opening  shapes 
the  emerging  fluid  (phase  A),  which  is  subject  to  small 
particle  stripping,  drag  deceleration,  and  amplifying  in- 
stabilities of  two  types — Taylor  instabilities  on  the  front 
surface  and  Helmholtz  instabilities  on  the  side  surfaces. 
The  small  particles  stripped  from  the  retardant  are  a 
relatively  insignificant  portion  of  the  total  mass.  The 
emerging  liquid  mass  also  expands  to  the  sides  and 
becomes  shorter  as  lateral  flow  replaces  some  of  the 
deceleration  (phase  B).  The  amplitude  of  the  Taylor  in- 
stabilities increases  until  air  pressure  on  the  front  surface 
causes  breakup  and  final  particle  size  distribution  (phases 
C  and  D).  Rheological  properties  and  airspeed  establish 
the  time  to  breakup  and  final  particle  size. 


STEADY- 
STATE 
PHASE 


TANK 


DOOR 
OPENING 


DOOR 
OPEN 


TANK  EDGE 
EFFECTS 


ONSET  OF 
TAYLOR 
INSTABILITIES 


ONSET  OF 
STRIPPING 


k  THIN  SECTION  BREAKUP  DUE  TO 
\  HELMHOLTZ  INSTABILITY 


THIN  SECTION 

BREAKUP 


■CONTINUED 
STRIPPING 


O  o  o 


POINT  D 

PATTERN  COMPOSITION 
AFTER  BREAKUP 


CONTINUED  THIN 
SECTION  BREAKUP 


^     BREAKUP  OCCURS 


Figure  2 — Schematic  representation  of  the  breakup  process. 
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In  the  case  of  one  or  more  compartments  located  aft  of, 
and  released  at  the  same  time  as,  a  front  compartment, 
the  forward  released  retardant  acts  as  a  wind  shield  for 
the  aft  releases,  thereby  delaying  the  breakup  until  the 
forward  release  is  nearly  at  its  final  particle  size.  At  that 
point,  the  aft  release  passes  through  the  retardant  cloud 
of  the  forward  release. 

The  algorithm  first  calculates  the  range  marginal  dis- 
tribution by  breaking  the  volume  discharged  from  the 
tank  into  small  increments  based  on  time  and  calculating 
each  one's  distribution  at  the  ground  level.  The  distribu- 
tions are  summed  to  produce  the  final  range  marginal 
distribution. 

To  develop  the  range  marginal  distribution  into  an 
isopleth,  a  normal  crossrange  marginal  is  assumed  and 
is  calculated  on  the  basis  of  altitude  and  retardant  type. 
An  amount  of  retardant  lost  to  evaporation  and  small 


particle  drift  is  discarded.  A  percentage  of  the  remaining 
retardant  (also  a  function  of  altitude  and  retardant  type) 
is  distributed  over  an  elliptical  area.  This  quantity  repre- 
sents the  fringe  of  the  pattern,  contains  low  coverage 
levels,  and  is  most  subject  to  wind  drift.  The  remaining 
70-90  percent  of  the  retardant  is  distributed  as  a  normal 
distribution.  The  sum  of  these  two  sections  produces  the 
simulated  ground  pattern  (fig.  3). 

In  order  to  calculate  the  most  efficient  time  interval 
between  single  tank  releases  that  make  up  a  trail  drop, 
a  subroutine  (PATADD)  is  used.  PATADD  shifts  and 
overlays  the  simulated  patterns,  determining  the  time 
interval  between  door  releases  for  each  coverage  level  that 
produces  the  longest  length  of  retardant  pattern  for  that 
coverage  level. 


VOLUME 


]TRAILOFF 


DOOR 
,  OPENING 
i  TIME  TO 
ESTABLISH 
FLOW 


(1) 


The  Program  Increments 
Flow  Rate 


(2)  Flys  Each  Packet 
to  Extinction 


(3) 


Distributes 
the  Fluid 
in  Space 


(4)  Sums  the  Distribution 
as  a  Ground  Marginal 


RANGE 


(5)  Adjust  for  Altitude 
Losses  (%  R'cvd) 


RANGE 

(6)  Removes  a  Uniform 
Fringe  Pattern 


 ^  ^\ 


SEMICROSS 
RANGE 


(7)  Develops  the  Two  Parts  of 
the  Pattern  as  a  Function 
of  Altitude 


Figure  3 — Simulation  overview  (Swanson  and  others  1975). 
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Inputs  to  PATSIM  Model 

The  simulation  is  designed  to  be  used  in  two  general 
modes:  (1)  the  tank  is  described  geometrically  and  the 
flow  rate  calculated  and  (2)  the  flow  rate  is  determined 
experimentally  and  input  directly  to  the  simulation. 
The  latter  has  proven  to  be  the  most  accurate  and  is  used 
to  produce  the  guidelines.  Experience  has  shown  that 
door-opening  rate  (which  is  reflected  in  the  flow  rate)  and 
ejection  velocity  from  the  tank  have  little  effect  on  the  ac- 
curacy of  the  prediction  and  are  therefore  not  used. 

The  inputs  used  in  PATSIM  for  guideline  production 
are:  (1)  experimental  tank  flow  data;  (2)  quantity  vs.  time; 
(3)  number  of  flow  data  entries;  (4)  total  gallons  in  tank; 
(5)  number  of  tanks  released;  (6)  retardant  type; 
(7)  aircraft  speed;  (8)  aircraft  altitude;  (9)  ADD,  a  flag  set 
if  trail  drops  are  to  be  calculated  using  the  PATADD 
subroutine  (table  1).  A  sample  input  is  contained  in 
appendix  B. 

PATSIM  Modification 

The  PATSIM  predicts  ground  response  patterns  satis- 
factorily for  most  airtankers  in  use.  But  the  ETAGS 
study  showed  that  modifications  to  PATSIM  were  neces- 
sary for  a  tank  with  very  clean  exit  geometry.  One  par- 
ticular tank,  installed  in  a  Conair  Aviation,  Ltd.,  DC-6, 
had  compartments  that  were  smooth,  with  no  reinforcing 
members  or  door  edges  that  could  cause  turbulence  in  the 
exiting  retardant.  Using  data  for  actual  ground  response 
obtained  over  a  test  grid,  it  was  found  that  an  erosion 
constant  Kl  and  a  shape  factor  DELZ  need  adjustment  to 
correctly  predict  the  down  range  marginal.  Similarly, 
SIGMA,  the  shape  factor  for  the  normal  distribution  in 
the  cross  range  marginal  also  needed  adjustment. 

After  using  various  values  for  these  constants  and  com- 
paring the  results  from  PATSIM  with  the  actual  pattern 
responses  from  drop  tests,  it  was  found  that  the  three 
constants  worked  best  for  generally  smooth  tanks  when 
converted  to  variables  that  were  a  function  of  drop  height 
and  the  number  of  compartments  dropped.  Only  one  such 
tank  has  been  tested.  PATSIM  was  developed  from 


ground  response  data  for  tanks  that  would  impart  turbu- 
lence to  the  exiting  retardant  due  to  internal  bracing,  door 
sills,  and  other  features. 

Program  Output 

The  program  output  consists  of  (1)  a  report  of  input 
data;  (2)  detailed  pattern  plot;  (3)  contour  interval  sum- 
maries; (4)  a  scaled  pattern  plot;  (5)  pattern  length  data; 
and  (6)  PATADD  output,  if  selected.  A  sample  output  is 
shown  in  appendix  B. 

INPUT  DATA 

An  echo  of  input  data  for  confirmation  of  data  and  iden- 
tification of  the  pattern. 

DETAILED  PATTERN  PLOT 

A  printout  of  the  pattern  on  a  28-column  by  56-row 
unsealed  array.  The  distance  between  columns  is  15  feet 
and  between  rows  is  30  feet. 

CONTOUR  INTERVAL  SUMMARIES 

1.  Table  of  the  amounts  of  retardant  in  each  coverage 
level  (gallons  per  100  square  feet)  0.0-0.5,  0.6-1.0,  1.1-2.0, 
2.1-3.0,  3.1-4.0,  4.1-5.0,  5.1-6.0,  6.1-7.0,  7.1-8.0,  8.1-9.0, 
9.1-10.0,  and  greater  than  10  gallons. 

2.  Total  amount  of  retardant  on  the  grid  in  gallons  and 
the  percent  of  total  gallons  dropped  that  is  recovered. 

3.  The  area  in  square  feet  covered  tabulated  by  the 
coverage  levels  and  the  total  area  covered. 

4.  A  scaled  pattern  schematic  with  hsted  coverage 
levels  and  a  table  of  pattern  lengths  tabulated  by  cover- 
age levels. 

5.  If  the  PATADD  option  was  selected,  a  table  of  the 
maximum  length  of  the  patterns  for  multiple  tank  drops 
and  the  delay  between  releases  required  to  obtain  the 
patterns  at  coverage  levels  equal  to  or  greater  than  0.5, 
1.0,  2.0,  3.0,  4.0,  6.0,  8.0,  10.0  gallons  per  100  square  feet 
(gpc). 


Table  1— PATSIM  inputs 

Variable  No. 

Name 

Description 

17 

GAL! 

Total  volume  of  tank  in  gallons 

18 

CFLAG 

Number  of  tanks  released 

24 

RDQV 

Flag  to  select  option  for  entry  of  experimental  flow  data  also 

the  number  of  entries  in  the  flow  data  table 

45 

RFLAG 

Retardant  type — 1  =  gum-thickened;  2  =  waterlike 

51 

VAGI 

Aircraft  speed 

2 

ALT 

Aircraft  altitude 

53 

ADD 

Flag  to  select  PATADD 

* 

* 

A  code  in  columns  1-8  indicating  the  number  of  patterns  to  add 
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GUIDELINE  PRODUCTION 
Selection  of  Inputs 

The  basic  input  will  be  the  volume  discharged  as  a 
function  of  time  (flow)  for  each  compartment  of  the  tank 
system  to  he  analyzed.  For  example,  a  six-tank  system 
might  have  its  six  compartments  in  a  row  across  the  tank, 
with  the  two  outside  tanks  differing  from  the  four  inner 
tanks  (fig.  4).  Flows  from  the  outer  and  inner  tanks 
would  be  expected  to  differ  even  though  their  volumes  are 
nearly  identical.  It  is  also  possible  that  due  to  capacities 
of  hydraulic  or  pneumatic  door-opening  systems  the  doors 
of  some  compartments  may  open  at  different  rates  than 
others.  This  produces  different  flows  from  the  same  com- 
partment when  the  compartment  is  released  alone  or  in 
conjunction  with  one  or  more  other  compartments.  The 
static  test  data  (flow  rates,  door-opening  time,  and  tank 
pressure)  should  include  tests  of  all  geometrically  dissimi- 
lar tanks  under  all  possible  release  conditions.  In  the 
case  of  the  tank  in  figure  5,  if  each  compartment  opened 
at  the  same  rate  for  1,  2,  3,  or  6  compartments  released  at 


0.1      0.2     0.3     0.4      0.5     0.6     0.7     0.8  0.9 


TIME  (sec) 
Figure  4 — P2V-7  tank  and  flow  rates. 


TANK  APPEARANCE  DURING  LAST  PORTION  OF  RELEASE 


SYNTHESIZED  FLOW 

200 1 


SIMILAR  TO  RESULTS  OBTAINED 
IN  29-INCH  UNDER-VENTED  TUBE. 


Figure  5 — Method  of  inferring  final  stages  of  flow 
from  an  undervented  tank. 


a  time,  two  different  flows  (outer  and  inner  compart- 
ments) would  be  needed  to  characterize  the  tank.  If  the 
flows  had  been  different  for  one  compartment  at  a  time 
and  six  compartments  at  a  time,  four  different  flows 
would  be  required.  The  best  method  of  comparing  flows 
is  to  plot  volume  discharged  vs.  time  from  static  tests  for 
any  compartments  that  could  be  different  (fig.  5).  In  cases 
where  a  difference  is  not  obvious,  PATSIM  can  be  applied 
using  the  flow  data  to  be  compared  as  inputs  to  determine 
if  pattern  differences  exist. 

In  the  case  of  inadequate  venting  of  a  compartment,  the 
flow  may  appear  to  halt  or  slow  momentarily  before  con- 
tinuing. Because  flow  measurements  are  taken  by  use  of 
a  float  on  the  surface  of  the  liquid,  it  would  appear  that 
the  flow  out  of  the  tank  was  halted.  But,  examination  of 
high-speed  films  and  testing  with  a  transparent  tube 
showed  that  the  tank  was  still  emptying  from  the  bottom. 
In  this  case,  the  flow  must  be  synthesized.  The  best 
approximation  of  flow  is  to  extrapolate  the  volume  dis- 
charged to  an  end  point  where  the  float  begins  to  free-fall 
(fig.  5).  With  undervented  tanks,  the  flow  and  thus  the 
groimd  pattern  can  exhibit  undesirable  effects  such  as 
thinning  or  gaps  in  coverage.  Whenever  possible,  adjust- 
able venting  systems  are  used  and  venting  is  fixed  after 
optimum  flow  rate  is  attained.  When  venting  adjust- 
ments cannot  be  adjusted  at  the  time  of  the  test,  vent  size 
is  set  on  the  basis  of  experience. 
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After  evaluation  of  the  flow  data,  typical  flows  are  se- 
lected for  use  in  the  computer  simulation  of  drop  patterns. 
If  there  is  a  large  amount  of  variation  in  flows  within 
replicate  static  tests,  a  mean  flow  should  be  selected.  But 
the  extremes  of  flow  data  should  also  be  evaluated  using 
the  simulation  to  assess  the  sensitivity  of  these  variations 
on  guideline  conclusions.  In  most  cases,  these  variations 
are  minor. 

The  computer  simulations  are  run  at  a  midrange 
aircraft  velocity  taken  from  the  velocities  listed  in  the 
approved  supplemental  flight  certificate  when  used  as  an 
airtanker.  This  range  is  typically  narrow  (110-135kt)  and 
ground  pattern  effects  are  relatively  small  compared  to 
drop  height,  drop  configuration,  or  release  intervals. 
Examination  of  the  sensitivity  of  patterns  to  aircraft  ve- 
locity suggests  that  differences  due  to  velocity  over  the 
range  of  typical  drops  is  less  than  might  be  expected. 
Higher  velocities  generally  reduce  peak  coverage  and 
increase  pattern  length  at  low  coverage  levels.  Slower 
aircraft  velocities  tend  to  increase  maximum  coverage 
levels  at  the  expense  of  the  lower  coverage  pattern 
lengths.  The  uncertainty  of  true  airspeed  at  the  time  of 
the  drop  is  generally  less  than  the  magnitude  of  pattern 
change  resulting  from  the  selection  of  a  specific  aircraft 
velocity.  Performance  in  an  actual  drop  may  be  somewhat 
better,  but  seldom  worse  than  the  simulation  shown  in 
the  guides.  Level  flight  over  flat  terrain  is  assumed  due 
to  the  very  large  number  of  graphs  needed  to  show  the 
effects  of  sloping  terrain.  This  effect  can  be  approximated 
by  drawing  a  sloping  line  through  the  pattern  footprint 
graphs  in  the  completed  guide. 

The  schedule  of  PATSIM  runs  should  include  runs  for 
each  different  type  of  compartment/drop  at  100-foot  inter- 
vals between  100  and  500  feet.  Where  identical  compart- 
ments are  released  sequentially  (trail  drop),  the  PATADD 
option  should  be  exercised.  Because  the  PATADD  option 
will  add  only  identical  compartments,  a  trail  drop  where 
nonidentical  compartments  are  released  requires  the 
application  of  the  separate  PATSUM  Program.  This  al- 
lows the  input  of  two  or  more  patterns  for  nonidentical 
compartments  obtained  from  previous  PATSIM  runs. 

Appendix  B  contains  the  data  from  the  computer  simu- 
lation which  consists  of: 

1.  Detailed  pattern  plots  used  in  generating  footprints 
and  for  entry  into  the  PATSUM  routine  for  summation  of 
dissimilar  trail  drops. 

2.  Pattern  area,  length,  and  volumes.  The  remainder  is 
for  inspection  only. 

3.  Maximum  line  length/tank  sequence  tables  from 
PATADD.  These  are  summarized  and  plotted  to  yield 
best  strategy  charts  and  are  presented  in  tabular  form  in 
the  guides.  (PATSUM  is  used  for  dissimilar  tanks.) 
Complete  instructions  for  PATSIM  use  may  be  found  in 
"Development  of  User  Guidelines  for  Selected  Retardant 
Aircraft"  (Swanson  and  others  1977). 

Detailed  Guideline  Preparation 

A  sample  guideline  is  included  as  appendix  C.  At  the 
beginning  of  the  guide  is  a  one-page  summary  of  the 
aircraft/tank  characteristics,  together  with  an  outline 


drawing  of  the  aircraft,  a  detailed  isometric  tank  drawing, 
a  summary  table  of  tank  system  characteristics,  and 
a  narrative  summary  of  the  aircraft/tank  system 
performance. 

The  narrative  summary  contains  an  aircraft  description 
and  tank  options,  characteristics  of  the  tank  and  how  they 
control  the  pattern,  a  summary  of  the  best  strategy 
charts,  noting  the  most  useful  drop  methods,  ground 
safety  precautions,  and  any  unique  characteristics. 

The  characteristics  chart  contains  the  total  tank  capac- 
ity, the  increment  contained  in  each  separately  releasable 
compartment,  evacuation  time  (total  time  to  empty),  the 
average  flow  rate  (capacity  of  tank  released/evacuation 
time),  and  peak  flow  rate.  Peak  flow  rate  is  calculated 
from  the  flow  data  used  as  input  to  the  simulation  as: 

peak  flow  rate  =  maximum  value  of  ^  from 
measured  flow  data 

when  DQ  is  the  volume  in  gallons  released  during  time 
interval  DT. 

Pattern  Coverage  Characteristics 
(Footprints) 

Pattern  footprints  display  pattern  performance  as  a 
continuous  function  of  drop  height  for  each  combination 
of  simultaneously  released  tanks,  door  option  (if  appli- 
cable), and  each  retardant  type.  The  procedure  explained 
below  is  illustrated  in  figure  6.  Pattern  length  and  initial 
down-range  impact  for  each  coverage  level  are  extracted 
from  the  central  pattern  section  (columns  14  or  15)  of  the 
detailed  pattern  plots.  All  coverage  levels  are  Unearly 
interpolated  to  the  nearest  foot  of  pattern  length.  At  zero 
drop  height,  a  length  representing  aircraft  travel  during 
the  time  of  tank  evacuation  is  calculated  to  be  used  as  the 
end  point  for  the  plotted  footprints  (speed  in  feet  per  sec- 
ond evacuation  time).  These  points  are  plotted  for  each 
coverage  level  1  through  5  and  a  smooth  curve  drawn  to 
form  the  footprints.  Footprint  half-width  for  each  cover- 
age level/drop  height  is  found  by  averaging  the  width 
from  each  row  of  the  detailed  pattern  plots.  These  half- 
widths  are  plotted  and  smoothed  similarly  to  the  pattern 
footprints.  The  values  for  length,  downrange  impact 
point,  and  width  are  calculated  and  tabulated  using  the 
summary  program  LDRWIDTH;  the  footprint  is  smoothed 
using  the  SMOOTHFOOT  program,  and  final  plots  are 
made  using  FINALFOOT. 

Scaled  silhouettes  of  the  side  and  head-on  views  of  the 
aircraft  are  placed  above  the  footprint  and  half-width 
graphs  to  indicate  the  relationships  of  the  graphs  to  air- 
craft flight  direction.  This  process  is  accomplished  using 
several  BASIC  programs. 

Also  included  on  the  pages  with  pattern  footprints  are 
a  brief  set  of  instructions  for  use  of  pattern  footprints  and 
a  table  of  recommended  coverage  levels  for  different  fuel 
types  (Deeming  1972;  Rothermel  and  Philpot  1974). 
LEDRWIDTH,  SMOOTHFOOT,  and  FINALFOOT  are 
listed  in  appendix  D. 
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Best  Strategy  Charts 

The  section  of  the  guide  containing  the  best  strategy 
charts  is  preceded  by  a  page  of  instructions  for  their  use. 
Best  strategy  charts  are  graphic  summaries  of  selected 
pattern  information  generated  by  using  the  PATADD 
option  to  PATSIM  (fig.  7).  These  summaries  are  also 
tabulated  and  presented  following  the  best  strategy 
charts.  Best  strategy  charts  are  generated  for  coverage 
levels  1  through  4  for  waterlike  and  gum-thickened 
retardants.  For  a  selected  coverage  level  and  retardant 
type,  the  length  of  line  at  drop  heights  of  100  to  500  feet 
for  each  possible  load  fi-action  are  plotted. 

For  a  tank  with  four  identical  compartments,  there  are 
six  possible  combinations: 

1.  Single-compartment  release  (ITxlR) 

2.  Two  compartments  released  simultaneously  (2TxlR) 

3.  Four  compartments  released  simultaneously  (4TxlR) 

4.  Two  compartments  released  sequentially  (lTx2R) 

5.  Four  compartments  released  sequentially  (lTx4R) 

6.  Two  two-compartment  sequential  releases  (2Tx2R) 

Each  of  these  combinations  is  regressed,  using  a  stan- 
dard multiform  regression  program,  and  roughly  plotted. 
From  this  rough  plot,  the  optimum  drop  for  each  load  size 
can  be  selected  for  the  final  best  trail  strategy  chart.  For 
example,  if  the  four-compartment  in  (lTx4R)  produces  a 


longer  line  length  than  a  four-compartment  salvo  (4TxlR), 
the  four-compartment  salvo  would  not  be  plotted  on  the 
final  chart  (fig.  7)  Gength  and  width  dimensions  for  the 
salvo  could  be  found  on  the  appropriate  pattern  footprint). 
Each  area  is  identified  with  the  proper  label  and  a  nomi- 
nal time  interval  between  releases  in  the  case  of  sequen- 
tial releases.  The  proper  time  interval  changes  wdth  drop 
height  and  airspeed  and  is  reported  to  0.1  second.  The 
value  of  the  interval  at  200  feet  for  the  drop  type  in  ques- 
tion is  chosen  as  being  most  representative.  A  minimum 
drop  height  safety  limit  of  150  feet  above  ground  deter- 
mined by  Forest  Service  policy  is  shown  on  the  bottom  of 
the  chart. 

The  accuracy  limit  curve  (Swanson  and  others  1975) 
is  based  on  a  Honeywell-derived  probability  of  hit  for  mili- 
tary ordnance  and  is  common  to  the  best  strategy  charts 
for  all  aircraft.  It  corresponds  to  a  curve  defined  by  set- 
ting the  probability  of  hit  at  0.6  for  no  wind  and  level 
flight.  The  curve  is  drop-height-dependent,  indicating  the 
length  of  line  below  which  a  hit  is  improbable.  Because 
the  curve  is  independent  of  the  pattern  width,  which  is 
relatively  constant  for  full -door  drops,  a  separate  curve  is 
generated  for  aircraft  equipped  with  a  trail  door  option, 
which  produces  a  narrower  retardant  pattern.  This  arbi- 
trary change  simply  moves  the  300-foot  length  limit  to  the 
500-foot  drop  height.  Best  strategy  charts  are  produced 
using  the  BASIC  programs:  BSPLOT,  BESTSTRAT.DAT. 
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Included  on  the  pages  containing  best  strategy  charts 
are  a  key  to  the  notation  on  the  charts,  the  table  of  recom- 
mended coverage  levels,  and  a  table  of  ground  safety 
notes.  The  ground  safety  table  shows  for  each  load  type 
the  height  below  which  the  retardant  mass  has  lost  for- 
ward momentum  and  broken  up  into  droplets,  falling  as 
a  heavy  rain.  The  equation  for  these  heights  was  derived 
from  films  of  test  drops  and  is  a  function  of  load  size  and 
peak  flow  rate:  Safe  drop  height  =  (100.73  +  0.0112  x  load 
size  +  0.0202  x  peak  flow  rate)  +  50.  The  arbitrary  50  feet 
are  added  as  a  safety  margin. 

The  remainder  of  the  guide  consists  of  detailed  tables 
of  the  aircraft's  performance  transcribed  and  rearranged 
from  the  PATADD  subroutine  and  technical  data  tables. 
The  technical  data  tables  contain  some  physical  dimen- 
sions of  the  tank  system  and  flow  data  derived  from  the 
static  test.  These  tables  potentially  enable  field  personnel 
to  identify  tank  changes  that  might  affect  the  guide  valid- 
ity or  to  judge  the  similarity  to  tank  systems  for  which 
guides  are  not  available.  Also  included  are  the  date  and 
full  identification  of  the  tanker  static-tested  and  any  other 
identical  airtankers. 

SUMMARY 

To  date  airtanker  performance  guidelines  have  been 
published  for  most  of  the  airtankers  currently  being  used. 
These  are  listed  in  table  2.  An  instruction  manual,  "Air- 
tanker  Performance  Guides:  General  Instruction  Manual" 
(Swanson  and  others  1976)  contains  detailed  instructions 
on  how  to  read  and  use  airtanker  performance  guides  that 
contains  general  information,  retardant  coverage  require- 
ments, pattern  footprints,  best  strategy  charts  and  de- 
tailed line  length  tables. 

A  need  was  identified  for  a  small,  simple,  inexpensive 
reference  that  could  be  used  in  real-time  to  identify  pri- 
mary airtanker  performance.  To  fill  this  need  airtanker 
performance  "slide  charts"  (retardant  coverage  computers) 
have  been  developed  along  with  an  instruction  manual  for 
their  use  (George  1981).  The  retardant  coverage  comput- 
ers are  a  slide-rule-type  device  that  can  quickly  show  the 
user  the  length  of  retardant  line  and  the  intertank  release 
interval  used  for  any  coverage  level  (1.0  -  5.0  +  gpc),  drop 
height  from  100  to  500  feet,  and  drop  type.  Retardant 
coverage  computers  have  been  produced  for  each  air- 
tanker covered  by  performance  guides. 

As  new  types  of  tank  and  gating  systems,  whether  new 
aircraft  types  or  systems  modified  to  improve  flexibility 
and  performance,  come  into  use,  new  or  updated  guides 
and  slide  charts  should  be  produced  in  order  to  keep  this 
valuable  information  current  and  useful. 

With  the  advent  of  computer-based  intervalometers  to 
select  drop  types  and  intervals  in  the  newer  airtankers, 
a  potential  exists  to  include  ground  pattern  performance 
data  in  the  intervalometer.  The  desired  coverage  level 
could  be  selected  and  the  computer  could  set,  or  indicate 
to  the  air  crew,  the  most  effective  drop  type  and  interval. 


Table  2 — List  ot  airtanker  performance  guidelines 


Aero  Union  tanker  B-17  (Request  by  title) 

Gilbertson  tanker  DC-6B  (Request  by  title) 

CDF/Hemet  Valley  tanked  S2F  PG-1 

CDF/Aero  Union  tanked  S2F  PG-2 

Raico  tanked  PV-2  PG-3 
Reader  Tank/Lynch  STOL  B-26 

(tanker  58)  PG-4 

Canadair  CL-21 5  PG-5 

Evergreen  Rosenbalm  tanked  B-17  PG-6 

Globe  tanked  B-1 7  PG-7 

Black  Hills  tanked  B-1 7  PG-8 

Aero  Union  tanked  DC-^/6/7  ( 1 ,800  gal)  PG-9 

Aero  Union  tanked  DC-4/6/7  (2,000  gal)  PG-10 

Aero  Union  tanked  DC-4/6/7  (2,200  gal)  PG-1 1 

Aero  Union  tanked  DC-4/6/7  (2,400  gal)  PG-1 2 

Aero  Union  tanked  DC-4/6/7  (2,600  gal)  PG-1 3 

Aero  Union  tanked  DC-4/6/7  (2,800  gal)  PG-1 4 

Aero  Union  tanked  DC-4/6/7  (3,000  gal)  PG-1 5 

W.A.I. G.  tanked  DC-4  PG-1 6 
Hemet  Valley/Aero  Union  tanked 

C119G-3E  PG-1 7 

Hawkins  &  Powers  tanked  CI  19G-3E  PG-18 

Hawkins  &  Powers  tanked  PB4Y-2  PG-1 9 

Black  Hills/Rosenbalm  tanked  P2V-5  PG-20 

Black  Hills/Rosenbalm  tanked  P2V-7  PG-21 

Central  Air  Services  tanked  DC-7  PG-22 

Transwest  tanked  DC-7  PG-23 
SIS-Q/Rosenbalm  tanked  DC-6/7 

(3,000  gal)  PG-24 

C-130  MAFFS  PG-25 

Black  Hills  tanked  P2V-5  PG-26 

Lynch  tanked  STOL  B-26  (tanker  01 )  PG-27 

Evergreen  P2V-5  PG-29 
SIS-Q/Rosenbalm  tanked  DC-6/7 

(2,000  and  2,450  gal)  PG-30 

TBM,  Inc.  F7F  PG-31 

Lynch  STOL  B-26  (tanker  57)  PG-32 

Central  Air  Services  DC-4  (2,000  gal)  PG-33 

Conair  DC-6B  PG-34 

TBM,  Inc.  C-123  PG-35 
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APPENDIX  A:  PATSIM  LISTING 

110  REM  PROGRAM  PATSIM ( INPUT , OUTPUT , TAPCONV=INPUT , TAPTERED=OUTPUT) 

120  REM  ('■')   COMMON/MIX/DST ( 1700 )  ,SUMD, RAT 

121  DIM  Dst[1700] 

130  REM  CO   COMMON  /ADD/  NPTS , VALIN ( 56 ) , ITC ( 8 ) 

131  DIM  Valin[56] , Itc[8] 

140  DIM  Timl [500] ,Vex[500] ,Qo[2500] ,Qdoto[2500] ,Qdotx[500] 

150  DIM  Zo[2500] ,Veo[2500] ,Vzo[2500] 

160  DIM  Dss[200] 

170  REM  (*)  REAL  K1,K2 

180  DIM  A [90] 

185  DIM  Itc[8] 

200  DIM  Ttim[100] ,Thim[100] 

205  Ctrl$  =  CHR$(12) 

210  IMAGE  8(4D.5D) 

220  IMAGE  8(5E) 

221  IMAGE  3(5E,4X) 

230  IMA  1A,20X,"S  TANDARD     PATTERN     SIMULATION  ",li,30X,"R  ETARDANT       PATTER  N",2L 
240  IMAGE  5X,5D,3X,6E 
250  IMAGE  20 (4A) 
260  IMAGE  1(1A) 
270  IMAGE  10X,80A 

280  IMAGE  2L,10X,"VELOCITY(KNOTS)=",8D,1L,10X,"ALTITUDE(FEET)=  " , 8D 

281  INPUT  PROMPT  "What  data  file?  ":FS 

282  OPEN  #1 :F$, "r" 

283  ON  EOF(l)   GOTO  2550    !    End  of  program 

284  OPEN  #7 : "prn  " , "w" 


300 

A  =  0 

! 

Inital ize 

310 

A[32] 

0    !  Bflag 

320 

A[18] 

1    !  Cflag 

330 

A[37] 

1    !  Alf 

340 

A[40] 

3.97    !  Kl 

350 

A[41] 

12    !  K2 

360 

A[A2] 

0.002    !  Dc 

370 

A[43] 

10    !  Tr 

380 

A[44] 

0.02    !  Tx 

390 

Irdqv 

0 

400 

A[45] 

1    !  Rflag 

410  REM  (")  CONTINUE 

411  DELETE  Tp , Leon , Conv , Sing , Lcmax , Dela 

412  DIM  Dst[1700] ,Qo[2500] ,Qdoto[2500]    !    , Qdotx [ 500 ] , TIMl [ 500 ] , VEX [ 500 ] 

413  DIM  Zo[2500] ,Veo[2500] ,Vzo[2500] 
420  PRINT  #7   USING  230:CtrI$ 

440  INPUT  #1:A$   !   Read  from  data  file 

441  IF  A$=""  THEN  2550 

450  AtstS  =  SEG$(A$,1,1)    |    IF  Atst$="*"  THEN  1350    !    Is  this  a  remark  line? 
460  IF  A$="   "  THEN  1370   !   Is  this  the  end  of  a  data  set? 
465  Al  =  VAL(AS)    !  Al  is  the  variable  code 
470  IF  AKO  THEN  1370 

480  INPUT  #1:A2    !   A2  is  the  value  of  the  variable 

510  PRINT  #7  USING  240:A1,A2 

530  J  =  Al    !   Renaming  the  variable  label 

540  IF  J<=0  THEN  1300 

550  A[J]  =  A2   !  Renaming  the  variable  value 

555  GOSUB  10000    !    Subroutine  to  change  the  names   in  the  A  array 
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556  IF  J<>5  THEN  560 

557  INPUT  #1:M,NS,C$    !   drop#,   comment,    file  for  storage 

560  IF  J0  53  OR  A[53]<=0  THEN  580    !   Not  Add 

561  FOR  K  =  1  TO  8    |    INPUT  #l:Itc[K]    |    NEXT  K   !   Which  patterns  are  computed? 
570  IMAGE  8 (ID) 

580  IF  JO30  THEN  620    !   Not  Ttab   (tabular  door  angle  vs  time) 

590  IF  A[J]=0  THEN  620    !   Ttab=0  deselected   (Tabular  door  angle  vs  time) 

600  PRINT  #7  USING  601: 

601  IMAGE  11, "DOOR  ANGLE  VS  TIME  DATA" , 7X , "TIME-SEC" , 6X , "ANGIE-DEG" 
610  INPUT  #l:Ttim[N] ,Thtim[N]    !   Door  angle  vs  time  data 

620  REM  ('■')  CONTINUE 

630  IF  J<>24  THEN  1300   !   Not  Rdqv  (Quantity  discharged  vs  time) 

640   IF  A[J]=0  THEN  1300    !   Rdqv=0  deselected   (Quantity  discharged  vs  time) 

650   IF  A[J]>0  THEN  780    !   Rdqv  selected   (Quantity  discharged  vs  time) 

660  REM  INPUT  Q  AND  THETA  VS.    TIME  AND  USE     TABLE  AS  IS 

670  Rdqv  =  -Rdqv 

680  Ant  =  Rdqv 

690  Nt  =  Ant   !   number  of  points  read  from  file 

700  FOR  N  =  1  TO  Nt    |    INPUT  #1 : Timl [N ] , Qdoto [ N ] , Veo [N ]    |    NEXT  N   !   Read  data 

710  IMA  "  QUANTITY  AND  DOOR  ANGLE  VS  TIME  DATA" , 11 , 2X, "TIME-SEC" , 8X, "Q-GAL" , 8X, "ANGLE-DEG' 
720  PRINT  #7  USING  710: 
730  FOR  N  =  1  TO  Nt 

7A0       PRI  #7  USI  220:Timl[N] ,Qdoto[N] ,Veo[N]    !   Quantity  disc,    and  theta  vs  time 
7A1  NEXT  N 
750  GOTO  1280 

760  REM  INPUT  Q  AND  THETA  VS.    TIME,    INTERPOLATE  VALUES  EVERY  TX  SECONDS, 
770  REM  AND  EXTRAPOLATE  THE  TABLE  TO  GALI  NUMBER  OF  GALLONS 
780  Nti  =  Rdqv 

790  FOR  N  =  1  TO  Nti    |    INPUT  #1 : Qo [ N ] , 2o [ N ] , Vzo [ N]    |   NEXT  N   !   Read  data 
800  PRINT  #7  USING  710: 
810  FOR  N  =  1  TO  Nti 

820  PRINT  #7  USING  221 : Qo [ N] , Zo [N ] , Vzo ( N]    !   Quantity  disc,    and  theta  vs  time 

821  NEXT  N 

830  Timl[l]   =  0   !   Beginning  the  extrapolation  between  Time , Quantity  points 
SAO  Qdoto[l]   =  0 
850  Veo[l]   =  0 
860  1=1 
870  I  =  I+l 

880  Timl[I]  =  Timl[I-l]+Tx 
890  FOR  N  =  1  TO  Nti 

900       ON  SGN(Timl[I]-Qo[N] )+2  GOTO  1060,990,910 

910  REM  ('')   CONTINUE   !   Case:   positive  values 

911  NEXT  N 

920  Timsav  =  Timl [I] 
930  Timl[I]  =  Qo[Nti] 
9A0  Qdoto[I]   =  Zo(Nti] 
950  Veo[I]  =  Vzo[Nti] 
960  Ant  =  I 
970  Nt  =  I 
980  GOTO  1090 

990  Qdoto[I]   =  ZofN]    !    Case:  zero 
1000  Veo[I]   =  0.5''(Vzo[N]+Vzo[N-l]  ) 
1010   IF  NoNti  THEN  870 
1020  Timsav  =  Timl[I]+Tx 
1030  Ant  =  I 
lOAO  Nt  =  I 
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lUiO  GOTO  loyo 

1059  !   Case:   negative  values 

1060  Qdoto[I]   =  (Timl[I]-Qo[N-l])/(Qo[N]-Qo[N-l])'V(Zo[N]-Zo[N-l])+Zo[N-l] 
1070  Veo[I]   =  0.5"(Vzo[N]+Vzo[N-l] ) 

1080  GOTO  870 

1090  IF  Qdoto[I]>=Gali-l.E-06  THEN  1280    !    extrapolation  finished? 

1100  Tfin  =  (Gali-Zo[Nti-l])/(Zo[Nti]-Zo[Nti-l])"(Qo[Nti]-Qo[Nti-l])+QorNti-l] 
1110  I  =  I+l 

1120  Timl[I]   =  Timsav 

1130  GOTO  1160 

1140  I  =  I+l 

1150  Timl[I]   =  Timl[I-l]+Tx 

1160  IF  Timl[I]>Tfin  THEN  1230 

1170  Qdoto[I]   =  (Gali-Zo[Nti]  )^''(  (Timl  [I] -Qo[Nti]  )  /  (Tf  in-Qo[Nti]  )  )+Zo[Nti] 

1180  Veo[I]   =  Vzo[Nti] 

1190  IF  Timl[I]<Tfin  THEN  1140 

1200  Ant  =  I 

1210  Nt  =  I 

1220  GOTO  1280 

1230  Timl[I]  =  Tfin 

1240  Qdoto[I]   =  Gali 

1250  Veo[I]   =  Vzo[Nti] 

1260  Ant  =  I 

1270  Nt  =  I 

1280  REM  CONTINUE   !   End  of  extrapolating  between  points 

1290  Irdqv  =  1    !    set  flag  quantity  dischaged  data  read  from  file 

1300  REM  (")  CONTINUE 

1310  GOTO  440 

1320  STOP 

1350  PRINT  #7  USING  270 :A$   !   Print  comment  line 

1360  GOTO  440 

1370  REM  (*)  CONTINUE 

1380  PRINT  #7  USING  280 : Vac i , Altl    !   End  of  reading  data  set 

1390  Kl  =  3.97 

1400  IF  Rflag=2  THEN  Kl  =  4 . 4   !   flag  Rflag=2  for  water-like  retardants 

1410  Ve  =  0.0 

1420  FOR  J  =  1  TO  1700    !    initalize  Dst 
1430       Dst[J]   =  0.0 

1440  REM  ('■')  CONTINUE 

1441  NEXT  J 

1450  FOR  I  =  1  TO  100    !    initalize  Dss 
1460       Dss[I]   =  0.0 

1470  REM   ('■')  CONTINUE 

1471  NEXT  I 
1510  Jr  =  Tr 
1520  Jr2  =  Jr/2 

1530  Vac  =  Vaci''1.689   !  knots  to  feet  per  second 

1540  Gal  =  Gali''Cflag  !   amount  of  retardant  released 

1550  K  =  1 

1560  Dt  =  0.02 

1580  G  =  32.17   !   acceleration  due  to  gravity 

1590  G2  =  2.0'''G 

1600  IF  Irdqv<l  THEN  1680   !   flag  Irdqv=l  then  data  read  from  file 

1610  Qdotx[l]  =  0 

1620  Vex[l]  =  0 

1630  FOR  N  =  2  TO  Nt   !   for  constant  ejection  veIocity=16  fps 

1640      Qdotx[N]  =  (Qdoto(N]-QdototN-l])*0.13368/(Timl[N]-Timl[N-l]) 
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1650       Vex[N]   =  16 

1660  REM  i*)  CONTINUE 

1661  NEXT  N 
1670  Irdqv  =  0 

1680  REM  (^)  CONTINUE 

1690   IF  A[54)<>0  THEN  410    !    flag  A(54)=l   stop  after  flow  value   is  determined 

1700  Dt  =  0.02 

1710  Z  =  0 

1720  DelE  =  4.0 

1730  Jk  =  -Jr2+1 

1740  Xout  =  1 

1750  Qsum  =  0 

1760   Zt  =  0 

1770  Dqout  =  Qdotx[l] 

1780  FOR  N  =  2  TO  Nt 

17  90       Jk  =  Jk+Jr 

1800       Dt  =  Timl[N]-Timl[N-l] 

1810       Qs  =  Qdotx[N]"Dt 

1820       Qsum  =  Qs+Qsum 

1830       IF  BflagOl  THEN  1850   !   Bflag=l  then  print  bulk  increment  data 
1840       PRINT  #7  USING  220 : Timl [N] , Vex [N] , Qdotx [N] , Qs , Qsum 
1850      REM  CONTINUE 

1860  REM  (")  CONTINUE 

1861  NEXT  N 

1870  IF  Bflag<>2.0  THEN  1890 

1880  GOTO  410 

1890  REM  ('•')  CONTINUE 

1900  Ant  =  Nt 

1910  Xout  =1.0 

1920  Xlast  =1.0 

1930  Qqq  =  0.0 

1940  FOR  Km  =  2  TO  Nt 

1950       Dt  =  Timl [Km] -Timl [Km-1]    !  time 
1960       Qsum  =  0.0 

1970      Ve  =  Vex[Km]    !   ejection  velocity 
1980       Qdot  =  Qdotx  [Km] 

1990       Q  =  Qdof'Df'Cf lag    !   quantity  discharged 
2000       IF  Q<=0  THEN  2310 

2010       Qqq  =  Qqq+Q   !   cummulative  quantity  discharged 

2020       Xout  =  1 

2030       Zt  =  0 

2040       Xx  =  0 

2050       Tstar  =  0 

2060       Tmt  =  Timl [Km] 

2070      Xlast  =  Tmf'Vac 

2080       FOR  Jj   =  2  TO  100 

2090  Zt  =  Zt+Delz 

2100  V  =  SQR(Ve^2+G2'''Zt) 

2110  Tstar  =  Deiz/V+Tstar 

2120  Ttt  =  Tstar+Tmt 

2130  Dqd  =  EXP (Kl ''Tstar )/ (Vac'>Q ) •*K2 

2140  Xx  =  Delz''Vac/V+Xx 

2150  Xout  =  Xx+Tmt"Vac 

2160  Xr  =  Xout 

2170  Nil  =  Xlast 

2180  NI2  =  Xr 

2190  Diff  =  Xr-Xlast+1 
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2200  Dqout  =  DeLz/V'Dqd 

2210  REM  DENS=DQD'*DELZ/ (V'DIFF) 

2220  IF  Dqout>Q-Qsum  THEN  Dqout  =  Q-Qsum 

2230  Dens  =  Dqout/Diff 

2240  FOR  N2  =  INT(Nll)   TO  INT(N12)    !   calculating  dst 

2250  Dst[N2]  =  Dst[N2]+Dens 

2251  NEXT  N2 

2260  Qsum  =  Qsum+Dqout 

2270  Xlast  =  Xout 

2280  IF  Qsum>=Q  THEN  EXIT  TO  2300 

2290  REM  C*)  CONTINUE 

2291  NEXT  Jj 

2300      REM  C')  CONTINUE 

2310  REM  ('■')  CONTINUE 

2311  NEXT  Km 
2320  Jm  =  0 
2330  Sumd  =  0 

23A0  FOR  N3  =  1  TO  1700 
2350       Sumd  =  Sumd+Dst[N3] 

2360  REM  CONTINUE 

2361  NEXT  N3 
2370  Jk  =  0 

2380  Rat  =  Gal /Sumd 

2390  FOR  Ij  =  1  TO  1700  STEP  25 

2A00       Jk  =  Jk+1 

2410      FOR  Jj  =  1  TO  25 

2420  Jm  =  Ij+Jj-1 

2430  Dss[Jk]   =  Dss[Jk]+Dst[Jm] 

24  40  REM  ('■')  CONTINUE 

2441       NEXT  Jj 

2450       Dss[Jk]  =  Dss[Jk]*Rat 

24  60      REM  ('■')  CONTINUE 

2461  NEXT  Ij 

24  70  Sumk  =  0 

2480  FOR  I  =  1  TO  68 

2490       Sumk  =  Dss[I]+Sumk 

2500  REM  ('■')  CONTINUE 

2501  NEXT  I 

2510  REM  CALL  MIXUN (GAL , VAGI , ALTl , RFLAG , FPAT , FRAN , FE ) 

2511  PRINT  "subroutine  MIXUN" 

2512  GOSUB  2560 

2520  IF  Add<=0  THEN  410    !   don't  calculate  line  length  table 

2530  REM  CALL  PATADD   (VAC I) 

2531  PRINT  "subroutine  PATADD" 

2532  GOSUB  6550 
2540  GOTO  410 

2550  CLOSE 

2551  REM 

2552  END 

2560  REM  SUBROUTINE  MIXUN ( GAL , VAC , ALTl , RFLAG , FPAT , FRAN , FE) 

2570  REM  (")   COMMON/MIX/DST(1700) , SUMD, RAT 

2580  REM  ('■')  COMMON  /ADD/  NPTS  ,  VALIN (  56)  ,  ITC (8) 

2590  DIM  Ud[28,56] ,Nd[28,56] ,Cd[28,56] 

2600  DIM  Mi[60] ,Mo[60] ,M1 [60] ,M2[60] 

2610  DIM  Cmi[60] ,Cmo[60] ,Cml[60] ,Cm2[60] 

2620  DIM  P[15] 

2630  DIM  Mii[60] ,Cmii[60] 


15 


2640  REM  (■")  REAL  ND  ,  LAMDA , MI ,  MO , Ml , M2  ,  CMII 

2650  Im  =  28 

2660  Jm  =  56 

2670  Delx  =  15 

2680  Dely  =  30 

2690  Alt2  =  Altl"Altl 

2700  IF  Rflag>l  THEN  2850 

2710  REM  C-R  PARAMETERS  FOR  PHOSCHEK  RETARDANT 

2720  IF  Vaci<=130  THEN  2750 

2730  Semx  =  AO  .  0+0  .  104048"'Altl-6  .  90477E-05'''Alt2 
2740  GOTO  2760 

2750  Semx  =  40  .  0+0  .  255714"Altl-0  .  0OO185714'''ALt2 

2760  REM  CONTINUE 

2770  Lamda  =  0 . 15+5 . E-05*Altl 

2780  IF  VacL<=130  THEN  2810 

2790  Sigma  =  8  .  0+0  .  0485715"'Altl-2  .  85715E-05'''Alt2 
2800  GOTO  2820 

2810  Sigma  =  8  .  0+0  .  0659048'''Altl-4  . 19048E-05'<Alt2 

2820  REM  (")  CONTINUE 

2830  Pctr  =  85.18-0.00535'''Altl 

2840  GOTO  2960 

2850  REM  (")  CONTINUE 

2860  REM  C-R  PARAMETERS  FOR  WATER  LIKE  RETARDANTS 

2870  IF  Altl>482  THEN  2910 

2880  Semx  =  27 . 5583+0 . 222214"Altl 

2890  Sigma  =  3 . 11609+0 . 0746535"Altl 

2900  GOTO  2930 

2910  Semx  =  135 

2920  Sigma  =  36  .  0427+0  .  00821018''Altl 

2930  REM   (>•')  CONTINUE 

2940  Lamda  =  0  . 15+5  .  E-05'"Altl 

2950  Pctr  =  75.18-0.02487"Altl 

2960  REM  (>■')  CONTINUE 

2970  REM 

2980  Areaf  =  100/ (Delx"Dely) 

2990  REM  FORM  MARGINAL  RANGE  WITH  CELL  SIZE  =  DELY 

3000  FOR  I  =  1  TO  Jm 

3010  Mi[I]   =  0 

3011  NEXT  I 
3020  Jk  =  0 

3030  FOR  Ij  =  1  TO  1680  STEP  30 

3040       Jk  =  Jk+1 

3050       FOR  Jj  =  1  TO  30 

3060  Jn  =  Ij+Jj-1 

3070  Mi[Jk]   =  Ml[ Jk]+Dst[Jn] 

3080  REM  CONTINUE 

3081  NEXT  Jj 

3090       Mi[Jk]   =  Mi[Jk]*Rat 

3100  REM  ('•')  CONTINUE 

3101  NEXT  Ij 

3110  FOR  I  =  1  TO  Jm 

3120  Mii[I]   =  Mi[I)"(Pctr/100.0) 

3121  NEXT  I 
3130  REM 

3140  REM  FORM  ELLITICAL  UNIFORM  GROUND  DISTRIBUTION  AND  MARGINAL 

3150  Xm  =  Im/2  .  0'''Delx 
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3160  FOR  I  =  1  TO  Im 
3170      FOR  J  =  1  TO  Jm 
3180  Ud[I,J]  =  0 

3190  Nd[I,J]  =  0 

3200  Cd[I,J]  =  0 

3201  NEXT  J 

3202  NEXT  I 
3210  Jml  =  1 

3220  FOR  I  =  1  TO  Jm 

3230       IF  Mitl]>0  THEN  EXIT  TO  3260 

3240       Jml  =  Jml+1 

3250  REM  ('■')  CONTINUE 

3251  NEXT  I 

3260  REM  (.■")  CONTINUE 

3270  FOR  I  =  Jml  TO  Jm 

3280       IF  Mi[I]<=0  THEN  EXIT  TO  3310 

3290       Jm2  =  I 

3300  REM  CONTINUE 

3301  NEXT  I 

3310  REM  ('■')  CONTINUE 
3320  A2  =  Semx'2 

3330  B2  =  ((Jm2-Jml+l)'-'Dely/2.0)^2 

3340  Ym  =   ( Jm2+Jml - 1 ) '^Dely  /  2  .  0 

3350  Nc  =  0 

3360  X  =  -Delx/2 

3370  FOR  I  =  1  TO  Im 

3380       X  =  X+Delx 

3390       Y  =  -Dely/2 

3400       FOR  J  =  1  TO  Jm 

3410  Y  =  Y+Dely 

3420  Ee  =  (X-Xm) ^2/A2+(Y-Ym) "2/B2 

3430  IF  Ee-1>0  THEN  3460 

3440  Ud[I,J]  =  1 

3450  Nc  =  Nc+1 

3460  REM  CONTINUE 

3461  NEXT  J 

3470  REM  (•■')  CONTINUE 

3471  NEXT  I 

3480  Fe  =  1.0/INT(Nc) 

3490  Ae  =  Lamda^Gal^CPctr/ 100 . 0) 

3500  Ca  =  Fe*Ae 

3510  FOR  I  =  1  TO  Im 

3520       FOR  J  =  1  TO  Jm 

3530  IF  Ud[I,J]=l  THEN  Ud[I,J]  =  Ca 

3540  REM  (*)  CONTINUE 

3541  NEXT  J 

3542  NEXT  I 

3550  FOR  J  =  1  TO  Jm 

3560       M1[J]   =  0 

3570       FOR  I  =  1  TO  Im 

3580  M1[J]   =  Ml[J]+Ud[I, J] 

3590  REM  CO  CONTINUE 

3591  NEXT  I 

3600  REM   (")  CONTINUE 

3601  NEXT  J 
3610  REM 

3620  REM  SUBTRACT  UNIFORM  MARGINAL  FROM  MACPHERSON  MARGINAL 
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3630  FOR  J  =  1  TO  Jm 

3640       Temp  =  MLi(J)-Ml[J] 

3650      ON  SGN(Temp)+2  GOTO  3660,3680,3680 
3660      M2[J]   =  0 
3670       GOTO  3690 
3680       M2[J]   =  Temp 

3690  REM  ('■')  CONTINUE 

3691  NEXT  J 
3700  REM 

3710  REM  DEVELOP  NORMAL  DISTRIBUTION  CONDITIONALS 

3720  PI  =  0.5 

3730  Im2  =  INT(Im/2) 

3740  X  =  0 

3750  FOR  I  =  1  TO  Im2 

3760      X  =  X+Delx 

3770       Xos  =  X/Sigma 

3778       GOTO  4800 

3780       Pn  =  Cdfn 

3790       IF  Xos>4  THEN  Pn  =  1 

3800       P[I]   =  Pn-Pl 

3810       PI  =  Pn 

3820  REM  (>■')  CONTINUE 

3821  NEXT  I 

3830  FOR  J  =  1  TO  Jm 

3840       IF  M2[J]<=0  THEN  3920 

3850       FOR  I  =  1  TO  Im2 

3860  II  =  Im2-I+1 

3870  12  =  I+Im2 

3880  Temp  =  P[I]'''M2[J] 

3890  Nd[Il , J]   =  Temp 

3900  Nd[I2,J]   =  Temp 

3910  REM   (■■')  CONTINUE 

3911  NEXT  I 

3920      Mo[J]   =  M1[J]+M2[J] 

3930  REM  ('•')  CONTINUE 

3931  NEXT  J 

3940  REM  COMBINE  UNIFORM  AND  NORMAL  DISTRIBUTIONS 

3950  Suml  =  0 

3960  Sum2  =  0 

3970  Sum3  =  0 

3980  Sum4  =  0 

3990  Sum9  =  0 

4000  FOR  J  =  1  TO  Jm 

4010       Suml  =  Suml+Mi[J] 

4020       Sum9  =  Sum9+Mii[J] 

4030       Sum2  =  Sum2+Ml[J] 

4040       Sum3  =  Sum3+M2[J] 

4050  Sum4  =  Sum4+Mo[J] 

4051  NEXT  J 
4060  Sum5  =  0 
4070  Suni6  =  0 
4080  Suj7i7  =  0 
4090  Sum8  =  0 
4100  SumlO  =  0 

4110  FOR  J  =  1  TO  Jm 
4120       Sum5  =  Sum5+Mi[J] 
4130       SumlO  =  SumlO+Mii[J] 


4140       Sum6  =  Sum6+Ml[J] 

4150       Sum?  =  Sum7+M2[J] 

/.1 60       Sum8  =  Sijm8+Mo[J] 

4170       Cmi[J]   =  Sum5/Suml 

4180       Cmii[Jj  =  Suml0/Sum9 

4190       Cml[J]   =  Sum6/Sum2 

4200       Cm2[J]   =  Sum7/Sum3 

4210       Cmo[J]   =  Sum8/Sum4 

4220       FOR  I  =  1  TO  Im 

4230  Ud[I,J]   =  Ud[I, J]*Areaf 

4240  Nd[I,J]  =  Nd[I,  Jl'-Areaf 

4250  Cd[I,J]   =  Ud[I, J]+Nd[I,J] 

4260  REM  (.■")  CONTINUE 

4261  NEXT  I 

4270  REM  CO  CONTINUE 

4271  NEXT  J 

4280  REM  OUTPUT  RESULTS 

4290  PRINT  #7  USING  4750:Ctrl$ 
4300  PRINT  #7  USING  4760: 
4305  PRINT  #7  USING  "4X,S": 

4310  FOR  I  =  1  TO  Im   I   PRINT  #7  USING  "4D,S":I    |   NEXT  I 

4315  PRINT  #7: 

4320  FOR  J  =  1  TO  Jm 

4330       PRINT  #7  USING  "3D,1X,S":J 

4340  FOR  I  =  1  TO  Im  I  PRINT  #7  USING  4780:Cd[I,J]    |   NEXT  I 

4341  PRINT  #7: 

4342  NEXT  J 

4350  IF  Fpat<=l  THEN  4630 

4360  REM  MULTIPLE  PATTERN  CODING 

4370  Npat  =  Fpat-1 

4380  Nran  =  Fran 

4390  Ncran  =  Fcran 

4400  FOR  I  =  1  TO  Im 

4410       FOR  J  =  1  TO  Jm 

4420  Nd[I,J]   =  Cd[I,J] 

4421  NEXT  J 
4  422  NEXT  I 
4430  Nc  =  0 
4440  Nr  =  0 

4450  FOR  N  =  1  TO  Npat 

4460       Nc  =  Nc+Ncran 

4470       Nr  =  Nr+Nran 

4  480       FOR  I  =  1  TO  Im 

4490  II  =  INT(I)-INT(Nc) 

4510  FOR  J  =  1  TO  Jm 

4520  Jl  =  J-Nr 

4530  IF  JKl  OR  Jl>Jm  THEN  4550 

4540  Cd[I,J]  =  Cd[I, J]+Nd[Il, Jl] 

4550  REM  (^')  CONTINUE 

4551  NEXT  J 

4560  REM  (*)  CONTINUE 

4561  NEXT  I 

4570  REM  (")  CONTINUE 

4571  NEXT  N 

4580  PRINT  #7  USING  4750: 
4590  PRINT  #7  USING  4640: 

4600  FOR  I  =  1  TO  Im   I   PRINT  #7  USING  4770:1    |   NEXT  I 


19 


4610  FOR  J  =  1  TO  Jm 

4620  FOR  I  =  1  TO  Im   I    PRINT  #7  USING  4780 ; J , Cd [ I . J ]    |   NEXT  I 

4621  NEXT  J 

4630  REM   (■■')  CONTINUE 

4640  IMAGE  4 7X , "MULTIPLE  PATTERN  DATA" 
4  650  REM  CALL  UNCOM(CD,GAL) 

4651  PRINT  "subroutine  UNCOM" 

4652  GOSUB  5250 

4655  IF  A[5]>0  THEN  GOSUB  20000 
4  660  REM  CALL  PPRINT(CD) 

4661  PRINT  "subroutine  PPRINT" 

4662  GOSUB  5990 
4670  Npts  =  0 

4680  FOR  J  =  1  TO  Jm 

4690       IF  Cd[14,J]=0  THEN  4730 

4700       Npts  =  Npts+1 

4710       Valin[J]   =  Cd[14,J] 

4720  REM   ('■')  CONTINUE 

4721  NEXT  J 

4730  REM   ('•')  CONTINUE 
4740  RETURN 
4  750  IMAGE  lA 

4760  IMAGE  4 7X , "DETAILED  PATTERN  DATA" 
4780   IMAGE  2D.1D,S 
4  7  90  END 

4800  REM  FUNCTION  CDFN(X) 
4810  DIM  Af [6] ,Df [6] ,Cf [6,10] 
4820  DIM  Ef[60] 

4830  REM   (")   EQUIVALENCE  (E(1),C(1,1)) 
4835  RESTORE  4840 

4840  DATA  0.625,1.25,2.0,2.45,3.5,4.62 

4841  FOR  J  =  1  TO  6    I    READ  Af[J]    |    NEXT  J 

4850  DATA  0.3,0.925,1.625,2.225,2.95,4.15 

4851  FOR  J  =  1  TO  6   |   READ  Df [J]    |   NEXT  J 

4860  DATA  0.00067982403291,-0.0012709753598,0.000679645257 

4865  DATA  -8 . 8500314297E-05 , 1 . 4791554152E-07 , 5 . 4879072878E-07 

4870  DATA  0 . 0051028640888 , -0 . 0014822649744 , -3 . 4589883732E-07 

4880  DATA  0 . 00056826070991 , -6 . 1965013125E-05 , -8 . 6656125337E-07 

4890  DATA  -0.00601608623  7  9,0.0087718965604,-0.00308834  01043 

4900  DATA  -0 . 000450 15531032 , 0 . 00022266710841 , 6 . 6877 7694 4 lE-07 

4910  DATA  -0.022725611373,-0.0020200415601,0.0072141082947 

4920  DATA  -0 . 0013872389725 , -0 . 00044390018986 , -6 . 4084653798E-07 

4930  DATA  0.033555772127,-0.034681674488,-0.001067703944 

4940  DATA  0 . 0055422185916 , 0 . 00062978750848 , 6 . 2471782478E-07 

4950  DATA  0.072703311825,0.047642403207,-0.024859750804 

4960  DATA  -0 . 010221129286 , -0 . 00066382630882 , -4 . 2105004094E-07 

4970  DATA  -0.14093895854,0.056850293056,0.057420315819 

4980  DATA  0 . 011850195831 , 0 . 00051265678121 , 2 . 0742133302E-07 

4990  DATA  -0.154  6891297,-0.22180615138,-0.065383  70517 

5000  DATA  -0.0088864030978,-0.00027657162532,-0.076856040218 

5010  DATA  0.5156304546,0.23979043073,0.040236129479 

5020  DATA  0 . 003993888570 1 , 9 . 3 751 4 1 3487E-05 , 1 . 8722020822E-08 

5030  DATA  0.164  31337971,0.4  0458836515,0.48922186659 

5040  DATA  0.4  9917416726,0.4  999848984  9,0.4  99999997  7  9 

5041  FOR  J  =  1  TO  60   |   READ  Ef[J]    |   NEXT  J 

5042  FOR  J  =  1  TO  10 

5043  FOR  K  =  1  TO  6 
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50UU  Cf[K,J]   =  Ef  [  (J-l)'''6+K] 

5045  NEXT  K 

5046  NEXT  J 

5050  Y  =  Xos'^0. 70710678119 
5060  Sgny  =  1 

5070  ON  SGN(Y)+2  GOTO  5100,5080,5120 

5080  Cdfn  =0.5 

5090  GOTO  3780 

5100  Sgny  =  -1 

5110  Y  =  -Y 

5120  FOR  K  =  1  TO  6 

5130       IF  Y-Af[K]<=0  THEN  EXIT  TO  5170 

5140  REM  (*)  CONTINUE 

5141  NEXT  K 
5150  Z  =  0.5 
5160  GOTO  5220 
5170  Y  =  Y-Df[K] 
5180  Z  =  Cf[K,l] 
5190  FOR  J  =  2  TO  10 
5200       Z  =  Z*Y+Cf[K,J] 
5210  NEXT  J 

5220  Cdfn  =  0.5+Sgny*Z 
5230  GOTO  3780 
5240  END 

5250  REM  SUBROUTINE  UNCOM( CD .GALLON) 
5260  REM 

5270  REM  COMPUTES  CONTOURS  FOR  SIMULATED  PATTERNS 

5280  REM 

5290  DIM  Cd[28,56] ,Btot [12] , Itot[12] 
5300  DIM  Cmean[12] 
5310  REM  (*)  REAL  ITOT 

5320  DATA  0.1,0.6,1.5,2.5,3.5,4.5,5.5,6.5,7.5,8.5,9.5,10. 
5330  RESTORE  5320 

5340  FOR  I  =  1  TO  12    |   READ  Cmean[I]    |   NEXT  I 

5350  Im  =  28 

5360  Jm  =  56 

5370  SumO  =  0 

5380  Suml  =  0 

5390  FOR  I  =  1  TO  12 

5400       Btot[I]   =  0 

5410  Itot[I]   =  0 

5411  NEXT  I 

5420  FOR  I  =  1  TO  Im 

5430      FOR  J  =  1  TO  Jm 

5440  Tot  =  Cd[I,J] 

5450  IF  Tot=0  THEN  5690 

5460  IF  Tot>0  AND  Tot<=0.2  THEN  5480 

5470  GOTO  5520 

5480  Btot[l]   =  Btot[l]+Tot 

5490  Itot[l]   =  Itot[l]+l 

5500  SumO  =  Suin0+Tot*4  .  5 

5510  GOTO  5690 

5520  IF  Tot>0.2  AND  Tot<=l  THEN  5540 

5530  GOTO  5580 

5540  Btot[2]   =  Btot[2]+Tot 

5550  Itot[2]   =  Itot[2]+l 

5560  SumO  =  Sum0+Tot*4.5 
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5570  GOTO  5690 

5580  FOR  K  =  1  TO  9 

5590  IF  Tor>K  AND  Tot<=K+l  THEN  5610 

5600  GOTO  5650 

5610  Btot[K+2]   =  Btot [K+2]+Tot 

5620  Itot[K+2]   =  Itot[K+2)+l 

5630  SumO  =  SumO+TofA  .  5 

56A0  EXIT  TO  5690 

5650  REM   (•")  CONTINUE 

5651  NEXT  K 

5660  Btot(12]   =  Btot[12]-(Tot 

5670  Itot(12]   =  Itot[12]+l 

5680  SumO  =  SumO+Tof*  .  5 

5690  REM   (")  CONTINUE 

5691  NEXT  J 

5692  NEXT  I 

5700  FOR  I  =  1  TO  12 

5710  Btot[Il   =  Btot[I]>*4.5 

5711  NEXT  I 

5720  Perr  =  SumO/Gai"100 
5730  FOR  I  =  1  TO  12 
5740       Itot[I]   =  Icot[I]"450 
5750       Suml  =  Suml+Itot [ I ] 


5751  NEXT  I 

5760  PRINT  #7  USING  5900: 
5765  PRINT  #7  USING  5910: 

5770  FOR  K  =  1  TO  11    |   PRINT  #7  USING  " 5D .  AD ,  3X  .  S"  :  Bt ot  [K ]    |   NEXT  K 

5775  PRINT  #7  USING  " IL ,  4X ,  3A ,  IL ,  50  .  iiD"  :  ">  10"  ,  Btot  [  12 ) 

5780  PRINT  #7  USING  5940 : SumO , Perr 

5790  Sum2  =  0 

5800  Sum3  =  0 

5810  FOR  I  =  1  TO  12 

5820       Sum2  =  Sum2+Btot [ I ] "Cmean ( I ] 

5830       Sum3  =  Sum3+Btoc[I] 

5840  REM   (")  CONTINUE 

5841  NEXT  I 

5850  Sujn2  =  Sum2/Sum3 

5860  PRINT  #7  USING  5960: Sum2 

5870  PRINT  #7  USING  5920: 

5879  PRINT  #7  USING  5930: 

5880  FOR  K  =  1  TO  11    |   PRINT  #7  USING  "9D , 3X , S" : Itot [K)    |   NEXT  K 

5881  PRINT  #7  USING   " IL , 4X , 3A , IL , 9D" : "> 10 " . It ot [ 12 ] 
5890  PRINT  #7  USING  5950: Suml 

5900  IMA  21 , 1 IX. "CONTOUR  INTERVALS  FOR  RETARDANT  DISTRIBUTION    (GALLGNS/100  SQUARE  FEET)   TOTAL  WITHIN  EACH  INTERVAL" 

5910   IMA  2X, "0,0-0.2  0.3-1.0  1.1-2.0  2.1-3.0  3.1-4.0  4.1-5.0  5.1-6.0  6.1-7.0  7.1-8.0  8.1 

5920  IMA  11 , 12X, "CONTOUR  INTERVALS  FOR  RETARDANT  DISTRIBUTION   (SQUARE-FOOT  COVERAGE  )   -  TOTAL  WITHIN  EACH  INTERVAL" 

5930  IKA  2X, "0.0-0. 2  0.3-1.0  1.1-2.0     .     2.1-3.0  3.1-4.0  4.1-5.0  5.1-6.0  6.1-7.0  7.1-8.0  8.1 

5940  IMA  2X, "TOTAL  AMOUNT  OF  RETARDANT  ON  GRID  =  ",7D.2D,"  GALLONS" , 1 1  OX , "PERCENT  OF  TOTAL  =  ",5D.3D 

5950   IMAGE  lOX, "TOTAL  AREA  COVERED   ",11D,"   SQUARE  FEET." 

5960  IMAGE  2X, "AVERAGE  PATTER-V  LEVEL  IS  ",4D.D,"  GALLONS  PER  100  SQV.^RE  FEET" 
5970  RETURN 
5980  END 

5990  REM  ("■)   SUBROUTINE  PPRINT(AA) 

6000  DIM  T[12] ,D[12] 

6001  DIM  SS[13] (132) ,P$[28] (132) 
6010  DATA  10,9,8,7,6,5,4,3,2,1,0.3,0 
6020  RESTORE  6010 

6030  FOR  I  =  1  TO  12    |    READ  7(1]    |    NEXT  I 

6040  DATA  "      10","        9","       8"."        7","        6"."       5","       4","        3","  2" 
6050  DATA  "       1" , "       +" , "        . " , " 
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6060 
6070 
6080 
6090 
6100 
6110 
6120 
6121 
6122 
6123 
612A 
6125 
6126 
6127 
6128 
6129 
6130 
6131 
6132 
6133 
6134 
6139 
6140 
6141 
6150 
6160 
6170 
6180 
6181 
6190 
6200 
6210 
6220 
6230 
6231 
6240 
6250 
6260 
6270 
6280 
6281 
6290 
6300 
6310 
6311 
6320 
6321 
6330 
6340 
6350 
6351 
6360 
6370 
6380 
6390 
6400 
6401 


RESTORE  6040 

FOR  I  =  1  TO  13    I   READ  S$[I]    |   NEXT  I 
Im  =  28 
Jm  =  56 
Y  =  -15 
Dy  =  30 

Ctrl$ 


PRINT 

a  7 

USING 

64  70  : 

PRINT 

if  1 

USING 

6480  : 

PRINT 

#7 

USING 

6481  : 

PRINT 

#7 

USING 

6482  : 

PRINT 

#7 

USING 

6483  : 

PRINT 

#7 

USING 

6484  : 

PRINT 

#7 

USING 

6485 : 

PRINT 

#7 

USING 

6486 : 

PRINT 

#7 

USING 

6487: 

PRINT 

#7 

USING 

6488: 

PRINT 

#7 

USING 

6489: 

PRINT 

#7 

USING 

6490: 

PRINT 

#7 

USING 

6491: 

PRINT 

#7 

USING 

6492: 

PRINT 

#7 

USING 

6493: 

PRINT 

#7 

USING 

"4X,S 

FOR  I 

1  TO  Im   1  PR 

PRINT 

#7 

USING 

"3L"  : 

Jst  = 

0 

FOR  J 

1  TO  Jm 

FOR 

I 

=  1  TO 

Im 

NEXT  I 


P$[I]   =  S$[13] 
NEXT  I 
1st  =  0 

FOR  I  =  1  TO  Im 
FOR  K  =  1  TO  12 

IF  Cd[I,J]>T[K]   THEN  EXIT  TO  6250 
REM  ('■')  CONTINUE 
NEXT  K 
GOTO  6280 
P$[I]   =  S$[K] 
1st  =  1 
Jst  =  1 

REM  (")  CONTINUE 
NEXT  I 

IF  Jst=l  AND  lst=0  THEN  EXIT  TO  6330 
Y  =  Y+Dy 

PRI  #7  USI  6510:J   |   FOR  II  =  1  TO  Im   |   PRI  #7  USI  6511:P$[I1]    |    NEX  II 
PRINT  #7  USING  6512 :Y 
REM  (")  CONTINUE 
NEXT  J 

REM  (")  CONTINUE 
FOR  I  =  1  TO  12 

D[I]  =  0 
NEXT  I 
I  =  14 

FOR  J  =  1  TO  Jm 
FOR  K  =  1  TO  12 

IF  Cd[I,J]>T[K]   THEN  D[K)   =  D[K]+Dy 
REM   (•■<)  CONTINUE 
NEXT  K 
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6410  REM  CO  CONTINUE 

6411  NEXT  J 

6420  PRINT  #7  USING  6520: 

6430  FOR  I  =  1  TO  12 
6440       K  =  12-1+1 

6450  PRINT  #7  USING  6530 : D [K] , T [K] 

6451  NEXT  I 
64  60  RETURN 


6A  70 

IMAGE 

lA 

6^80 

IMAGE 

40X, 

'S  C  A 

L  E 

D 

P 

A 

T  T  E  R  N 

P  L  0  T  " ,1L 

5481 

IMAGE 

38X, "CELL 

SIZE 

IS  15 

FEET  WIDE  BY 

30  FEET 

LONG" 

6  A  82 

IrlAljii 

38X, 

-  0  . 

0 

TO 

0 

2 

GALLONS/ 100 

SQUARE 

FEET" 

6483 

IMAGE 

38X, 

'  +  -- 

-  0. 

3 

TO 

1 

0 

GALLONS/ 100 

SQUARE 

FEET" 

6484 

IMAGE 

38X, 

'   1  -- 

-  1. 

1 

TO 

2 

0 

GALLONS/ 100 

SQUARE 

FEET" 

6485 

IMAGE 

38X, 

'  2 

-  2. 

1 

TO 

3 

0 

GALLONS/ 100 

SQUARE 

FEET" 

6486 

IMAGE 

38X, 

'  3 

-  3. 

1 

TO 

4 

0 

GALLONS/ 100 

SQUARE 

FEET" 

6487 

IMAGE 

38X, 

'   4  -- 

-   4  . 

1 

TO 

5 

0 

GALLONS/100 

SQUARE 

FEET" 

6488 

IMAGE 

38X, 

'  5  -- 

-  5. 

1 

TO 

6 

0 

GALLONS/ 100 

SQUARE 

FEET" 

6489 

IMAGE 

38X, 

'  6  -- 

-  6. 

1 

TO 

7 

0 

GALLONS/ 100 

SQUARE 

FEET" 

6490 

IMAGE 

38X, 

'  7  -- 

-   7  . 

1 

TO 

8 

0 

GALLONS/100 

SQUARE 

FEET" 

6491 

IMAGE 

38X, 

'   8  -- 

-  8. 

1 

TO 

9 

0 

GALLONS/ 100 

SQUARE 

FEET" 

6492 

IMAGE 

38X, 

'  9  -- 

-  9. 

1 

TO 

10. 

GALLONS/100 

SQUARE 

FEET" 

6493   IMAGE  38X,"   10    GREATER  THAN  10  GALLONS/100  SQUARE  FEET" 

6510  IMAGE  3D,1X,S 

6511  IMAGE  4A,S 

6512  IMAGE  5D , "  FEET",4L 

6520  IMAGE  48X, "PATTERN  LENGTH  DATA     " , IL 

6530  IMAGE  35X,5D,   "FEET  ABOVE  ",4D.D,"  GALLONS/100  SQUARE  FEET" 
6540  END 

6550  REM  ('M   SUBROUTINE  PATADD(VACI) 

6551  DELETE  Ds t , Qo , Qdoto , Zo , Veo , Vzo   !   QDOTX , TIMl , VEX 
6560  DIM  Tp[2000] ,Lcon[9] ,Conv[9]    !  TP[13440] 

6570  DIM  Sing[1680] 

6580  DIM  Lcmax[9] ,Dela[9] 

6590  REM  ('•')   COMMON   / ADD / NPTS  ,  VALIN (  56 )  ,  ITC ( 8 ) 
6600  DATA  0.5,1,2,3,4,6,8,10,12 
6610  RESTORE  6600 

6620  FOR  I  =  1  TO  9    |   READ  Conv[I]    |    NEXT  I 

6630  IMA  lA, 36X, "MAXIMUM  LINE  LENGTHS  WITH  CORRESPONDING  TANK  OPENING  DELAYS   " , 2L 

6640  IMA  iL,21X,"  NO.   OF  ",8(2X,"  "),1L,21X,"  TANKS     ",8("       FEET  SEC" 

6660   IMAGE  2X 

6670  Knots  =  Vaci 

6680  Valin[Npts+l]   =  0 

6690  Valin[Npts+2]  =  0 

6700  Jfrst  =  0 

6710  REM  WRITE  HEADINGS 

6720  PRINT  #7  USING  6630:Ctrl$ 

6729  PRINT  #7  USING  "34X,S": 

6730  FOR  I  =  1  TO  8    |    PRINT  #7  USING  "2X, 3D . D , 4A , S " : Conv [ I ] , "  GPC"    |    NEXT  I 

6731  PRINT  #7  USING  6640: 
6740  Lalt  =  0 

6750  REM  LENGTH  OF  SUBCELL 
6760  Mcell  =  5 

6770  REM  NUMBER  OF  SUBCELLS  TO  START 
6780  Mstcl  =  3 

6790  REM  NUMBER  OF  SUBCELLS  TO  A  CELL 
6800  Mlncl  =  6 
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6810  Fps  =  Knots''!  .  68781 

6820  REM  DELT  IS  LENGTH     (SEC)  OF  INTERVAL 

6830  Delt  =  0.1 

6840  Fpdel  =  Fps'^Delt 

6850  REM  COMPUTE  NUMBER  OF  CELLS  PER  INTERVAL 

6860  Idcell  =  INT ( Fpde I /Mce 1 1+0 . 5 ) 

6870  REM  EXPAND  SINGLE  PATTERN 

6880  Qnext  =  0 

6890  Lastd  =  0 

6900  Nextd  =  Mstci 

6910  Nto  =  Npts+2 

6920  FOR  I  =  1  TO  Nto 

6930       Qlast  =  Qnext 

6940      Qnext  =  Valin[I] 

6950       Slope  =  (Qnext-Qlast ) / (Nextd-Lastd) 

6960       Jx  =  Lastd+1 

6970       FOR  Ix  =  Jx  TO  Nextd 

6980  Sing[Ix]   =  Qlast+( Ix-Lastd) "Slope 

6981  NEXT  Ix 

6990      Lastd  =  Nextd 

7000  Nextd  =  Nextd+Mlncl 

7001  NEXT  I 

7010  Lsp  =  INT(Npts'>Mlncl+Mlncl/2) 
7020  REM  MAXIMUM  NUMBER  OF  INTERVALS 
7030  Nintv  =  INT (Lsp / Idcel 1+1 ) 
7040  REM  SINGLE  PATTERN  IS  EXPANDED 
7050  REM 

7060  REM  ROUTINE  TO  ADD  TOGETHER  PATTERNS 

7070  Ntank  =  0 

7080  FOR  Jj  =  1  TO  8 

7090       Ntank  =  Ntank+1 

7100       FOR  Kc  =  1  TO  9 

7110  Dela[Kc]   =  0 

7120  Lcmax[Kc]  =  0 

7121  NEXT  Kc 

7130       IF  Itc[Jj]<>l  THEN  7410 
7140       FOR  Jr  =  1  TO  Nintv 
7150  J  =  Jr-1 

7160  Idisp  =  INT(J*Idcell) 

7170  Dly  =  J-^Delt 

7180  Maxl  =  INT(Lsp+(Ntank-l)'-'Idisp) 

7190  FOR  K  =  1  TO  Maxl 

7200  Tp[K]  =  0 

7201  NEXT  K 
7210  Idat  =  0 

7220  FOR  K  =  1  TO  Ntank 

7230  FOR  Kl  =  1  TO  Lsp 

7240  Kk  =  INT(K1+Idat) 

7230  Tp[Kk]  =  Tp[Kk]+Sing[Kl] 

7251  NEXT  Kl 

7260  Idat  =  INT(Idat+Idisp) 

7261  NEXT  K 

7270  REM  ROUTINE  TO  COMPUTE  LENGTHS 

7280  FOR  I  =  1  TO  9 

7290  Lcon[I]  =  0 

7291  NEXT  I 

7300  FOR  I  =  1  TO  Maxl 
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7310  FOR  J  =  1  TO  9 

7320  IF  Tp[I]>=Conv[ J]   THEN  Lcon[J]  =  INT (Leon [ J] +Mcel 1 ) 

7321  NEXT  J 

7330  REM  (")  CONTINUE 

7331  NEXT  I 

7340  FOR  K  =  1  TO  9 

7350  IF  Lcmax[K]>=Lcon[K]   THEN  7380 

7360  Lcmax[K]   =  INT(Lcon[K]) 

7370  Dela[K]   =  Dly 

7  380  REM   (")  CONTINUE 

7381  NEXT  K 

7390  REM   ('")  CONTINUE 

7391  NEXT  Jr 

7400  PRINT  #7  USING  " 22X , 7X , ID , 4X , S" : Ntank 

7401  FOR  I  =  1  TO  8    I    PRI  #7  USI   " 2X , 4D , IX , 2D . D , S" : Lcmax [ I ] , Dela [ I ] 

7402  PRINT  #7 : 

7410  REM   ('■')  CONTINUE 

7411  NEXT  Jj 

7420  REM  END  OF  THIS  ALTITUDE  AND  VELOCITY 

7430  RETURN 

7440  END 

10000  Gali  =  A[17] 

10010  Cflag  =  A[18] 

10020  Rdqv  =  A[24] 

10030  Convex  =  A[25] 

10040  Bflag  =  A[32] 

10050  Alf  =  A[37] 

10060  Kl  =  A[40] 

10070  K2  =  A[41] 

10080  Delt  =  A[42] 

10090  Tr  =  A[43] 

10100  Tx  =  A[44] 

10110  Rflag  =  A[45] 

10120  Fpat  =  A[46] 

10130  Fran  =  A[47] 

10140  Fcran  =  A[48] 

10150  Vaci  =  A[51] 

10160  Altl  =  A[52] 

10170  Add  =  A[53] 

10180  R  =  A[62] 

10190  RETURN 

20000  OPEN  #5:C$, "f" 

20010  WRITE  #5:M,N$,K1, Vaci, Altl, K2,Delz, Sigma, Cd, Gali 

20020  CLOSE 
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APPENDIX  B:  SAMPLE  PATSIM  INPUT-OUTPUT 
B-1  Sample  input  file 

PATSIM  INPUT 

17 
500 
18 
1 

51 

125 

52 

200 

45 

1 

53 

1 

1 

1 

0 

1 

0 

0 

0 

0 

24 
10 


0 

0 

0 

0 

1 

2.68 

0 

0 

2 

7.98  0 

0 

3 

13.22 

0 

0 

4 

23.52 

0 

0 

5 

36 

0 

0 

6 

52.22 

0 

0 

7 

68.92 

0 

0 

8 

84.68 

0 

0 

9 

100.24 

0 

*  DEMPSAY  DC-4  4"  OPEN  200'   125  KT  GUM 

B-2  Sample  output 

STANDARD    PATTERN  SIMULATION 

RETARDANT  PATTER 


17  5.000000E+02 

18  1.000000E+00 

51  1.250000E1-02 

52  2.000000E+02 
45  l.OOOOOOE+00 

53  I.OOOOOOE^OO 
24  1.000000E*01 

QUANTITY  AND  DOOR  ANGLE  VS  TIME  DATA 

TIME-SEC  Q-GAL  ANGLE-DEG 


O.OOOOOE+00 

O.OOOOOE+00 

O.OOOOOE+00 

l.OOOOOE 

01 

2.68000E+00 

O.OOOOOEfOO 

2.00000E 

01 

7.98000E+00 

O.OOOOOE+00 

3.00000E 

01 

1.32200E*01 

O.OOOOOE+OO 

4.00000E 

01 

2.35200E+01 

O.OOOOOE+00 

5.00000E 

01 

3.60000E+01 

O.OOOOOE+00 

6.00000E 

01 

5.22200E+01 

O.OOOOOE+00 

7.00000E 

01 

6.89200E+01 

O.OOOOOE+00 

8.00000E 

01 

8.46800E+01 

O.OOOOOE+00 

9.00000E 

01 

1.00240E+02 

O.OOOOOE+00 

DEMPSAY  DC-4  4" 

OPEN  200'  125 

VEL0CITY(KNOTS)= 

125 

ALTITUOE(FEET)=  200 
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DETAILED  PATTERN  DATA 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

1 

0. 

.0 

0. 

.0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

0 

0. 

,0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

1 

0. 

1 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

2 

0. 

.0 

0. 

.0 

0. 

,0 

0. 

0 

0. 

,0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

2 

0. 

3 

0. 

3 

0. 

2 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

3 

0. 

.0 

0. 

,0 

0. 

.0 

0. 

0 

0. 

.0 

0. 

0 

0. 

,0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

1 

0. 

2 

0. 

3 

0. 

3 

0. 

2 

0. 

1 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

4 

0. 

.0 

0. 

,0 

0, 

.0 

0. 

0 

0. 

.0 

0. 

0 

0. 

.0 

0. 

0 

0. 

0 

0. 

0 

0. 

1 

0. 

1 

0. 

,1 

0. 

2 

0. 

2 

0. 

1 

0. 

1 

0. 

1 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

5 

0. 

.0 

0. 

,0 

0. 

,0 

0. 

0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

0 

0. 

0 

0. 

0 

0. 

1 

0. 

1 

0. 

,2 

0. 

,4 

0. 

4 

0. 

2 

0. 

1 

0. 

1 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

6 

0. 

.0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

.0 

0. 

,0 

0. 

0 

0. 

,0 

0. 

0 

0. 

1 

0. 

1 

0. 

3 

0. 

,5 

0. 

,5 

0. 

3 

0. 

1 

0. 

1 

0. 

,0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

7 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

0 

0. 

0 

0. 

1 

0. 

1 

0. 

2 

0. 

,4 

0. 

.6 

0. 

,6 

0. 

4 

0. 

2 

0. 

1 

0. 

1 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

0. 

0 

8 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

0 

0. 

,0 

0. 

1 

0. 

1 

0. 

2 

0. 

,5 

0. 

,9 

0. 

,9 

0. 

,5 

0. 

2 

0. 

1 

0. 

,1 

0. 

0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

0 

0. 

0 

0. 

,0 

0. 

0 

0. 

0 

9 

0. 

,0 

0. 

0 

0. 

,0 

0. 

0 

0. 

.0 

0. 

0 

0. 

.0 

0. 

0 

0. 

,0 

0. 

1 

0. 

1 

0. 

3 

0. 

,7 

1 . 

,1 

1 , 

.1 

0. 

,7 

0. 

,3 

0. 

1 

0. 

,1 

0. 

0 

0. 

,0 

0. 

0 

0. 

,0 

0. 

,0 

0. 

0 

0. 

,0 

0. 

,0 

0. 

.0 

10 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

0 

0. 

,0 

0. 

0 

0. 

,0 

0. 

,0 

0. 

.0 

0. 

1 

0. 

,1 

0. 

,3 

0. 

.7 

1 . 

.1 

1 . 

,1 

0. 

.7 

0. 

,3 

0. 

1 

0. 

.1 

0. 

,0 

0. 

.0 

0. 

.0 

0. 

,0 

0. 

.0 

0. 

,0 

0. 

.0 

0. 

.0 

0, 

.0 

11 

0. 

.0 

0. 

.0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

1 

0, 

,1 

0. 

,2 

0. 

,6 

1 . 

.0 

1 . 

,0 

0. 

,6 

0. 

.2 

0. 

1 

0. 

,1 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

.0 

0. 

,0 

0. 

.0 

12 

0. 

,0 

0. 

.0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

.0 

0. 

,1 

0. 

,1 

0, 

,2 

0. 

.6 

1 . 

,0 

1 . 

.0 

0. 

.6 

0. 

,2 

0. 

.1 

0. 

,1 

0. 

,0 

0, 

.0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

.0 

0. 

.0 

0, 

.0 

13 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,1 

0. 

1 

0. 

,1 

0. 

,2 

0, 

.6 

1 . 

,0 

1 . 

.0 

0. 

,6 

0. 

,2 

0. 

,1 

0. 

,1 

0. 

,1 

0. 

.0 

0. 

,0 

0. 

.0 

0. 

,0 

0. 

,0 

0. 

.0 

0. 

.0 

0. 

.0 

14 

0. 

,0 

0. 

.0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,1 

0. 

,1 

0. 

.1 

0. 

,2 

0. 

.6 

1 . 

.0 

1 . 

.0 

0, 

.6 

0. 

.2 

0. 

,1 

0. 

.1 

0. 

,1 

0, 

.0 

0. 

.0 

0. 

.0 

0. 

.0 

0. 

.0 

0, 

.0 

0. 

.0 

0, 

,0 

15 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

.0 

0. 

.0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,1 

0, 

,1 

0. 

.1 

0. 

.2 

0. 

.6 

1 . 

.0 

1 . 

.0 

0. 

,6 

0. 

,2 

0. 

,1 

0. 

.1 

0. 

.1 

0. 

.0 

0. 

.0 

0. 

.0 

0. 

.0 

0. 

.0 

0, 

.0 

0. 

,0 

0, 

.0 

16 

0. 

,0 

0. 

,0 

0. 

.0 

0. 

,0 

0. 

.0 

0. 

,0 

0. 

.0 

0. 

,0 

0. 

.1 

0. 

,1 

0. 

,1 

0. 

.2 

0. 

,6 

1 . 

.0 

1 . 

.0 

0. 

.6 

0. 

,2 

0. 

,1 

0. 

.1 

0. 

.1 

0, 

,0 

0. 

,0 

0, 

.0 

0. 

.0 

0, 

.0 

0, 

.0 

0, 

.0 

0, 

.0 

17 

0. 

,0 

0. 

.0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,1 

0. 

,1 

0. 

,1 

0. 

,2 

0. 

,6 

1 . 

.0 

1 . 

,0 

0. 

,6 

0. 

,2 

0. 

,1 

0. 

.1 

0. 

.1 

0, 

.0 

0, 

.0 

0, 

.0 

0. 

.0 

0, 

.0 

0, 

.0 

0, 

.0 

0, 

.0 

18 

0. 

,0 

0. 

.0 

0. 

,0 

0. 

0 

0. 

.0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

.0 

0, 

,1 

0. 

,1 

0. 

,2 

0, 

.6 

1. 

.0 

1 . 

.0 

0. 

.6 

0. 

,2 

0. 

.1 

0. 

.1 

0. 

,0 

0, 

.0 

0, 

.0 

0, 

.0 

0, 

,0 

0. 

.0 

0, 

.0 

0, 

.0 

0, 

.0 

19 

0. 

.0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

.0 

0. 

0 

0. 

,0 

0. 

,0 

0. 

.0 

0. 

,1 

0. 

,  1 

0. 

,2 

0, 

.6 

1 , 

.0 

1 , 

.0 

0. 

.6 

0. 

,2 

0. 

,1 

0. 

.1 

0. 

.0 

0, 

.0 

0, 

.0 

0, 

.0 

0, 

.0 

0. 

.0 

0, 

.0 

0, 

.0 

0 

.0 

20 

0. 

.0 

0. 

.0 

0. 

,0 

0. 

,0 

0. 

.0 

0. 

,0 

0. 

.0 

0. 

.0 

0. 

.0 

0. 

.1 

0. 

,1 

0, 

,2 

0, 

.6 

1 , 

.0 

1 , 

.0 

0. 

.6 

0. 

,2 

0. 

.1 

0, 

.1 

0. 

.0 

0, 

.0 

0, 

.0 

0, 

.0 

0, 

,0 

0. 

,0 

0 

.0 

0, 

.0 

0 

.0 

21 

0. 

.0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

.0 

0. 

,0 

0, 

.0 

0. 

,1 

0. 

.1 

0. 

.2 

0, 

.6 

1 . 

.0 

1 , 

.0 

0, 

.6 

0. 

.2 

0. 

.1 

0, 

,1 

0. 

.0 

0, 

.0 

0, 

.0 

0, 

.0 

0, 

.0 

0, 

.0 

0 

.0 

0 

.0 

0 

.0 

22 

0. 

,0 

0. 

.0 

0. 

.0 

0. 

0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

,0 

0, 

.0 

0. 

,1 

0. 

,1 

0. 

,2 

0, 

.6 

1 , 

.0 

1 , 

.0 

0, 

.6 

0. 

.2 

0. 

.1 

0, 

.1 

0, 

.0 

0 

.0 

0, 

.0 

0, 

.0 

0, 

,0 

0, 

.0 

0 

.0 

0 

.0 

0 

.0 

23 

0. 

,0 

0. 

.0 

0. 

.0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

.0 

0. 

,0 

0. 

.0 

0. 

,  1 

0. 

,1 

0. 

.2 

0, 

.6 

1 , 

.0 

1, 

.0 

0. 

.6 

0. 

.2 

0. 

.1 

0, 

,1 

0, 

.0 

0 

.0 

0 

.0 

0, 

.0 

0, 

.0 

0, 

.0 

0 

.0 

0, 

.0 

0 

.0 

24 

0. 

.0 

0. 

,0 

0. 

,0 

0. 

,0 

0. 

.0 

0. 

,0 

0, 

,0 

0. 

,0 

0. 

.0 

0. 

,0 

0. 

,1 

0. 

,2 

0, 

.6 

1 , 

.1 

1 , 

.1 

0, 

.6 

0. 

.2 

0. 

.1 

0, 

.0 

0. 

.0 

0 

.0 

0, 

.0 

0, 

.0 

0, 

.0 

0, 

.0 

0 

.0 

0, 

.0 

0 

.0 

25 

0. 

.0 

0. 

.0 

0. 

,0 

0. 

.0 

0. 

.0 

0. 

,0 

0, 

,0 

0. 

,0 

0, 

.0 

0. 

,0 

0. 

,1 

0. 

.2 

0, 

.6 

1 , 

.1 

1 , 

.1 

0, 

.6 

0. 

.2 

0. 

.1 

0, 

.0 

0, 

.0 

0 

.0 

0 

.0 

0, 

.0 

0 

.0 

0, 

.0 

0 

.0 

0 

.0 

0 

.0 

26 

0. 

.0 

0. 

,0 

0. 

.0 

0. 

,0 

0, 

.0 

0. 

.0 

0, 

.0 

0. 

.0 

0, 

.0 

0. 

,0 

0. 

.1 

0, 

.2 

0, 

.6 

1 , 

.0 

1, 

.0 

0, 

.6 

0, 

.2 

0. 

,1 

0 

.0 

0, 

.0 

0 

.0 

0 

.0 

0 

.0 

0 

.0 

0, 

.0 

0 

.0 

0 

.0 

0 

.0 

27 

0. 

,0 

0. 

.0 

0, 

.0 

0. 

.0 

0, 

.0 

0. 

.0 

0, 

.0 

0, 

.0 

0 

.0 

0, 

.0 

0. 

.0 

0, 

.2 

0 

.5 

0, 

.8 

0, 

.8 

0, 

.5 

0, 

.2 

0, 

.0 

0, 

.0 

0, 

.0 

0 

.0 

0 

.0 

0 

.0 

0 

.0 

0 

.0 

0 

.0 

0 

.0 

0 

.0 

28 

0. 

.0 

0. 

,0 

0, 

.0 

0. 

,0 

0, 

.0 

0. 

.0 

0, 

.0 

0, 

,0 

0, 

.0 

0. 

,0 

0. 

.0 

0. 

.1 

0, 

.3 

0, 

.5 

0, 

.5 

0, 

.3 

0, 

.1 

0. 

.0 

0, 

.0 

0, 

.0 

0 

.0 

0 

.0 

0 

.0 

0 

.0 

0 

.0 

0 

.0 

0 

.0 

0 

.0 

29 

0, 

.0 

0. 

.0 

0, 

.0 

0, 

.0 

0, 

.0 

0, 

.0 

0, 

.0 

0, 

.0 

0, 

.0 

0. 

.0 

0. 

.0 

0, 

.0 

0 

.1 

0, 

.2 

0, 

.2 

0, 

.1 

0, 

.0 

0, 

.0 

0, 

.0 

0, 

.0 

0 

.0 

0 

.0 

0 

.0 

0 

.0 

0, 

.0 

0 
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CONTOUR  INTERVALS  FOR  RETARDANT  DISTRIBUTION  {GALLONS/100  SQUARE  FEET)  TOTAL  WITHIN  EACH  INTERVAL 
0.0-0.2         0.3-1.0         1.1-2.0    .    2.1-3.0         3.1-4.0         4.1-5.0         5.1-6.0         6.1-7.0         7.1-8.0         8.1-9.0  9.1-10 
47.7202         212.1845         160.6453  0.0000  0.0000  0.0000  0.0000  0.0000  0.0000  0.0000 

TOTAL  AMOUNT  OF  RETARDANT  ON  GRID  =  420.55  GALLONS 

AVERAGE  PATTERN  LEVEL   IS        0.9  GALLONS  PER  100  SQUARE  FEET 

CONTOUR   INTERVALS  FOR  RETARDANT  DISTRIBUTION  (SQUARE-FOOT  COVERAGE  )  -  TOTAL  WITHIN  EACH  INTERVAL 
0.0-0.2         0.3-1.0         1.1-2.0    .    2.1-3.0         3.1-4.0         4.1-5.0         5.1-6.0         6.1-7.0         7.1-8.0         8.1-9.0  9.1-10.0 
93600  47700  15300  00000000 

TOTAL  AREA  COVERED  156600  SQUARE  FEET. 

SCALED      PATTERN  PLOT 
CELL  SIZE  IS  15  FEET  WIDE  BY  30  FEET  LONG 
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195  FEET 


29 


8  ....++++....  225  FEET 

9  ....♦11+....  255  FEET 

10  ....+11+....  285  FEET 

11  ....+11+....  315  FEET 

12  ....+11+....  345  FEET 

13  ....+11+....  375  FEET 

14  ....+11+....  405  FEET 

15  ....+11+....  435  FEET 

16  ....+11+....  465  FEET 

17  ....+11+....  495  FEET 

18  ....+11+....  525  FEET 

19  ....+11+....  555  FEET 

20  ....+11+....  585  FEET 

21  ....+11+....  615  FEET 


30 


22 


+      11  + 


645  FEET 


23  ....+11+....  675  FEET 

24  ....+11+....  705  FEET 

25  ....+11+....  735  FEET 

26  ....++++....  765  FEET 

27  ....++++....  795  FEET 

28  ....++++....  825  FEET 

29    855  FEET 


PATTERN  LENGTH  DATA 
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MAXIMUM  LINE  LENGTHS  WITH  CORRESPONDING  TANK  OPENING  DELAYS 


0.5  GPC  1.0  GPC  2.0  GPC  3.0  GPC  4.0  GPC  6.0  GPC  8.0  GPC  10. 'I  GPC 

1  670    0.0  525  0.0  0  0.0  0  0.0  0  0.0  0  0.0  0  0.0  0  0.0 

2  1390  3.6  1145  3.1  525  0.0  0  0.0  0  0.0  0  0.0  0  0.0  0  0.0 
4           2830    3.6  2385  3.1  845  1.6  590  0.0  525  0.0  0  0.0  0  0.0  0  0.0 
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APPENDIX  D:  BASIC  PROGRAMS 


PATSUM 
LLDELAY.TAB 
LDR  WIDTH 
BESTSTRAT.DAT 
BS  Plot 

SMOOTHFOOT 
FINAL  FOOT 


Adds  dissimilar  tanks. 

Tabulates  line  length  and  delay  data. 

Calculates  line  length  down  range  and  pattern  width  from  PATSIM. 
Data  portion  of  best  strategy  plot. 
Plots  best  strategy  graphs. 

Smooths  footprint  data  prior  to  the  use  of  FINALFOOT. 
Plots  footprint  graphs. 


PATSUM 


PATSUM  is  used  to  add  nonequal  patterns  to  simulate  trail  drops. 

100  PRINT   "patadd  for  2  DISSIMILAR  tanki 

105  SET  DIA  COL  3    |    SET  TEX  COL  14    |    SET  TEX  SIZ   4.0    |    SET  TEX  STY  0 

106  T$  =  CHR$(27) 

107  T$  =  T$  Sc  "MQl"    I    T$  =  "echo  "   &  T$ 

108  CALL  T$ 

109  INPUT  PROMPT  "what  is  the  title?  ":Tit$ 

110  INPUT  PROMPT  "what   is  the  airspeed?  ":Knots 
120  DIM  Valinl[56] ,Valin2[56] ,Itc[16] 

130  DIM  Tp[13440] ,Lcon[9] ,Conv[9] ,Aa[8] ,Singl[1680] ,Lcmax[9] ,Xdel[9] 
135  DIM  Sing2[1680] ,Singx[1680] 

140  DATA  0.5,1,2,3,4,6,8,10,12    |   RESTORE  140 

141  Itc  =  0    I    Itc[l]    =  1    I    Itc[2]   =  1    I    Itc[4]   =  1    I    Itc[8]   =  1    |    Itc[16]   =  1 
150  FOR  I  =  1  TO  9   I  READ  Conv[I]    |   NEXT  I 

160  PRINT  "ENTER  DATA  FOR  DROP  1  CENTER  (ROW  14)  ROW  ENTER  -1  TO  STOP" 
170  Valinl  =  0    |  Valin2  =  0 

180  FOR  I  =  1  TO  56    |   PRINT  "ENTER  VALUE  FOR  ROW  ":I 

190       INPUT  Q   I   IF  Q<0  THEN  EXIT  TO  200   |  Valinl [I]   =  Q  |   NEXT  I 

200  Nptsl  =  I-l 

210  PRINT  "ENTER  DATA  FOR  DROP  2  CENTER  (COL  14)   COLUMN     ENTER  -1  TO  STOP" 

220  FOR  I  =  1  TO  56    |   PRINT  "ENTER  VALUE  FOR  ROW  " ; I 

230       INPUT  Q   I    IF  Q<0  THEN  EXIT  TO  240    |    Valin2[I]   =  Q   |    NEXT  I 

240  Npts2  =  I-l 

250  Jfirst  =  0 

260  Lalt  =  0 

261  lalt  =  1 

270  Mcell  =  5    !   LENGTH  OF  SUBCELL 
280  Mstcl  =  3    !   NUMBER  OF  SUBCELLS  TO  START 
290  Mlncl  =  6    !   NUMBER  OF  SUBCELLS  TO  A  CELL 
300  Fps  =  Knots'^  . 68781 

310  Delt  =0.1    !   LENGTH   (SEC)   OF  INTERVAL 
320  Fpdel  =  Fps'''Delt 

330  Idcell  =  Fpdel/Mcell+0.5    !   COMPUTE  NUMBER  OF  CELLS  PER  INTERVAL 

340  REM  EXPAND  SINGLE  PATTERN 

350  Ntol  =  Nptsl+2   I   Nto2  =  Npts2+2 

355  Qnextl  =  0    |   Qnext2  =  0    |   Xlastdl  =  0    |   Xlastd2  =  0    |   Xnextdl  =  Mstcl    |   Xnextd2  =  Mstcl 
360  FOR  I  =  1  TO  Ntol 
370      Qlastl  =  Qnextl 
380      Qnextl  =  Valinl [I] 

390       Slopel  =  (Qnextl-QIastl) / (Xnextdl-Xlastdl) 
400       Jxl  =  Xlastdl+1 

410       FOR  J  =  Jxl  TO  Xnextdl    |    Singl[J]   =  Qlastl  +  (  J-Xlastdl )  •''Slopel    |    NEXT  J 

420      Xlastdl  =  Xnextdl 

430       Xnextdl  =  Xnextdl+Mlncl 
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440  NEXT  I 

450  Lspl  =  Nptsl'''Mlncl+Mlncl/2 

460  Nintvl  =  Lspl / IdceLl+1 

470  REM  SINGLE  PATTERN  1   IS  EXPANDED 

480  FOR  I  =  1  TO  Nto2 

490  Qlast2  =  Qne}tt2 

500  Qnext2  =  Valin2[I] 

510  Slope2  =  (Qnext2-Qlast2) / (Xnextd2-Xlastd2) 

520  Jx2  =  Xlastd2+1 

530  FOR  J  =  Jx2  TO  Xnextd2    |    Sing2[J]   =  Qlast2+( J-Xlastd2)'''Slope2    |   NEXT  J 

540  Xlastd2  =  Xnextd2 

550  Xnextd2  =  Xnextd2+Mlncl 

560  NEXT  I 

570  Lsp2  =  Npts2"MlncH-Mlncl/2 

580  Nintv2  =  Lsp2/Idcell+1 

590  REM  SINGLE  PATTERN  2  IS  EXPANDED 

1000  REM  ROUTINE  TO  ADD  TOGETHER  PATTERNS 

1001  CLEAR 

1002  Xpace  =  1 

1005  MOVE  0,95   |   SET  TEXT  COLOR  1   |   TEXT  Tit$ 

1010  Ntankl  =  0 

1015  Xdel  =  0   I   Lcmax  =  0 

1020  FOR  Jj  =  1  TO  16 

1030  Ntankl  =  Ntankl+1 

1040  Xdel  =  0    I   Lcmax  =  0 

1050  IF  Itc[Jj]<>l  THEN  2000 

1051  K  =  1 

1052  GOSUB  5000 

1060  FOR  Jr  =  1  TO  Nintvx 

1065  SET  DIALOG  COLOR  13    |    PRINT  Jr;"   ";Nintvx   |    SET  DIALOG  COLOR  3 

1070  J  =  Jr-1 

1080  Idisp  =  J"'Idcell 

1090  Dly  =  JADelt 

1100  Maxl  =  Lspx+(Ntankl-l)AIdisp 

1110  Tp  =  0 

1120  Idat  =  0 

1130  FOR  K  =  1  TO  Ntankl 

1135  GOSUB  5000 

1140  FOR  Kl  =  1  TO  Lspx 

1150  Kk  =  KL+Idat 

1160  Tp[Kk]   =  Tp[Kk]+SLngx[Kl] 

1170  NEXT  Kl 

1180  Idat  =  Idat+Idisp 

1190  NEXT  K 

1200  REM  ROUTINE  TO  COMPUTE  LENGTHS 

1205  SET  TEXT  COLOR  13 

1210  Leon  =  0 

1220  FOR  I  =  1  TO  Maxl 

1230  FOR  J  =  1  TO  9 

1240  IF  Tp[I]>=Conv[J]   THEN  Lcon[J]   =  Leon [ J ] +Mcel 1 

1250  NEXT  J 

1260  NEXT  I 

1270  FOR  K  =  1  TO  9 

1280  IF  Lcmax[K]>=Leon[K]   THEN  1500 

1290  Lcmax[K]  =  Leon[K] 

1300  Xdel[K]   =  Dly 
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1500  NEXT  K 

1506  IF  Ntankl=l  THEN  EXIT  TO  1520 

1510       NEXT  Jr 

1520       IF  Lalt=Ialt  THEN  1900 
1530       IF  Jfirst=0  THEN  1800 
1540      REM  PRINT 
1800       Jfirst  =  1 

1810  MOVE  0,95-Xpace"10    |    SET  TEXT  COLOR  13 

1811  TEXT  USING  "D , A, 2D , 9 ( 4D) " : lalt , "x" , Ntankl , Lcmax 

1812  MOVE  0,90-Xpace"10    |   SET  TEXT  COLOR  3   |   TEXT  USING  " 4X , 9 (X, D . D) " :XdeI 
1820       Lalt  =  lalt 

1830       GOTO  1950 

1900  MOVE  0,95-XpaceAlO    |    SET  TEXT  COLOR  13 

1901  TEXT  USING  "2X, 2D , 9 ( AD ) " : Ntankl , Lcmax 

1902  MOVE  0  ,  90-Xpace'''10    |    SET  TEXT  COLOR  3    |    TEXT  USING  "  4X ,  9  (X ,  D  .  D )  "  :  Xdel 
1950      Xpace  =  Xpace+1 

2000  NEXT  Jj 

2010  REM  END  OOF  LTHIS  ALTITUDE  AND  VELOCITY 

4999  END 

5000  SELECT  CASE  K 

5010   CASE  1;4;5;8;9:12:13;16 

5020       GOTO  6000 

5030  CASE  ELSE 

5040      GOTO  7000 

5050  END  SELECT 

6000  Nintvx  =  Nintvl 

6010  Lspx  =  Lspl 

6020  Singx  =  Singl 

6030  SET  DIALOG  COLOR  14    |   PRINT  "#1  JJ=";Jj    |    SET  DIALOG  COLOR  3 

6040  RETURN 

7000  Nintvx  =  Nintv2 

7010  Lspx  =  Lsp2 

7020  Singx  =  Sing2 

7030  SET  DIALOG  COLOR  4    |   PRINT  "#2  jj=";Jj    |    SET  DIALOG  COLOR  3 
7040  RETURN 

LLDELAY.TAB 

This  program  stores  prints  and  stores  the  data  for  the  detailed  line  length/opening  delay  tables,  located  near  the  end  of 
the  guideline. 

1  GOTO  3000 

4  GOTO  100 

8  GOTO  700 
12  GOTO  4000 
16  GOTO  800 
20  GOTO  1000 
24  GOTO  5000 
28  GOTO  6000 
100  CLEAR 

110  DIM  Dl[50,8] ,D2[50,8] ,Z1[50] ,Z2[50] 
120  Dl  =  0 
130  D3  =  0 
140  D2  =  0 
150  D4  =  0 

160  DATA  0.5,1,2,3,4,6,8,10 
170  DIM  Z[8] 
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190  FOR  I  =  1  TO  8 
200       READ  Z[I] 


210 

NEXT 

I 

220 

DATA 

1  . 

1  .  1 

•  2  , 

,  1  ■ 

3,1, 

.6,1 

.  1 

230 

DATA 

1  . 

2  .  1 

.  3 

1  . 

6,1. 

1,1 

.  2 

24  0 

DATA 

1 . 

3,1 

.6, 

,  1 . 

1,1. 

2,1 

.  3 

250 

DATA 

1  . 

6,1 

■  1  , 

,  1  • 

2,1. 

.3,1 

.  6 

260 

DATA 

2. 

1,2 

■  2, 

,  2  . 

3,2. 

.1,2 

■  2  , 

2  . 

3  ,  2 

.1,2. 

2,2. 

3,2. 

1 

270 

DATA 

2. 

2,2 

.3, 

,  2  . 

1,2. 

2,2 

■  3  , 

3  . 

1 ,  3 

.2,3. 

1,3. 

2,3. 

1 

280 

DATA 

3  . 

2,3 

•  1 , 

,  3  . 

2,3 

.1,3 

■  2  , 

6  . 

1 ,  6 

.1,6. 

1,6. 

1,6. 

1 

310 

RESTORE 

220 

320 

FOR  I 

1 

TO 

50 

330 

READ 

Z1[I] 

340 

NEXT 

I 

350 

DATA 

1  , 

1,1,1. 

.2, 

2,2, 

,2,3 

,3, 

3, 

3,4 

,4,4, 

4,5, 

5,5, 

5  , 

360 

DATA 

2, 

2,2 

,  3  , 

.3, 

3,4. 

,4,4 

,5, 

5, 

5,1 

,1,2, 

2,3, 

3,4, 

4 

370 

DATA 

5, 

5,1 

,2, 

.3, 

4,5 

380 

RESTORE 

350 

390 

FOR  I 

■  1 

TO 

50 

4  00       READ  V3 

410       Z2[I]   =  V3"100 

420  NEXT  I 

440  FOR  I  =  1  TO  8 

450       FOR  J  =  1  TO  50 

460  PRINT  Z1[J];"     ":Z[I]:"   gpc  ";Z2[J];"   feet  "; 

470  INPUT  D1[J,I] 

480       NEXT  J 

4  90  NEXT  I 

500  PRINT  "DELAYS" 

510  FOR  I  =  1  TO  8 

520       FOR  J  =  1  TO  50 

530  PRINT  Z1[J]:"     ";Z[I]:"   gpc   ";Z2[J]:"  FEET  "; 

540  INPUT  D2[J,I] 

550      NEXT  J 
560  NEXT  I 

570  PRINT  "DATA  ENTERED" 

580  PRINT  "DO  YOU  WANT  TO  STORE  THE  DATA  ON  A  FILE?  "; 

590  INPUT  Z$ 

600  0$  =  SEGS(Z$,1,1) 

610  IF  0$="Y"  THEN  700 

620  END 

700  PRINT  "ENTER  TITLE-"; 

710  INPUT  Q$ 

720  PRINT  "WHAT  FILE?  "j 

730  INPUT  F$ 

740  OPEN  #1 :F$, "F" 

750  PRINT  #1:F$,Q$,D1,D2 

760  PRINT   "DATA  FOR:    ";QS;"  ON  HARD  DISC  ON  FILE  ":F$ 
770  CLOSE 
780  END 

800  PRINT   "WHAT  FILE?  "; 

810  INPUT  FS 

811  DELETE  Dl , D2 , Z 1 , D3 , D4 , Z2 

812  DIM  Dl [50,8] ,D2[50,8] 
820  OPEN  #1:F$,"R" 

830  INPUT  #1:F$,Q$,D1,D2 
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840 

PRINT 

"DATA  ENTERED  FROM  FILE  ";F$;" 

"  ;QS 

850 

END 

1000 

OPEN 

ff41  :  "printf  lie"  ,  "f  " 

1010 

PRINT 

#41 :Q$ 

1020 

Z4  = 

0 

1030 

C  =  1 

1040 

S$  = 

1050 

T$  = 

1  1" 

1                      1  1 

1060 

DATA 

1,2,3,6,1,2,3,5,1,2,3,6,1,2,3,6, 

1,2,3,6 

1065 

DATA 

T        001       Oil        OTI  O'iT 

1,2,3,1,2,3,1,2,3,1,2,3,1,2,3 

1070 

DATA 

1,2,1,2,1,2,1,2,1,2,1,1,1,1,1 

1090 

P$  = 

1100 

RESTORE  1060 

1110 

DELETE  Z3 

1120 

DIM  Z3[50] 

1130 

FOR  I 

=  1  TO  50 

llAO 

READ  Z3[I] 

1150 

NEXT 

I 

1160 

PRINT 

#41  :  " 

1170 

FOR  I 

=  26  TO  129 

1180 

PRINT  #41 : "_" : 

1190 

NEXT 

I 

1200 

PRINT 

#41 : "_" 

1210 

PRINT 

#41  :  " 

1" 

1220 

PRINT 

#41 : "                     LEVEL   1  1 

1230 

PRINT 

#41:"       LEVEL  2     1       LEVEL  3  | 

LEVEL  4 " ; 

1240 

PRINT 

#41:  " 

1  II 
1 

1250 

PRINT 

#41 : S$;T$;P$ 

1260 

PRINT 

#41:  " 

1  11 

1270 

PRINT 

#41:  " 

COVERAGE ( GALLONS / " 

1280 

PRINT 

#41 : "100   SQ.FT. ) 

1" 

1290 

PRINT 

#41 :SS;T$;P$ 

1300 

PRINT 

#41  :  " 

1  "  ' 

1310 

PRINT 

#41   USING  2790 : 0 . 5, " 1 " ,1 , " 1 " ,2, 

ill,       ti  I  It       II  1            1  II       II  1 

1    >  3 ,    1    , 4 ,    1    , 6 ,    1    , 8 ,  1 

1320 

PRINT 

#41:"        10.0  1" 

1330 

S$  = 

1  1 

1          -  1 

1340 

PRINT 

#41 :S$;T$:P$ 

1350 

PRINT 

#41 : " |compartments|drop|no.of 

1     MAX   1  DELAY |     MAX   | DELAY | " 

1360 

PRINT 

#41:"      MAX    1  DELAY       MAX    | DELAY 

MAX     DELAY |     MAX   | DELAY | " ; 

1370 

PRINT 

#41:"     MAX    1  DELAY       MAX    | DELAY 

1380 

PRINT 

#41:"|      RELEASED      |    HT    | PATTERNS | LENGTH   ( SEC ) | LENGTH | ( SEC ) | " 

13  90 

PRINT 

#41 : "LENGTH   (SEC) | LENGTH | (SEC) | 

LENGTH   (SEC) | LENGTH | (SEC) | " ; 

14  00 

PRINT 

#41 : "LENGTH| (SEC) |lENGTH| (SEC) 1 

1  H  1  u 

PRINT 

#41  :  "  1      AT  A  TIME    |    FT  | 

1    (FT)    1            1    (FT)  1 

1420 

PRINT 

#41:"    (FT)    1            1    (FT)  | 

(FT)    1               (FT)    1  1": 

1430 

PRINT 

#41:"    (FT)    1            1    (FT)    1  1 

1440 

S$  = 

1450 

T$  = 

14  60 

P$  = 

14  70 

DELETE  V$,U$ 

1480 

DIM  VS(132) ,US(132) 

1490 

U$  = 

S$  &  T$ 

1500 

U$  = 

US  &  P$ 

1510 

PRINT 

#41 :U$ 

1520 
1530 

SS  = 
V$  = 

"1                              |.-._|  j... 

SS   5.  T$ 
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1540  V$  =  V$   S.  P$ 
1550  FOR  I  =  1  TO  50 

1560        IF   1=3  OR   1=27  OR  1=47  OR  1=61   OR  1=71   THEN  2380 

1570       IF  1=8  OR  1=31  OR  1=49  OR  1=63  OR  1=72  THEN  2400 

1580       IF  1=13  THEN  2440 

1590       IF  1=35  THEN  2460 

1600       IF  1=53  THEN  2480 

1610       IF  1=65  THEN  2500 

1615       IF  1=73  THEN  2520 

1620  IF  1=18  OR  1=39  OR  1=56  OR  1=67  OR  1=74  THEN  2640 
1630       IF  1=23  OR  1=43  OR  1=59  OR  1=69  OR  1=75  THEN  2760 

2070  IF  Z3(I]<>12  THEN  2075 

2071  PRINT  #41:"|  |  |        ";Z3[I];"  "; 

2072  GOTO  2080 

2075       PRINT  #41:"|  |  |  ":Z3[I]:"  "; 

2080  FOR  J  =  1  TO  8 

2081  IF  Q$<>"BLANK"  THEN  2090 

2082  IF  D1[I,J]=0  THEN  2095 

2090  PRINT  #41  USING  "A, 6D , A, 3D . D , S" : " | " , Dl [ I , J] , " | " , D2 [ I , J] 

2093  GOTO  2100 

2095  PRINT  #41  USING  "A , 6X , A , 5X , S" : " | " , " | " 

2100       NEXT  J 

2110       PRINT  #41 : " I " 

2115       REM  NEXT  IS  LONG  LINE 

2120       IF  1=25  OR  1=45  OR  1=60  OR  1=70  OR  1=75  THEN  2830 
2125      REM  NEXT  TWO  ARE  SHORT  LINES 

2130       IF  1=5  OR  1=10  OR  1=15  OR  1=20  OR  1=29  OR  1=33  OR  1=37  THEN  2810 

2140  IF  1=41  OR  1=48  OR  1=51  OR  1=54  OR  1=57  OR  1=62  THEN  2810 

2141  IF  1=64  OR  1=66  OR  1=68  THEN  2810 
2340  NEXT  I 

2350  Z4  =  1 

2360  CALL   "Ipr  printfile" 

2361  END 
2370  NEXT  I 

2380  PRINT  #41:"|  |    100|  ";Z3[I];"  ": 

2390  GOTO  2080 

2400  PRINT  #41:"|  |    200 1  ";Z3[I];"  "; 

2410  GOTO  2080 

2420  ON  C  GOTO  2440,2460,2480,2500,2520,2540,2560,2580,2580,2580 
2440  PRINT  #41: "I  ONE  |    300 |  ":Z3[I];"  ": 

2450  GOTO  2590 

2460   PRINT  #41:"  I  TWO  |    300  |  ";Z3[I];"  ": 

2470  GOTO  2590 

2480  PRINT  #41:" I  FOUR  |    300 |  ";Z3[I];"  "; 

2490  GOTO  2590 

2500  PRINT  #41: "I  SIX  |    300 |  ";Z3[I];" 

2510  GOTO  2590 

2520  PRINT  #41: "I        TWELVE        |    300 |  ":Z3[I];"  "; 

2530  GOTO  2590 
2580  GOTO  2860 
2590  C  =  C+1 
2600  GOTO  2080 

2610  ON  C-1  GOTO  2640,2660,2680,2700,2620,2720,2740,2620,2620,2620 
2620  PRINT  #41:"|    (ALL  COMB)    |    400|  ";Z3[I]:"  "; 

2630  GOTO  2080 

2640  PRINT  #41:"|  |    400|  ":Z3[I];"  "; 

2650  GOTO  2080 
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2660 
2670 
2680 
2690 
2700 
2710 
2720 
2730 
2740 
2750 
2760 
2770 
2780 
2781 
2790 
2800 
2810 
2820 
2830 
28A0 
2850 
2860 
2870 
2880 
2890 
2900 
2910 
2920 
2930 
3000 
3010 
3020 
3030 
30A0 
3050 
3060 
3065 
3066 
3080 
3090 
3100 
3110 
3120 
3130 
3140 
3150 
3160 
AOOO 
4010 
4020 
4021 
4030 
4040 
4042 
4044 
4050 
4052 


400 


400 


400 


400 


500 


PRINT  #41  :  "  I  (I)  I  400 

GOTO  2080 

PRINT  #41 : "  I  (0,1) 
GOTO  2080 

PRINT  #41  :  "  I  (1,1) 
GOTO  2080 

PRINT  #41 : " I  (O+I) 
GOTO  2080 

PRINT  #41 : " I  (I+I) 
GOTO  2080 
PRINT  #41 : " I 
GOTO  2080 

CALL  "Ipr  printfile" 
END 

IMAGE  7(4X,2D.D,4X.A) ,S 
IMAGE  A,5D,A,D.D,S 
PRINT  #41 : V$ 
GOTO  2340 
PRINT  #41 :U$ 
IF  24=1  THEN  2360 
GOTO  2340 

ON  C-7  GOTO  2880,2900,2920 
GOTO  2590 

TWO 


';Z3[I]  ; 
':Z3[I]  ; 
':Z3[I]  ; 
':Z3[I]  ; 
':Z3[I]  : 
':Z3[I) ; 


300 


':Z3[I] 
':Z3[I] 
':Z3[I] 


PRINT  #41 ; " I 
GOTO  2590 

PRINT  #41:" I        THREE  |    300 | 

GOTO  2590 

PRINT  #41  :  "  I       SIX  I    300 1 

GOTO  2590 
CLEAR 

PRINT   "USER  DEFINABLE  KEYS:" 
1 .     ENTER  DATA" 
STORE  DATA" 
MAKE  CORRECTION" 
RETRIVE  DATA" 

PRINT  LINE  LENGTH-DELAY  TABLE" 
SUBTRACT  WATER  FROM  RETARDANT  &  PRINT  TABLE" 
PRINT  TABLE  W/O  LINES   (FOR  PRINTERS)" 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
END 

FOR  I  =  1  TO  30 
FOR  J  =  1  TO  8 

D1[I,J]    =  100>''I+J 
D2[I,J]   =  1+J/lO 
NEXT  J 
NEXT  I 
GOTO  570 
END 
CLEAR 

PRINT  "ENTER  THE  ROW  NUMBER  THEN  THE  GPC 
INPUT  C1,C2 

IF  C1=0  AND  C2=0  THEN  580 
IF  C2=6  THEN  4070 
IF  C2<>8  THEN  4050 
C2  =  7 
GOTO  4070 

IF  C2O10  THEN  4060 
C2  =  8 
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4054 
A060 
A070 
4080 
4090 
4095 
5000 
5005 
5010 
5020 
5030 
5040 
5050 
5060 
5070 
5080 
5090 
5100 
5110 
5120 
5125 
5130 
6000 
6010 
6020 
6030 
6040 
6050 
6060 
6070 
6080 
6100 
6110 
6120 
6130 
6140 
6150 
6160 
6170 
6210 
6220 
6230 
6240 
6250 
6260 
6270 
6280 
6290 
6300 
6310 
6320 
6330 
6350 
6360 


GOTO  4070 

C2  =  INT(C2+1) 

PRINT  "THE  PRESENT  VALUES  ARE    :    " : Dl [ CI , C2 ] : "  ";D2[C1,C2] 

PRINT  "ENTER  NEW  VALUES  ": 

INPUT  Dl [CI ,C2] ,D2[C1 ,C2] 

GOTO  4010 

CLEAR 

DELETE  W1,R1,D1,D2 

DIM  Wl[50,8) ,R1[50,8] ,D2[50,8] ,D1[50,8] 

PRINT  "WHAT  FILE  FOR  WATER?  "; 

INPUT  F$ 

OPEN  #1 :F$, "R" 

READ  #1:F$,W$,W1 

PRINT  "WHAT  FILE  FOR  RETARDANT?  "; 

INPUT  G$ 

OPEN  #2:G$, "R" 

READ  #2:G$,R$,R1 

Dl  =  Rl-Wl 

D2  =  0 

Q$  =  "RETARDANT-WATER" 

CLOSE 

GOTO  1000 

REM 

OPEN  #41 : "PRINTFILE" , "W" 

PRINT  #41: 

Z4  =  0 

C  =  1 

S$  =  "" 

T$  =  "" 

DATA  1,2,3,6,1,2,3,6,1,2,3,6,1,2,3,6,1,2,3,6,1,2,3,1,2,3,1,2,3 
DATA  1,2,3,1,2,3,1,2,1,2,1,2,1,2,1,2,1,1,1,1,1 
P$  =  "" 
RESTORE  6070 
DELETE  Z3 
DIM  Z3[50] 
FOR  I  =  1  TO  50 
READ  Z3[I] 


NEXT  I 

PRINT 

#41 

PRINT 

#41 

PRINT 

#41 

11 

PRINT 

#41 

LEVEL  1  " 

PRINT 

#41 

LEVEL  2             LEVEL  3 

PRINT 

#41 

tt 

PRINT 

#41 

PRINT 

#41 

PRINT 

#41 

PRINT 

#41 

"100  SQ.FT.) 

PRINT 

#41 

PRINT 

#41 

PRINT 

#41 

USING  7800:0.5,"  ",1,"  ",2," 

"  ,  3  ,  " 

PRINT 

#41 

10.0 

PRINT 

#41 

PRINT 

#41 

"  COMPARTMENTS  DROP  NO. OF 

MAX 

LEVEL  4"; 


COVERAGE ( GALLONS / 


MAX  DELAY 
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6370 
6380 
6390 
6A00 
6A10 
6A20 
6430 
6440 
6500 
6520 
6540 
6560 
6562 
6570 
6580 
6590 
6600 
6610 
6620 
6630 
6640 
6650 
6660 
6670 
6680 
6690 
6700 
6710 
6720 
6730 
6740 
6750 
6760 
6770 
6780 
6790 
6800 
6810 
6820 
6830 
6840 
6850 
6860 
6870 
6880 
6890 
6900 
6910 
6920 
6930 
6940 
6950 
6960 
6990 
7000 
7110 
7120 


MAX  DELAY 
MAX  DELAY 
RELEASED 


MAX     DELAY       MAX  DELAY 


MAX  DELAY 
MAX     DELAY  " 

HT     PATTERNS  LENGTH   (SEC)   LENGTH  (SEC) 


LENGTH  (SEC)  LENGTH  (SEC)  LENGTH  (SEC)  LENGTH  (SEC) 
LENGTH   (SEC)   LENGTH   (SEC)  " 

AT  A  TIME       FT  (FT) 

(FT)  (FT)  (FT) 

(FT)  (FT) 


(FT) 
(FT) 


PRINT  #41 : " 
PRINT  #41 : " 
PRINT  #41 :  " 
PRINT  #41 : " 
PRINT  #41 : " 
PRINT  #41 : " 
PRINT  #41 : " 
PRINT  #41 : " 
U$  =  "" 
PRINT  #41: 
V$  = 

FOR  I  =  1  TO  50 

IF  1=3  OR  1=36  OR  1=46  OR  1=22  THEN  6800 
IF  1=6  OR  1=25  OR  1  =  38  OR  1=47  THEN  6820 
IF  1=10  THEN  6850 
IF  1=28  THEN  6870 
IF  1=40  THEN  6890 
IF  1=48  THEN  6910 

IF  1=14  OR  1=31  OR  1=42  OR  1=49  THEN  6990 
IF  1=18  OR  1=34  OR  1=44  OR  1=50  THEN  7110 
PRINT  #41:"  ";Z3[I] ;"  "; 

FOR  J  =  1  TO  8 

IF  Q$<> "BLANK"   THEN  6680 

IF  D1[I,J]=0  THEN  6700 

PRINT  #41  USING  "A , 6D , A , 3D . D , S" : "   ",D1[I,J],"  ",D2[I,J] 
GOTO  6710 

PRINT  #41  USING  " A , 6X , A , 5X , S" : "   ","  " 
NEXT  J 

PRINT  #41  :  "  " 

IF  1=20  OR  1=35  OR  1=45  OR  1=50  THEN  7180 

IF  1=4  OR  1=8  OR  1=12  OR  1=16  OR  1=23  OR  1=26  OR  1=29  THEN  7160 
IF  1=32  OR  1=37  OR  1=39  OR  1=41  OR  1=43  THEN  7160 

NEXT  I 

Z4  =  1 

END 

NEXT  I 

PRINT  #41:"  100  ";Z3[I] ;"  "; 

GOTO  6650 

PRINT  #41  :  "  200  " : Z3[I]  ; "       "  : 

GOTO  6650 

ON  C  GOTO  6850,6870,6890,6910,2520,2540,2560.6930,6930,6930 


ONE 


300 


TWO 


THREE 


SIX 


300 


300 


300 


;Z3[I] ;" 
;Z3[I]  ;" 
;Z3[I] ;" 
;Z3[I]  ;" 


PRINT  #41 : " 
GOTO  6940 
PRINT  #41 : " 
GOTO  6940 
PRINT  #41 : " 
GOTO  694  0 
PRINT  #41 : " 
GOTO  6940 
GOTO  7210 
C  =  C+1 
GOTO  6650 

ON  C-1  GOTO  6990,7010,7030,7050,6970,7070,7090,6970,6970,6970 
PRINT  #41 : "  400  " ; Z3[I] ; "        " ; 

GOTO  6650 

PRINT  #41 : "  500  " ; Z3[I] ; "        " ; 

GOTO  6650 
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7130  CALL  "Ipr  PRINTFILE" 

7131  END 

71A0  IMAGE  7(4X,2D.D,4X,A) ,S 
7150  IMAGE  A,5D,A,D.D,S 
7160  PRINT  #A1:VS 
7170  GOTO  6760 
7180  PRINT  #41 :U$ 
7190  IF  Z4=l  THEN  6780 
7200  GOTO  6760 

7210  ON  C-7  GOTO  7230,7250,7270 
7220  GOTO  6940 

7230  PRINT  #41      I  TWO         |    300 |  ";Z3[I];"  "; 

7240  GOTO  6940 

7250  PRINT  #41; "I       THREE         |    300|         ":Z3[I];"  "; 
7260  GOTO  6940 

7270  PRINT  #41: "I        SIX  |    300 |  ";23[I];"  "; 

7280  GOTO  6940 

7800  IMAGE  7 ( 4X , 2D . D , 4X , A) , S 

LDRWIDTH 

This  program  requires  input  of  data  from  the  detailed  pattern  section  of  the  PATSIM  output  and  calculates  the  line 
length,  downrange  position  and  width  for  1,  2,  3,  4,  and  5  gpc  to  be  used  in  the  FOOTPRINTS.  Because  the  detailed  pat- 
tern section  of  the  PATSIM  output  is  bilaterally  symetrical  about  rows  14-15,  only  data  from  the  left  side  are  entered.  The 
program  uses  a  linear  interpolation  for  its  calculations  and  starts  at  1  gpc.  Data  more  than  one  space  to  the  left  or  above 
and  below  1-gpc  values  are  not  used  and  need  not  be  entered. 

8  GOTO  50000 
20  GOTO  50000 

100  REM  CALCULATE  THE  AVERAGE  WIDTH,  DOWNRANGE  AND  LENGHT  FROM  PATSIM 
105  K7  =  0 

110  PRINT  "ENTER  THE  ID  FOR  THIS  DROP--"; 
120  INPUT  1$ 

130  DIM  D[40,14] ,Ans[40,5] ,Start[2] ,Finish[2] 
140  Ans  =  0 
150  D  =  0 
160  Start  =  0 
170  Finish  =  0 

180  PRINT  "DO  YOU  WISH  TO  RECALL  A  FILE?  IF  YES  ENTER  NAME  ELSE  CR  "; 

190  INPUT  F$ 

200  IF  F$=""  THEN  250 

205  DIM  X$(300) 

210  OPEN  #1;F$,"R",X$ 

220  READ  #1:I$,D 

230  GOTO  4  90 

240  PRINT  "ENTER  -1  TO  TERMINATE  A  COLUMN" 
250  FOR  I  =  1  TO  7 
260       FOR  J  =  1  TO  4  0 

270  PRINT  "WHAT  IS  THE  VALUE  No.": J;"  FOR  COL.    " : 1+7 : "  ?  "; 

280  INPUT  Test 

290  IF  Test<0  THEN  EXIT  TO  320 

300  D[J,I]   =  Test 

310       NEXT  J 

320  NEXT  I 

330  FOR  I  =  1  TO  40 

340       D[I,8]   =  D[I,7] 

350       D[I,91   =  DrT,61 
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360  D[I,10]  =  D[I,5] 
370  D[I,11]  =  D[I,4] 
380  D[I,12]  =  D[I,3] 
390  D(I,13]  =  D[I,2] 
400  D[I,14]  =  D[I,1] 
410  NEXT  I 

420  PRINT  "WHAT  FILE  FOR  STORAGE?  "i 

430  INPUT  S$ 

450  OPEN  #1 : S$, "F" ,XS 

460  WRITE  #1 : I$,D 

470  CLOSE 

480  PRINT  "DATA  SAVED  ON  DISC  FILE  ":S$ 

490  Count  =  1 

500  Start  =  0 

520  FOR  Gpc  =  1  TO  5 

530       FOR  Row  =  1  TO  14 

540  IF  Start[l]>0  THEN  590 

550  IF  D[ Count, Row] <Gpc  THEN  630 

560  IF  Start [1]<>0  THEN  630 

570  Start [1]   =  Row 

580  Start [2]  =  D[ Count, Row] 

590  IF  D [ Count , Row] >=Gpc  THEN  630 

600  Finish[l]   =  Row-1 

610  Finish[2]   =  D [ Count , FLnish[ 1 ] ] 

620  EXIT  TO  640 

630       NEXT  Row 

640       IF  Start[l]=0  THEN  790 

650       W  =   (FinLsh[l]-Start[l]  )''15 

660       IF  Start[l]=l  THEN  710 

670       Wtop  =   (Start[2]-Gpc)/ (Start[2]-D[Count,Start[l]-l])*15 
680       Wbot  =  (Flnish[2] -Gpc) / (FLnLsh[2] -D[Count,FLnish[l]+l] )'*15 
690       W  =  W+Wtop+Wbot 
700      GOTO  740 

710       Wtop  =  (Start[2]-Gpc)/Start[2]^''15 
720       Wbot  =  (Finish[2]-Gpc)/FLnish[2]>''15 
730       W  =  W+Wtop+Wbot 
740       Ans [Count , Gpc]   =  W 

750       REM  print    "*";gpc;"    ";w;"    ";Start(l);"    "; finish( 1 );"•*" ;wtop ; "  ";wbot 

760       Start  =  0 

770       Row  =  0 

780  NEXT  Gpc 

790  Count  =  Count+1 

800  IF  Count=41  THEN  820 

810  GOTO  520 

820  DIM  Som[5] 

830  Som  =  0 

840  FOR  Gpc  =  1  TO  5 

850       FOR  I  =  1  TO  40 

860  Som[Gpc]  =  Som [Gpc ] +Ans [ I , Gpc ] 

870       NEXT  I 

880  NEXT  Gpc 

890  CLEAR 
1000  DIM  Ave[5] 
1010  De  =  41 
1020  Ave  =  0 
1030  Countl  =  0 
1040  Count2  =  0 
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1050 

Count 3  =  0 

1060 

Count4  =  0 

1070 

Counts  =  0 

1080 

P$  =  "5(1X,4D.D,2X)" 

1090 

FOR  I  =  1  TO  40 

1100 

IF  Ans[I,l]=0  THEN  1200 

1110 

Countl  =  Countl+1 

1120 

IF  Ans[I,2]=0  THEN  1200 

1130 

Count2  =  Count2+l 

1140 

IF  Ans[I,3]=0  THEN  1200 

1150 

Counts  =  Count3+l 

1160 

IF  Ans[I,4]=0  THEN  1200 

1170 

Count4  =  Count4+l 

1180 

IF  Ans[I,5]=0  THEN  1200 

1190 

Counts  =  Count5+l 

1200 

NEXT  I 

1210 

Ave[l]  =  Som[l] /Countl 

1220 

IF  Som[2]=0  THEN  1380 

1230 

Ave [2]  =  Som[2] /Count2 

12A0 

IF  Som[3]=0  THEN  1380 

1250 

Ave[3]  =  Som[ 3] /Counts 

1260 

IF  Som[4]=0  THEN  1380 

1270 

Ave[4]  =  Som[A] /CountA 

1280 

IF  Som[5]=0  THEN  1380 

1290 

Ave[5]  =  Som[ 5] /Counts 

1380 

GOTO  10000 

10000 

DIM  D14[40] ,L[5,8] ,L9[5] ,M2[8] 

10010 

D14  =  0 

10020 

FOR  I  =  1  TO  40 

10030 

D14[I]   =  D[I,7] 

10040 

NEXT  I 

10050 

L  =  0 

10060 

L9  =  0 

10070 

M2  =  0 

10080 

FOR  J  =  1  TO  5 

10090 

C  =  1 

10100 

K  =  1 

10110 

IF  D14(K]>=J  THEN  10150 

10120 

K  =  K+1 

10130 

IF  K=41  THEN  10370 

10140 

GOTO  10110 

10150 

A  =  K"30 

10160 

0  =  D14[K] 

10170 

IF  Kol  THEN  10200 

10180 

U  =  0 

10190 

GOTO  10210 

10200 

U  =  D14[K-1] 

10210 

A  =  A-30/ (0-U)'"(0-J) 

10220 

L[J,C]   =  A 

10230 

C  =  C+1 

10240 

K  =  K+1 

10250 

IF  K=41  THEN  10310 

10260 

IF  D14[K]<J  THEN  10300 

10270 

K  =  K+1 

10280 

IF  K=41  THEN  10370 

10290 

GOTO  10260 
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10300  B  =  (K-D'-SO 

10310  0  =  D14[K-1] 

10320  U  =  D14[K] 

10330  B  =  B+30/ (O-U)A(O-J) 

10340  L[J,C]   =  B 

10350  C  =  C+1 

10360  GOTO  10120 

10370  L9[J]  =  L9[J]+(L[ J,2] -L[J,1] )+(L[J,4] -L[ J, 3] )+(L[J,6]-L[J,5] ) 

10380  L9[J]  =  L9[J]+L[J,8]-L[J,7] 

10390  NEXT  J 

20000  REM  PRINT  SECTION 

20001  Far  =  LEN(XS) 

20002  XS  =  SEG$(X$,200,Far-200) 

20003  X$  =  EDIT$(X$, "b") 
20020  P$  =  "2X,2(4D,X) ,S" 
20022  DELETE  FILE  "printfile" 
20025  OPEN  #2 : "printfile" , "f" 

20030  PRINT  #2:1$:"      " ; X$ 

20031  PRINT  #2: 

20032  PRINT  #2: 

20040  PRINT    #2  :  "^^'"'1*^1''^''' A'>A'''AAAAi''l'"*fi'fVfi'fA'*fT'cycVfyryryfVfyfycyfVfiVVfVfVfVfVfVfyffrVrAVn'rT'rffyrVry(rfVf>'("  . 

20050  PRINT  #2:  "'■<>■"■"■""■"■"■""'*" 

20060  PRINT  #2:"       ":"  ": 

20070  PRINT  #2:"         1  2  3  4  5" 

20080  PRINT  #2  :  "  Avvft,vftVcVrVrDRVcVcv,D+L'''''"'''''DR''"''>''D+L''"""''DR''"'"''D+L'-''-"'"''DR'""''D+L^  ; 

20090  PRINT  #2:  "''DR'"*''D+L''"-'" 

20100  R$  =  "4X,4D,4X,S" 

20110  FOR  J  =  1  TO  5 

20120  ON  J  GOTO  20130,20270,20400,20530,20660 

20130  PRINT  #2: "LI  "; 

20140  FOR  II  =  1  TO  5 

20150  PRINT  #2  USING  P$ : L [ II , 1 ] , L [ II , 2] 

20160  NEXT  II 

20170  PRINT  #2: 

20180  PRINT  #2: "LEN      " ; 

20190  FOR  Count  =  1  TO  5 

20200  PRINT  #2  USING  R$ : L [ Count , 2 ] -L [ Count , 1 ] 

20210  NEXT  Count 

20220  PRINT  #2: 

20230  PRINT  #2: "  

20240  PRINT  #2:  "  " 

20250  NEXT  J 

20260  GOTO  20660 

20270  PRINT  #2: "L2       "  ; 

20280  FOR  II  =  1  TO  5 

20290  PRINT  #2  USING  P$ : L [ I 1 , 3 ] , L [ II , 4 ] 

20300  NEXT  II 

20310  PRINT  #2: 

20320  PRINT  #2  :  "LEN     "  ; 

20330  FOR  Count  =  1  TO  5 

20340  PRINT  #2  USING  RS : L [ Count , 4 ] -L [ Count , 3 ] 

20350  NEXT  Count 

20360  PRINT  #2: 

20370  PRINT  #2:  "  " 

20380  PRINT  #2:  "  " 

20390  GOTO  20250 

20400  PRINT  #2:  "L3       "  ; 
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20410  FOR  II  =  1  TO  5 

20420       PRINT  #2  USING  P$ ; L [ II , 5 ] , L [ II , 6 ] 

20430  NEXT  II 

20440  PRINT  #2: 

20450  PRINT  #2: "LEN  ": 

20460  FOR  Count  =  1  TO  5 

20470       PRINT  #2  USING  R$ : L [ Count , 6 ] -L [ Count , 5 ] 
20480  NEXT  Count 
20490  PRINT  #2: 

20500  PRINT  #2: "  

20510  PRINT  #2:  "  " 

20520  GOTO  20250 

20530  PRINT  #2:"L4  "; 

20540  FOR  II  =  1  TO  5 

20550      PRINT  #2  USING  P$ : L [II , 7 ] ,L [ II , 8] 

20560  NEXT  II 

20570  PRINT  #2: 

20580  PRINT  #2: "LEN  "; 

20590  FOR  Count  =  1  TO  5 

20600       PRINT  #2  USING  R$ : L [Count , 8] -L [Count , 7 ] 
20610  NEXT  Count 
20620  PRINT  #2: 

20630  PRINT  #2: "  

20640  PRINT  #2:  "  " 

20650  GOTO  20250 

20660  Q$  =  "4X,3D.D,3X,S" 

20670  PRINT  #2: "WIDTH  "; 

20680  FOR  I  =  1  TO  5 

20690       PRINT  #2  USING  Q$:Ave[I] 

20700  NEXT  I 

20710  PRINT  #2: 

30000  PRINT  #2: "ROW  8  9         10         11         12         13  14" 

30010  FOR  I  =  1  TO  22 

30020       PRINT  #2  USING  "2D,3X,S":I 

30030       FOR  J  =  1  TO  7 

30040  PRINT  #2  USING  "X , 2D . D , X , S" : D [ I , J ] 

30050       NEXT  J 

30060       PRINT  #2: 

30070  NEXT  I 

30080  CLOSE  #2 

30085  CALL  "Ipr  printfile" 

30090  GOTO  100 

50000  REM  correct  data 

50005  K7  =  1 

50010  PRINT  "WHAT  FILE?  "; 
50020  INPUT  F$ 

50030  DIM  D[40,14] ,Ans[40,5] ,Start[2] ,Finish[2] 

50031  Ans  =  0 

50032  Start  =  0 

50033  D  =  0 

50034  Finish  =  0 

50040  OPEN  #1:F$,"R",XS 
50050  READ  #1:I$,D 
50055  CLOSE 
50060  CLEAR 

50070  PRINT  "DATE  FOR:    ":!$:"   IN  MEMORY" 

50071  PRINT  "FNTER  P  FOR  PRINT  ELSE  CR" ; 
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50072  INPUT  Q$ 

50073  IF  Q$="P"  THEN  490 

50080  PRINT  "TO  EXIT  CORRECTION  FOR  THIS  FILE  ENTER  -1  FOR  COLUMN  AND  ROW" 
50090  PRINT  "ENTER  COLUMN  THEN  ROW  "; 
50100  INPUT  Column, Row 
50105  Column  =  Column-7 
50110  IF  Column<0  THEN  50150 

50120  PRINT  "PRESENT  VALUE:    "; D [Row, Column] ; "  NEW  VALUE?  "; 
50130  INPUT  D[Row, Column] 
501A0  GOTO  50090 

50150  REM  STORAGE  AND  FILLOUT  SECTION 


50160 

FOR  1=1 

TO  4  0 

50170 

D[I,8] 

=  D[I,7] 

50180 

D[I,9] 

=  D[I,6] 

50190 

D[I,10] 

=  D[I,5) 

50200 

D[I,11] 

=  D[I,4] 

50210 

D[I , 12] 

=  D[I,3] 

50220 

D[I , 13] 

=  D[I,2] 

50230 

D[I,14] 

=  D[I,1] 

50240  NEXT  I 

50250  PRINT  "PRESS  RETURN  TO  RESTORE  DATA  ON  THE  SAME  FILE  " ; 
50260  INPUT  Q$ 

50270  OPEN  #1:F$,"F",X$   !  ?? 
50280  WRITE  #1:I$,D 
50290  CLOSE 

50300  PRINT  "DO  YOU  WANT  TO  PRINT  THE  NEW  DATA?   (Y  OR  N  ) 

50310  INPUT  Q$ 

50320  IF  Q$="N"  THEN  50340 

50330  GOTO  490 

50340  END 

BESTSTRAT.DAT/BS  Plot 

The  line  length  versus  height  data  for  each  drop  type  is  smoothed  by  using  a  standard  simple  regression  that  deter- 
mines the  best  fit  of  the  data  to  one  of  eight  curve  forms,  Y  =  AX;Y  =  A+BX;  Y  =  Ae^^;  Y  =  ViA+BX);  Y  =  A+B/X; 
Y  =  A+BlogX;  Y  =  AX'^;  and  Y  =  X1{A+BX).  Y  =  drop  height  and  X  =  line  length.  These  equations  along  with  the  type  of 
drop  are  entered  into  BESTSTRAT.DAT  lines  1000  to  2000  and  the  data  are  plotted.  After  selection  of  the  type  of  drop  for 
each  volume,  the  remaining  equations  are  set  to  0  (lines  1000-2000)  and  a  final  plot  made. 

1  GO  TO  100      ! PRINT  MENU 

2  REM  WATER  LIKE  3  GPC  CONAIR  DC-6B  TANKER  454 
4  GO  TO  2010      ! CALCULATE  STRATIGIES  AND  PLOT 

8  GO  TO  3270      IREPLOT  WITH  PREVIOUSLY  ENTERED  PARAMETERS 
12  GO  TO  8000      ! PRINT  A  LABEL  ON  THE  4662  PLOTTER 

16  T9=0     !PLOT  CALCULATED  DATA 

17  GO  TO  2010 

20  GO  TO  9000     ! PRINT  ON  PLOTTER  XTxXR  LABEL 
24  GO  TO  11000      !DRAW  LEFT  ARROWHEAD 
28  GO  TO  10000      ! DRAW  RIGHT  ARROWHEAD 

44  APPEND  "BSPLOT.A"  AT  2000      ! APPEND  PLOT  PROGRAM 

45  V6=0 

46  GO  TO  4 
100  CLEAR 
102  V6=0 

110  PRINT  "1.   CALCULATE  DATA" 

120  PRINT  "2.   DRAW  AXES  WITH  PREVIOUSLY  ENTERED  PARAMETERS" 
130  PRINT  "3.    ENTER  AND  PRINT  A  LABEL" 

140  PRINT  "4.    INPUT  PARAMETERS  FOR  AND  DRAW  GRAPH  FROM  CALCULATED  DATA" 
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150  PRINT  "5.   PRINT  XTxXR  LABEL" 

160  PRINT  "6.   DRAW  LEFT  POINT  ARROWHEAD" 

170  PRINT  "7.   DRAW  RIGHT  POINT  ARROWHEAD" 

180  PRINT  "11.   APPEND  PLOT  PORTION  FRON  4907  DISC" 

280  END 

990    ! PROGRAM  EQUATIONS  AND  ID  STRINGS  FOR  STRATIGIES  X=0  FOR  NO  ENTRY 

1000  AS="1X1  " 

1010  X=0 

1020  RETURN 

1030  B$="1X2  " 

1040  X=0 

1050  RETURN 

1060  C$="1X4" 

1070  X=3614  .  14759953"'J"-0  .  599571655A13 
1080  RETURN 

1090  D$="1X12"      ! EQUATION  AT  1110  IS  BASIS  FOR  NO  CAPABILITY  AREA 
1095  IF  J>376  THEN  1251 

1100  X=l 496.  071864 37^''EXP( -0.00347487486045* J) 
1110  RETURN 
1120  ES="1X6" 

1130  X=1295. 28917901-189. 3 A9621044*LOG( J) 

1140  RETURN 

1150  F$="2X1" 

1160  X=0 

1170  RETURN 

1180  G$="2X2" 

1190  REM  X=779. 924741195-119. 17387092*LOG(J) 

1191  X=0 
1200  RETURN 
1210  H$="2X3" 

1220  REM  X=12637. 4777927'' J^-0. 741958999899 

1221  X=0 
1230  RETURN 
1240  I$="2X6" 

1250  RETURN 

1251  X=l/ (5.694  73125E-4+5.052039738E-6"J) 
1260  RETURN 

1270  J$="4X1" 

1280  REM  X=236.  666666667-0.  675*'J 

1281  X=0 
1290  RETURN 
1300  K$="4X2" 

1310  REM  X=587  .  447659418'''EXP(-0  .  00286385476043'''J) 

1311  X=0 
1320  RETURN 
1330  L$="4X3" 

1340  REM  X=873 . 435667466"EXP(-0 . 00244759118039"J) 

1341  X=0 
1350  RETURN 
1360  M$="6X1" 

1370  REM  X=864 . 121558318-135. 493350465''LOG(J) 

1371  X=0 
1380  RETURN 
1390  N$="6X2" 

1400  REM  X=751.760972274'''EXP(-0.0024352907181''»J) 

1401  X=0 
1410  RETURN 
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1420  0S="12X1" 

1A30  REM  X=1762. 81484802-292.  712597148'''LOG(J) 
1431  X=0 
1440  RETURN 

2000  REM  TARGET  LINE  FOR  APPENDING  PLOT  PORTION 

2000  REM  PLOT  PORTION  OF  BEST  STRATEGY  CHART 

2010  DELETE  A , B , C , D , El , E2 , E3 , E4 , E5 , E6 , E7 , E8 , E9 , Fl , F2 

2020  DIM  A(26) ,B(26) ,0(26) ,D(26) ,E1(26) ,E2(26) ,E3(26) 

2030  DIM  E4(26) ,E5(26) ,E6(26) ,E7(26) ,E8C26) ,E9(26) 

2040  DIM  Fl(26) ,F2(26) ,Height(26) 

2050  FOR  1=1  TO  15 


2060  Cl=l 

2070  FOR  J=150  TO  500  STEP  14 

2080  GOSUB  I  OF  1000,1030,1060,1090,1120,1150,1180,1210,1240,1270 

2090  IF  Kll  THEN  2110 

2100  GOSUB  I-IO  OF  1300,1330,1360,1390,1420 

2110  GOSUB  I  OF  2190,2220,2240,2260,2280,2300,2320,2340,2360,2380 

2120  IF  Kll  THEN  2140 

2130  GOSUB  I-IO  OF  2400,2420,2440,2460,2480 

2140  C1=C1+1 

2150  NEXT  J 


2160  NEXT  I 
2170  T9=0 
2180  GO  TO  3000 
2190  A(C1)=X 
2200  HeLght(Cl)=J 
2210  RETURN 
2220  B(C1)=X 
2230  RETURN 
2240  C(C1)=X 
2250  RETURN 
2260  D(C1)=X 
2270  RETURN 
2280  E1(C1)=X 
2290  RETURN 
2300  E2(C1)=X 
2310  RETURN 
2320  E3(C1)=X 
2330  RETURN 
2340  E4(C1)=X 
2350  RETURN 
2360  E5(C1)=X 
2370  RETURN 
2380  E6(C1)=X 
2390  RETURN 
2400  E7(C1)=X 
2410  RETURN 
2420  E8(C1)=X 
24  30  RETURN 
2440  E9(C1)=X 
2450  RETURN 
2460  F1(C1)=X 
24  70  RETURN 
2480  F2(C1)=X 
24  90  RETURN 

3000  PRINT  "ENTER  MAXIMUM  LINE  LENGTH 
3010  INPUT  G 
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3020  IF  G>0  THEN  30A0 

3030  STOP 

30A0  C5=0 

3050  IF  T9=0  THEN  3130     IREPLOT  WITH  SAME  PARAMETERS 

3060  IF  G<=3000  THEN  3100 

3070  WINDOW  G-1500,G, 0,500 

3080  VIEWPORT  10,140,6,92.5 

3090  GO  TO  3160 

3100  WINDOW  G/2,G, 0,500 

3110  G2=G*(130/1500) /2 

3120  GO  TO  3150 

3130  WINDOW  0,G, 0,500 

3140  G2=G'''(130/1500) 

3150  VIEWPORT  10 ,G2+10 , 6 , 92 . 5 

3160  PRINT  "DEVICE  No.   FOR  PLOTTED  OUTPUT?  "; 

3170  INPUT  P 

3180  IF  T9>0  THEN  3270     IREPLOT  WITH  SAME  PARAMETERS 

3190  PRINT  "ENTER  1  FOR  LABELS  ON  BS  LINES  OR  0  FOR  NO  LABELS  "; 

3200  INPUT  Q7 

3210  PRINT  "COVERAGE  LEVEL? 

3220  INPUT  C2 

3230  IF  Pol  THEN  3270     !  IF  NOT  PLOTTER 

3240  PRINT  "SPEED? 

3250  INPUT  S 

3260  PRINT  @1,32:"BY"!S 

3270  PAGE 

3280  PRINT  @1, 17:1. 125, 2. 576 

3290  MOVE  (1P:0,0 

3300  DRAW  (aP:G,0 

3310  DRAW  @P:G,500 

3320  DRAW  gP: 0,500 

3330  DRAW  dP : 0 , 0 

3340  MOVE  @P: 0,150 

3350  DRAW  (aP:G,150 

3360  FOR  1=0  TO  G  STEP  100 

3370         MOVE  @P:I,0 

3380         RMOVE  @P : 0 , 5 

3390         RDRAW  @P:0,-5 

3400  NEXT  I 

3410  FOR  1=0  TO  G  STEP  200 

3420         MOVE  dP: 1+21,0 

3430         PRINT  @P:"J_H_H_H_":I 

3440  NEXT  I 

3450  MOVE  (aP:G/2,0 

3460  IF  T9>0  THEN  3480 

3470  PRINT  @P: "J_J_J_H_H_H_H_H_H_H_H_H_H_H_H_H_DESIRED  LINE  LENGTH  (FEET)" 

3480  FOR  1=100  TO  500  STEP  100 

3490         MOVE  @P:0,I 

3500         RMOVE  @P:10,0 

3510         RDRAW  @P:-10,0 

3520         RMOVE  @P:0,-5 

3530        PRINT  @P:"H_H_H_":I 

3540  NEXT  I 

3550  IF  P=32  THEN  3630 

3560  MOVE  @P: 0,150 

3570  PRINT  @P : "H_H_H_H_H_" ; 

3580  PRINT  @1 ,25:90 
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3590  IF  T9>0  THEN  3610 

3600  PRINT  @P:"DROP  HEIGHT  (FEET  ABOVE  GROUND)" 
3610  PRINT  @1,25:0 
3620  IF  T9>0  THEN  3880 
3630  MOVE  @P:G/2,70 
3640  SET  DEGREES 

3650  PRINT  @P: "H_H_H_H_H_H_H_H_H_H_AIRCRAFT  SAFETY  LIMIT" 

3660  DIM  Safex(lO) ,Safey(10)      ! DRAW  SAFETY  LIMIT 

3670  Safex(l)=0 

3680  Safex(2)=G 

3690  Safex(3)=G 

3700  Safex(6)=G/2-lA0 

3710  Safex(A)=0 

3720  Safex(5)=0 

3730  Safex(7)=G/2+lA5 

37A0  Safex(8)=G/2+lA5 

3750  Safex(9)=G/2-lA0 

3760  Safex(10)=G/2-lA0 

3770  Safey(l)=0 

3780  Safey(7)=75+10 

3790  Safey(2)=0 

3800  Safey(3)=150 

3810  Safey(A)=150 

3820  Safey(5)=0 

3830  Safey(6)=75-10 

38A0  Safey(7)=75-10 

3850  Safey(8)=75+10 

3860  Safey(9)=75+10 

3870  Safey(10)=75-10 

3880  HATCH  ROTATE  -A5 

3890  HATCH  SPACE  A 

3900  HATCH  ALIGN  10,31.95 

3910  HATCH  @P : Saf ex, Saf ey 

3920  MOVE  @P: 0,500 

3930   IF  T9>0  THEN  3950 

39A0  PRINT  @P:"K_  COVERAGE  LEVEL:  ";C2 

3950  IF  V6=l  THEN  3980      !V6=1   IF  HIT  PROBABILITY  FILE  HAS  BEEN  READ 
3960  GOSUB  7000 
3970  CLOSE 

3980  PRINT  "ENTER  1  FOR  LINE  HP  OR  0  FOR  FULL  HP"; 

3990  INPUT  Q9 

AOOO  IF  Q9=l  THEN  A050 

AOlO  HATCH  ALIGN  10,31.95 

A020  HATCH  ROTATE  0 

A030  HATCH  SPACE  0.7 

AOAO  HATCH  @P : Hx , Hy 

A050  MOVE  @P:Hx(l) ,Hy(l) 

A060  DRAW  @P:Hx,Hy 

5000  MOVE  @P:A(1) ,HeLght(l)      ! DRAW  STRATEGY  CURVES 

5010  DRAW  @P: A, Height 

5020  IF  Q7=0  THEN  5050 

5030  MOVE  @P:A(15) ,Height(15) 

50A0  PRINT  (aP:"H_";AS 

5050  REM  GO  TO  19100 

5060  MOVE  @P:B(1) ,Hei6ht(l) 

5070  DRAW  @P:B, Height 

5080  IF  Q7=0  THEN  5110 
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5090  MOVE  @P:B(1A)  , Heighten) 

5100  PRINT  @P:"H_":B$ 

5110  REM  GO  TO  19100 

5120  MOVE  @P:C(1)  .Heightd) 

5130  DRAW  @P:C, Height 

5140  IF  Q7=0  THEN  5170 

5150  MOVE  @P:C( 13) .Height (13) 

5160  PRINT  @P: "H_" ;C$ 

5170  REM  GO  TO  19100 

5180  MOVE  @P:D(1) ,Height(l) 

5190  DRAW  @P:D. Height 

5200  IF  Q7=0  THEN  5230 

5210  MOVE  @P:D(12) ,Height(12) 

5220  PRINT  @P:"H_":D$ 

5230  REM  GO  TO  19100 

5240  MOVE  @P : El (1) .Height (1) 

5250  DRAW  gP : El , He ight 

5260  IF  Q7=0  THEN  5290 

5270  MOVE  (aP:El(ll)  ,Height(ll) 

5280  PRINT  @P:"H_":ES 

5290  REM  GO  TO  19100 

5300  MOVE  @P:E2(1) ,Height(l) 

5310  DRAW  @P:E2, Height 

5320  IF  Q7=0  THEN  5350 

5330  MOVE  @P:E2(10) ,Height(10) 

5340  PRINT  @P:"H_";F$ 

5350  REM  GO  TO  19100 

5360  MOVE  @P:E3(1) ,Height(l) 

5370  DRAW  dP : E3 , Height 

5380  IF  Q7=0  THEN  5410 

5390  MOVE  @P:E3(9) ,Height(9) 

5400  PRINT  @P: "H_" ;G$ 

5410  REM  GO  TO  19100 

5420  MOVE  @P:E4(1) ,Height(l) 

54  30  DRAW  @P : E4 , He ight 

5440  IF  Q7=0  THEN  5470 

5450  MOVE  (aP:E4(8)  ,Height(8) 

5460  PRINT  @P:"H_";H$ 

5470  REM  GO  TO  19100 

5480  MOVE  @P:E5(1) ,Height(l) 

5490  DRAW  @P:E5,Height 

5500  IF  Q7=0  THEN  5530 

5510  MOVE  @P:E5(7) ,Height(7) 

5520  PRINT  @P:"H_";I$ 

5530  REM  GO  TO  19100 

5540  MOVE  @P:E6(1) ,Height(l) 

5550  DRAW  (aP:E6,Height 

5560  IF  Q7=0  THEN  5590 

5570  MOVE  (aP:E6(6)  ,Height(6) 

5580  PRINT  @P:"H_":J$ 

5590  REM  GO  TO  19100 

5600  MOVE  @P:E7(1) ,Height(l) 

5610  DRAW  @P:E7,Height 

5620  IF  Q7=0  THEN  5650 

5630  MOVE  @P:E7(5) ,Height(5) 

5640  PRINT  @P:"H_":K$ 

5650  MOVE  @P:E8(1) ,Height(l) 
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5660  DRAW  (3P  :  E8  ,  He  ight 

5670  IF  Q7=0  THEN  5700 

5680  MOVE  @P:E8(4) .Height (4) 

5690  PRINT  (aP:  "H_":L$ 

5700  MOVE  @P:E9(1) .Height (1) 

5710  DRAW  @P:E9. Height 

5720  IF  Q7=0  THEN  5750 

5730  MOVE  @P:E9(8) ,Height(8) 

5740  PRINT  @P:"H_":M$ 

5750  MOVE  @P:F1(1) .Height(l) 

5760  DRAW  @P:Fl.Height 

5770  IF  Q7=0  THEN  5800 

5780  MOVE  gP:Fl(17) .Height(17) 

5790  PRINT  @P:"H_";N$ 

5800  MOVE  @P:F2(1) .Height(l) 

5810  DRAW  @P:F2. Height 

5820  IF  Q7=0  THEN  6000 

5830  MOVE  @P:F2(16) ,Height(16) 

5840  PRINT  @P:"H_":0$ 

6000  DIM  Capx(62) ,Capy(62)      !DRAW  NO  CAPABILITY  SECTION 
6010  Count=l 

6020  FOR  J=150  TO  700  STEP  10 

6030         GOSUB  1090 

6040  Capx(Count)=X 

6050         Capy (Count)=J 

6060  Count=Count+l 

6070  NEXT  J 

6080  Capx(55)=G 

6090  Capy(55)=500 

6100  Capx(56)=G 

6110  Capy(56)=150 

6120  Capx(57)=Capx(l) 

6130  Capy(57)=150 

6140  Startx=(G+D(20) ) /2 

6150  Starty=Height(20) 

6160  Capx(58)=Startx 

6170  Capy(58)=Starty 

6180  Capx(59)=Startx+160 

6190  Capy (59)=Starty 

6200  Capx(60)=Startx+160 

6210  Capy(60)=Starty+35 

6220  Capx(61)=Startx 

6230  Capy(61)=Starty+35 

6240  Capx(62)=Startx 

6250  Capy(62)=Starty 

6260  HATCH  SPACE  4 

6270  HATCH  ROTATE  45 

6280  HATCH  ALIGN  130  /  1500''G+10  .  31  .  95 
6290  HATCH  @P:Capx,Capy 

6300  MOVE  @P: (G+D(20) )/2+80.Height(20)+21 

6310  PRINT  dP: "H_NOJ_H_H_H_H_H_H_CAPABILITY" 

6320  T9=l 

6330  GO  TO  3000 

6340  END 

7000  DATA  0.12.8,25.1.35.3,43.7.50.6.56.6.61.1.66.69.9.73.6.77.2.80.5 
7005  Z=103 

7010  DATA  83. 7, 86. 6, 89. 3, 91. 8. 94. 1. 96. 2. 98. 2, 100. 1,101. 8. 103. 4. 104. 9 
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7  011  DATA  106. 3, 107. 7, 109, 110. 3, 111-5, 112. 6, 113. 7,11^.7,115. 8, 116. 9, 118 

7012  DATA  119.1,120.2,121.3,122.5,123.7,125,126.3,127.8,129.2,130.7 

7013  DATA  132. 2, 133. 8, 135. 4, 137. 1,138. 8,1 AO. 5, 142. 1,1 A3. 7,145.4,147.1 

7014  DATA  14  9.1,151.1,153.1,155.1,156.8,158.3,159.9,161.5,163.2,165.1 

7015  DATA  167.1,169,171,173.1,175.3,177.6,180.1,182.7,185.4,188.1,190.9 

7016  DATA  193.9,196.9,200,203.1,206.4,209.8,213.3,217,220.9,224.9,229.1 

7017  DATA  233.3,237.5,241.7,245.9,250.1,254.6,259.3,264.1,269,274.2 

7018  DATA  279.4,284.8,290.1,295.5,0,0 

7020  IF  Pol  THEN  7030 

7021  GIN  @1 iPausel ,Pause2 
7030  DIM  Hx(Z) ,Hy(Z) 
7035  RESTORE  7000 

7040  FOR  1=1  TO  Z 
7050         READ  Hx(I) 
7055  Hy(I)=(I-l)"5 
7060  NEXT  I 

7080  Hy(102)=500 

7081  Hy(103)=0 
7090  V6=l 
7100  RETURN 

8000  PRINT  "??     ":      !PRINT  A  KEYBOARD  ENTERED  LABEL  ON  PLOTTER 
8010  INPUT  L$ 

8020  PRINT  (31,17:1.125,2.576 
8030  PRINT  @1 :L$; 
8040  GO  TO  100 

9000  PRINT  "#  OF  TANKS-THEN  #  OF  REPEATS-THEN  #  OF  SECONDS?  "; 

9010  INPUT  L,M,N$      IPRINT  TXR  LABEL  ON  PLOTTER 

9020  GIN  (aP:Pl,P2 

9030  PRINT  (11,32:  "BY";  50 

9040  IF  N$="0"  THEN  9090 

9050  PRINT  @1, 17:1. 125, 2. 576 

9060  PRINT  @1  :  "H_H_"  ;L;  "Tx"  ;M:  "RJ_H_H_H_H_"  :N$;  "s"  ; 
9070  MOVE  @P:P1,P2 
9080  END 

9090  PRINT  (51  :  "H_H_"  ;L;  "Tx"  ;M;  "R"  ; 
9100  MOVE  @P:P1 ,P2 
9110  END 

10000  GIN  @1:X8,X9      ! DRAW  RIGHT  ARROW 

10010  PRINT  @1, 17:1. 125, 2. 576 

10020  FOR  1=0  TO  6 

10030         DRAW  @1 :X8-15+I,X9+6-I 

10040         DRAW  @1 :X8-15+I ,X9-6+I 

10050         DRAW  @1:X8,X9 

10060  NEXT  I 

10070  END 

11000  GIN  (ai:X8,X9      !  DRAW  LEFT  ARROW 
11010  PRINT  01,17:1.125,2.576 
11020  FOR  1=0  TO  6 
11030         DRAW  @1 :X8+15-I ,X9+6-I 
11040  DRAW  @1 :X8+15-I ,X9-6+I 

11050  DRAW  @1:X8,X9 

11060  NEXT  I 
11070  END 
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SMOOTHFOOT 

This  program  uses  data  from  the  program  LDRWIDTH  and  aligns  the  footprint  data  for  input  to  FINALFOOT.  Data  for 
each  gpc  is  entered  and  presented  graphically  on  the  computer  for  modification  to  eliminate  crossing  gpc  lines,  etc. 

1  GOTO  60000 

2  GOTO  63000 

3  END 

4  GOTO  100 
8  GOTO  1000 

12  GOTO  50000 

16  GOTO  51000 

20  GOTO  30000 

24  GOTO  10000 

28  GOTO  20000 

32  GOTO  40000 

36  GOTO  55000 

40  GOTO  1080 

76  GOTO  61000 

80  GOTO  62000 
100  REM  data  entry  section 
110  DELETE  L1,L2,L3,L4,L5,F1,F2,F3,F4,F5 

120  DIM  Ll[40,2],L2t40,2],L3[40,2],L4[40,2],L5[40,2],Fl[40,2],F2[40,2] 

130  DIM  F4[40,2] ,F5[40,2] ,F3[40,2] 

140  Fl  =  0 

150  F2  =  0 

160  F3  =  0 

170  F4  =  0 

180  F5  =  0 

190  LI  =  0 

200  L2  =  0 

210  L3  =  0 

220  L4  =  0 

230  L5  =  0 

240  CLEAR 

250  PRINT   "TITLE?  "; 
260  INPUT  T$ 

270  T$  =  "CONAIR  DC-6B  TANKER  454   "   S,  T$ 
280  FOR  I  =  1  TO  40 

290       PRINT  "ENTER  LEADING  EDGE  FEET  THEN  ALTITUDE  FOR  1  GPC  "; 

300       INPUT  L1[I,1] ,L1[I,2] 

310       IF  L1[I,1]=0  THEN  330 

320  NEXT  I 

330  FOR  I  =  1  TO  40 

340       PRINT  "ENTER  LEADING  EDGE  FEET  THEN  ALTITUDE  FOR  2  GPC  "; 

350       INPUT  L2[I,1] ,L2[I,2] 

360       IF  L2[I,1]=0  THEN  380 

370  NEXT  I 

380  FOR  I  =  1  TO  4  0 

390       PRINT  "ENTER  LEADING  EDGE  FEET  THEN  ALTITUDE  FOR  3  GPC  "; 

400       INPUT  L3[I,1] ,L3[I,2] 

410       IF  L3 [I , 1]=0  THEN  430 

420  NEXT  I 

430  FOR  I  =  1  TO  40 

440       PRINT  "ENTER  LEADING  EDGE  FEET  THEN  ALTITUDE  FOR  4  GPC  "; 
450       INPUT  L4[I,1] ,L4[I,2] 
460       IF  L4[I,1]=0  THEN  480 
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I^^Q  NEXT  I 

hSO  FOR  I  =  1  TO  40 

490       PRINT  "ENTER  LEADING  EDGE  FEET  THEN  ALTITUDE  FOR  5  GPC  " 

500       INPUT  L5[I,1] ,L5[I,2] 

510       IF  L5[I,1]=0  THEN  530 

520  NEXT  I 

530  FOR  I  =  1  TO  40 

540       PRINT  "ENTER  TRAILING  EDGE  FEET  THEN  ALTITUDE  FOR  1  GPC 

550       INPUT  F1[I,1] ,F1[I,2] 

560       IF  F1[I,1]=0  THEN  580 

570  NEXT  I 

580  FOR  I  =  1  TO  40 

590       PRINT  "ENTER  TRAILING  EDGE  FEET  THEN  ALTITUDE  FOR  2  GPC 

600       INPUT  F2[I,1] ,F2[I,2] 

610       IF  F2[I,1]=0  THEN  630 

620  NEXT  I 

630  FOR  I  =  1  TO  4  0 

640       PRINT  "ENTER  TRAILING  EDGE  FEET  THEN  ALTITUDE  FOR  3  GPC 

650       INPUT  F3[I,1]  ,F3[I,,2] 

660       IF  F3[I,11=0  THEN  680 

670  NEXT  I 

680  FOR  I  =  1  TO  4  0 

690      PRINT  "ENTER  TRAILING  EDGE  FEET  THEN  ALTITUDE  FOR  4  GPC 

700       INPUT  F4[I,1] ,F4[I,2] 

710       IF  F4[I,1]=0  THEN  730 

720  NEXT  I 

730  FOR  I  =  1  TO  40 

740      PRINT  "ENTER  TRAILING  EDGE  FEET  THEN  ALTITUDE  FOR  5  GPC 
750       INPUT  F5[I,1] ,F5[I,2] 
760       IF  F5[I,1]=0  THEN  780 
770  NEXT  I 

780  PRINT  "DATA  ENTERED" 

790  GOTO  60000 
1000  PRINT 
1030  PRINT  T$ 

1040  PRINT  "MAXIMUM  LENGTH?  "j 

1050  INPUT  M 

1055  SET  CLIP  OFF 

1060  SET  VIEWPORT  10,120,10,80 

1070  SET  WINDOW  0,M, 0,500 

1080  SET  LINE  COLOR  15 

1090  SET  TEXT  COLOR  15   |   SET  TEXT  STYLE  0   |   SET  TEXT  SIZE  12 

1091  CLEAR 
1095  MOVE  0,520 
1100  TEXT  T$ 

1110  PLOT  AXIS  50,100 

1120  FOR  I  =  500  TO  0  STEP  -100 

1130       MOVE  0,1+5 

1140       TEXT  ABS(I-500) 

1150       MOVE  0,1 

1160       DRAW  M,I 

1170  NEXT  I 

1180  FOR  I  =  0  TO  M  STEP  50 

1190      MOVE  1-10,0+5 

1200       TEXT  I 

1210  NEXT  I 

1230  SET  LINE  COLOR  2 
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1240  MOVE  0,500 

1250  FOR  I  =  1  TO  40 

1260       IF  L1(I,1]=0  THEN  1290 

1270       DRAW  L1[I,1] ,500-Ll[I,2] 

1280  NEXT  I 

1290  MOVE  0,500 

1300  FOR  I  =  1  TO  40 

1310       IF  F1[I,1]=0  THEN  1350 

1320       DRAW  F1[I,1] ,500-Fl[I,2] 

1330  NEXT  I 

1350  SET  LINE  COLOR  3 

1360  MOVE  0,500 

1370  FOR  I  =  1  TO  40 

1380       IF  L2[I,1]=0  THEN  1410 

1390       DRAW  L2[I,1] ,500-L2[I,2] 

1400  NEXT  I 

1410  MOVE  0,500 

1420  FOR  I  =  1  TO  40 

1430       IF  F2[I,1]=0  THEN  1460 

1440       DRAW  F2[I,1] ,500-F2[I,2] 

1450  NEXT  I 

1460  SET  LINE  COLOR  6 

1480  MOVE  0,500 

1490  FOR  I  =  1  TO  40 

1500       IF  L3[I,1]=0  THEN  1530 

1510       DRAW  L3[I,1] ,500-L3[I.2] 

1520  NEXT  I 

1530  MOVE  0,500 

154  0  FOR   I   =   1   TO   4  0 

1550       IF  F3[I,1]=0  THEN  1580 

1560       DRAW  F3[I,1] ,500-F3[I,2] 

1570  NEXT  I 

1580  SET  LINE  COLOR  14 

1600  MOVE  0,500 

1610  FOR  I  =  1  TO  40 

1620       IF  L4[I,1]=0  THEN  1650 

1630       DRAW  L4[I,1] ,500-L4[I,2] 

164  0  NEXT  I 

1650  MOVE  0,500 

1660  FOR  I  =  1  TO  40 

1670       IF  F4[I,1)=0  THEN  1700 

1680       DRAW  F4[I,1] ,500-F4[I,2] 

1690  NEXT  I 

1700  SET  LINE  COLOR  15 

1720  MOVE  0,500 

1730  FOR  I  =  1  TO  40 

1740       IF  L5[I,1]=0  THEN  1770 

1750       DRAW  L5[I,11 ,500-L5(I,2] 

17  60  NEXT  I 

1770  MOVE  0,500 

1780  FOR  I  =  1  TO  40 

1790       IF  F5[I,1]=0  THEN  1820 

1800       DRAW  F5[I,1] ,500-F5[I,2] 

1810  NEXT  I 

1820  END 

10000  REM  CHANGE  LEADING  EDGE 

10010  DELETE  L 


10020  DIM  L[A0,2] 

10030  HOME 

10040  PRINT 

10050  PRINT  "GPC?  "; 

10060  INPUT  G 

10070  ON  G  GOSUB  10180 , 10200 , 10220 , 102A0 , 10260 

10080  FOR  I  =  1  TO  40 

10090       PRINT  USING  " 4D , X , S" : L [ I , 1 ] 

10100  NEXT  I 

10110  PRINT 

10120  PRINT  "WHICH  ENTRY?  "; 
10130  INPUT  E 

lOlAO  PRINT  "NEW  VALUE?  "; 
10150  INPUT  L[E,1] 

10160  ON  G  GOTO  10280,10300,10320,10340,10360 

10170  END 

10180  L  =  LI 

10190  RETURN 

10200  L  =  L2 

10210  RETURN 

10220  L  =  L3 

10230  RETURN 

10240  L  =  L4 

10250  RETURN 

10260  L  =  L5 

10270  RETURN 

10280  LI  =  L 

10290  END 

10300  L2  =  L 

10310  END 

10320  L3  =  L 

10330  END 

10340  L4  =  L 

10350  END 

10360  L5  =  L 

10370  END 

20000  REM  CHANGE  TRAILING  EDGE 

20010  DELETE  F 

20020  DIM  F[40,2] 

20030  HOME 

20040  PRINT 

20050  PRINT  "GPC?  "i 

20060  INPUT  G 

20070  ON  G  GOSUB  20180,20200,20220,20240,20260 

20080  FOR  I  =  1  TO  40 

20090       PRINT  USING  " 4D , X , S " : F [ I , 1 ] 

20100  NEXT  I 

20110  PRINT 

20120  PRINT  "WHICH  ENTRY?  "j 
20130  INPUT  E 
20140  PRINT  "NEW  VALUE?  "; 
20150  INPUT  F[E,1] 

20160  ON  G  GOTO  20280,20300,20320,20340,20360 
20170  END 
20180  F  =  Fl 
20190  RETURN 
20200  F  =  F2 
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20210 
20220 
20230 
20240 
20250 
20260 
20270 
20280 
20290 
20300 
20310 
20320 
20330 
20340 
20350 
20360 
20370 
30000 
30005 
30010 
30020 
30030 
30040 
30050 
30060 
30070 
30075 
30076 
30080 
30081 
30090 
30091 
30100 
30101 
30110 
30111 
30120 
30130 
30140 
30150 
30160 
30170 
30180 
30190 
30200 
30210 
30220 
30230 
30240 
30250 
30260 
30270 
30280 
30290 
30300 
30310 
30320 


RETURN 
F  =  F3 
RETURN 
F  =  F4 
RETURN 
F  =  F5 
RETURN 
Fl  =  F 
END 

F2  =  F 
END 

F3  =  F 
END 

F4  =  F 
END 

F5  =  F 
END 

REM  PRINT  OUT  SECTION 
OPEN  #1 : "printf ile" , "f " 
DELETE  K1,K2,L,F 

DIM  Kl[5] ,K2[5] ,L[40,2] ,F[40,2] 
Kl  =  0 
K2  =  0 

P$  =  "X,4D,2X,4D,2X,S" 
CLEAR 

PRINT  T$    I    PRINT  #1:T$ 

PRINT    I    PRINT  #1: 

PRINT    I    PRINT  #1: 

PRINT  "  1  2 

PRINT  #1:"  1  2 

PRINT  "  5" 

PRINT  #1 : "  5" 

PRINT  " 


PRINT  #1  : 
PRINT  " 


PRINT  #1:"  " 

FOR  I  =  1  TO  40 

IF  L5[I,1]=0  AND  L4[I,1]=0  AND  L3[I,1]=0  AND  L2[I,1]=0  THEN  30150 
GOTO  30160 


IF  L1[I,1]=0  THEN  EXIT  TO  30230 
PRINT  USING  P$:H[I.l]  ,L1[I,2] 
PRINT  USING  P$:L2[I,1] ,L2[I,2) 
PRINT  USING  P$:L3[I,1] ,L3[I,2] 
PRINT  USING  PS:L4[I,1] ,L4[I,2] 
PRINT  USING  PS:L5(I,1] ,L5[I,2] 
PRINT  I  PRINT  #1 : 
NEXT  I 

FOR  I  =  1  TO  5 

ON  I  GOSUB  10180,10200,10220,10240,10260 
ON  I  GOSUB  20180,20200,20220,20240,20260 
FOR  J  =  1  TO  40 

IF  L[J,1]=0  THEN  EXIT  TO  30290 
NEXT  J 

IF  J=l  THEN  30310 
K1[I]   =  L[J-1,1] 
FOR  J  =  1  TO  40 

IF  F[J,1]=0  THEN  EXIT  TO  30340 


PRINT  #1  USING  P$:L1[I,1] ,L1(I,2] 
PRINT  #1  USING  PS:L2(I,1] ,L2[I,2] 
PRINT  #1  USING  P$:L3[I,1] ,L3[I,2] 
PRINT  #1  USING  P$:L4[I,1] ,L4[I,2] 
PRINT  #1  USING  P$:L5[I,1] ,L5[I,2] 
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30330       NEXT  J 

30340       IF  J=l  THEN  30360 

30350       K2[I]   =  F[J-1,1] 

30360  NEXT  I 

30370  Kl  =  K2+K1 

30380  Kl  =  Kl/2 

30390  FOR  I  =  1  TO  5 

30400       PRINT  USING  P$:K1[I],0    |   PRINT  #1  USING  PS:K1[I],0 

30410  NEXT  I 

30420  PRINT    |   PRINT  #1: 

30430  FOR  I  =  40  TO  1  STEP  -1 

30440       IF  F5[I,1]=0  AND  F4[I,1]=0  AND  F3[I,1]=0  AND  F2[I,1]=0  THEN  30460 

30450       GOTO  30470 

30460       IF  F1[I,1]=0  THEN  30530 


PRINT  #1  USING  P$:F1[I,1] ,F1[I,2] 
PRINT  #1  USING  P$:F2[I,1] ,F2[I,2] 
PRINT  #1  USING  P$:F3[I,1] ,F3[I,2] 
PRINT  #1  USING  P$:F4[I,1] ,F4[I,2] 
PRINT  #1  USING  P$:F5[I,1] ,F5[I,2] 


30470  PRINT  USING  P$ : Fl [ I , 1 ] . Fl ( I , 2] 
30480  PRINT  USING  PS : F2 [ I , 1 ] , F2 [ I , 2 ] 
30490  PRINT  USING  P$ : F3 [ I , 1 ] , F3 [ I , 2 ] 
30500  PRINT  USING  P$ : F4 [ I , 1 ] , [ I , 2 ] 
30510  PRINT  USING  P$ : F5 [ I , 1 ] , F5 [ I , 2] 
30520  PRINT  |  PRINT  #1: 
30530  NEXT  I 

30540  CALL  "Ipr  printfile" 
30550  END 

40000  PRINT  "WHICH  GPC?  "; 
40010  INPUT  G 

40020  ON  G  GOTO  1230,1350,1460,1580,1700 
40030  END 

50000  REM  SAVE  DATA  ON  DISC 
50010  PRINT  "FILE  NAME?  "j 
50020  INPUT  F$ 
50050  OPEN  #1:F$,"F" 

50060  WRITE  #1 : T$ , LI , L2 , L3 , L4 , L5 ,F1 , F2 , F3 , F4 , F5 
50070  CLOSE 
50080  END 

51000  REM  RECALL  STORED  DATA 

51010  DELETE  LI , L2 , L3 , L4 , L5 , Fl , F2 , F3 , F4 , F5 

51020  DIM  L1[40,2],L2[40,2],L3[40,2],L4[40,2],L5[40,2],F1[40,2],F2[40,2] 
51030  DIM  F4[40,2] ,F5[40,2] ,F3[40,2] 
51040  PRINT  "FILE  NAME  TO  BE  RECALLED?  "; 
51050  INPUT  F$ 
51060  OPEN  #1:F$,"R" 

51070  READ  #1 : T$ , LI , L2 , L3 , L4 , L5 , Fl , F2 , F3 , F4 , F5 
51080  CLOSE 
51090  END 

55000  PRINT  "WHICH  GPC?  "; 
55010  INPUT  G 

55020  PRINT  "WHAT  HEIGHT?  "; 
55030  INPUT  H 

55040  PRINT  "WHAT  AMOUNT?  "; 
55050  INPUT  A 

55060  ON  G  GOTO  55100,55200,55300,55400,55500 
55070  PRINT  "REENTER" 
55080  STOP 

55100  L1[H/100,1]  =  L1[H/100,1]+A 
55110  F1[H/100,1]  =  F1[H/100,1]+A 
55120  END 

55200  L2[H/100,1]  =  L2 [H/ 100 , 1 ] +A 


77 


55210 

F2[H/100,1]   =  F2[H/100,1]+A 

55220 

END 

55300 

L3[H/100,1]   =  L3[H/100,1]+A 

55310 

F3[H/100,1]   =  F3[H/100,1]+A 

55320 

END 

55400 

L4[H/100,1]   =  L4[H/100,1]+A 

55410 

F4[H/100,1]   =  F4[H/100,1]+A 

55420 

END 

55500 

L5[H/100,1]   =  L5[H/100,1]+A 

55510 

F5[H/100,1]  =  F5[H/100,1]+A 

55520 

END 

60000 

SET  DIALOG  COLOR  3 

60005 

PRINT  "USER  DEFINABLE  KEY  1: 

ENTER  DATA" 

60010 

PRINT  "  2: 

DRAW  DATA  SET" 

60020 

PRINT  "  3: 

STORE  DATA" 

60030 

PRINT  "  4: 

RECALL  DATA" 

60040 

PRINT  "  5; 

PRINT  RESULTS" 

60050 

PRINT  "  6: 

CHANGE  LEADING  EDGE" 

60060 

PRINT  "  7: 

CHANGE  TRAILING  EDGE" 

60070 

PRINT  "  8: 

DRAW  ANY  ONE  GPC" 

60075 

PRINT   "  9: 

CHANGE  BOTH  LEADING  AND  TRAILING" 

60080 

PRINT   "  10: 

REDRAW  DATA  SET" 

60090 

PRINT   "  19: 

RECALL  A  PATDAT  FILE" 

60100 

PRINT  "  20: 

SAVE  THE  RECALLED  PATDAT  FILE" 

60110 

END 

61000 

PRINT  "FILE?  ": 

61010 

DELETE  M1,M2 

61020 

DIM  Ml[40] ,M2[8] 

61030 

INPUT  F$ 

61040 

F$  =  "PATDAT"   k  F$ 

61050 

OPEN  #1 :F$, "R" 

61060 

READ  #1:M1,M2 

61070 

CLOSE 

61075 

CLEAR 

61079 

FOR  J  =  0  TO  30  STEP  10 

61080 

FOR  I  =  1  TO  10 

61081 

PRINT  USING  "3D.3D,S":M1[I+J] 

61082 

NEXT  I 

61083 

PRINT 

61084 

NEXT  J 

61090 

SET  WINDOW  0,40,0,15 

61100 

SET  VIEWPORT  10,120,10,50 

61110 

MOVE  0 , 0 

61120 

DRAW  4  0,0 

61130 

MOVE  0,0 

61140 

FOR  I  =  1   TO  4  0 

61150 

DRAW  I, Ml [I] 

61160 

NEXT  I 

61170 

MOVE  0,1 

61180 

DRAW  30,1 

61190 

SET  TEXT  COLOR  2   |   TEXT  1 

61200 

MOVE  0,2 

61210 

DRAW  30,2 

61220 

SET  TEXT  COLOR  3    |    TEXT  2 

61230 

MOVE  0,3 

61240 

DRAW  30,3 

61250 

SET  TEXT  COLOR  6   |   TEXT  3 
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61260  MOVE  0,4 
61270  DRAW  30,4 

61280  SET  TEXT  COLOR  14   |   TEXT  4 
61290  MOVE  0,5 
61300  DRAW  30,5 

61310  SET  TEXT  COLOR  15    |   TEXT  5 
61320  MOVE  12,15 
61330  TEXT  F$ 
61340  END 

62000  OPEN  #1 :F$, "F" 
62010  WRITE  #1:M1,M2 
62020  CLOSE 
62030  END 

FINAL  FOOT 

This  program  accepts  data  for  footprint  and  width  and  draws  the  footprint  camera  ready  on  a  plotter.  Inputs  are  width 
for  hne  lengths  of  1,  2,  3,  4,  and  5  gpc  at  100-foot  drop  height  intervals,  tension,  and  downrange  position  of  the  footprints. 

Tension  is  the  tightness  of  the  smoothing  of  the  data  through  the  points,  the  lower  the  value  entered,  the  looser  the  fit, 
for  most  guidelines  a  value  of  1.0  appears  to  be  best. 

The  length  portion  of  the  footprints  for  each  gpc  always  starts  at  0,0  then  the  value  at  100  feet,  200,  etc.,  to  500  feet; 
then  starts  up  the  right  side  of  the  footprint  at  500  feet,  400  feet,  etc.  The  next  to  last  point  is  0  feet  and  the  distance  the 
aircraft  flies  while  the  tanks  are  emptying  (evacuation  time  x  speed  in  feet  per  second),  with  the  last  point  being  0,0. 

Two  files,  SIDE  and  END,  containing  a  digitized  outline  of  the  side  and  end  view  of  the  aircraft  must  also  be  available 
to  the  program.  Each  is  organized  as  N,  X,  X^„„  X^):^,,  Y^„„  Y^.  where  N  is  the  number  of  X,  Y  data  points. 

1  GO  TO  2000 
4  GO  TO  100 
8  GO  TO  1000 
12  GO  TO  1090 
16  GO  TO  1900 

20  K9=0 

21  VIEWPORT  G8,G9,20,60 

22  WINDOW  0,90,0,500 

23  GO  TO  2500 

24  GO  TO  2300 

28  PRINT  "WHAT  GPC?   " ; 

29  INPUT  Cgpc 

30  GO  TO  Cgpc  OF  7500,7600,7700,7800,7900 
32  GO  TO  5500 

36  GO  TO  1700 
40  GO  TO  8200 
44  GO  TO  6200 
48  GO  TO  8000 
52  GO  TO  60000 
56  GO  TO  63000 
60  GO  TO  65000 
64  GO  TO  64000 

80  PRINT  "ARRAY?" ; 

81  INPUT  Ard 

82  GO  TO  Ard  OF  51000,52000,53000,54000,55000 

89  END 

90  DELETE  AS 

91  DIM  A$(300) 

92  CALL  "MOUNT", 0,A$ 

93  PRINT  AS 

94  END 
100  PAGE 
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110  PRINT  "INPUT  WIDTHS:     GIVE  WIDTH  FOR  EACH  CONCENTRATION  AT  EACH" 

120  PRINT  "ALTITUDE.      IF  NO  WIDTH  FOR  GIVEN  ALTITUDE  ENTER  0." 

130  DIM  Xl(6) ,X2(6) ,X3(6) ,X/i(6) ,Y1(6) ,Y2(6) ,Y3(6) ,Y4(6) ,Y9(6) 

140  DIM  Cl(6) ,C2(6) ,H1(6) ,H2(6) ,H3(6) ,HA(6) ,P(20) 

150  P9=19 

160  X1(1)=0 

170  X2(l)=0 

180  X3C1)=0 

190  X4(l)=0 

200  H1(1)=0 

210  IMAGE  "INPUT  WIDTHS  FOR  ",D,"  GPC" 
220  PRINT  USING  210:1 
230  FOR  1=2  TO  6 

240         PRINT   (I-l)'nOO;"  FEET  "; 
250         INPUT  XI (I) 
260  NEXT  I 

270  PRINT  USING  210:2 
280  FOR  1=2  TO  6 

290         PRINT   (I-l)'''100;  "  FEET  "; 
300         INPUT  X2(I) 
310  NEXT  I 

320  PRINT  USING  210:3 
330  FOR  1=2  TO  6 

3A0  PRINT    (I-l)'-'lOO;  "   FEET  "; 

350         INPUT  X3(I) 
360  NEXT  I 

370  PRINT  USING  210 : 4 
380  FOR  1=2  TO  6 

390         PRINT   d-D-nOO;"  FEET  "; 
400         INPUT  X4(I) 
AlO  NEXT  I 

420  PRINT  USING  210:5 
4  30  FOR  1=2  TO  6 

440         PRINT   d-D-nOO;"  FEET  "; 

450         INPUT  H1(I) 

4  60  NEXT  I 

4  70  Xl=Xl/2 

480  X2=X2/2 

490  X3=X3/2 

500  X4=X4/2 

510  Hl=Hl/2 

520  FOR  1=1  TO  6 

530         Y9(I)  =  (I-1)'''100 

54  0  NEXT  I 

550  Yl=500-Y9 

560  Y2=Y1 

570  Y3=Y1 

580  Y4=Y1 

590  H2=Y1 

600  PRINT  "INPUT  TENSION  (0<T<100)": 
610  INPUT  P(2) 
620  PAGE 

630  DELETE  X5 , Y5 , X6 . Y6 , X7 , Y7 , X8 , Y8 , X9 , Y9 , A2 
640  DIM  A2(40,2) 
650  FOR  J=l  TO  5 
660  K=0 

670         FOR  1=1  TO  40 
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680 


PRINT  "ENTER  DISTANCE  DOWN  RANGE  THEN  DROP  HEIGHT  FOR 


:  J: 


GP 


C 


690 


INPUT  A2(I,1) ,A2(I,2) 


700 


A2(I,2)=500-A2(I,2) 


710 


IF  A2(I,2)>501  THEN  7A0 


720 


K=K+1 


730        NEXT  I 

740         GOSUB  J  OF  6670,6800,6910,7020,7130 
750  PAGE 
760  NEXT  J 
770  END 
1000  PAGE 

1010  PRINT  "WHAT  IS  THE  LENGTH  OF  THE  GRAPH?  ": 
1020  INPUT  L 

1030  PRINT  "NUMBER  OF  COMPARTMENTS?  "; 
1040  INPUT  T$ 

1050  PRINT  "NUMBER  OF  GALLONS?  "; 
1060  INPUT  L2 

1070  REM  PRINT  "AFTER  CPT  LABEL?  "; 
1080  REM  INPUT  L$ 

1090  PRINT  "ENTER  1  FOR  PLOTTER,    32  FOR  SCREEN"; 
1100  INPUT  D 
1110  Q9=D 

1120  IF  Dol  THEN  1200 

1130  PRINT  "ENTER  PLOTTER  SPEED:    " ; 

1140  INPUT  By 

1150  PRINT  (21,32:  "BY", By 

1200  PAGE 

1210  REM  G8=L/6.299212598A3+25 

1211  G8=18+L'm2/1000 

1220  G9=G8+8. 74285714286" (1300/ 800) 

1221  G9=90/1000'ni2+G8 
1230  GOSUB  6370 

1250  WINDOW  0,1000,0,500 

1260  PRINT  (51,17:1.125,2.576 

1270  VIEWPORT  12,124,20,60 

1280  FOR  1=500  TO  0  STEP  -100 

1290         MOVE  (30:0,1 

1300         RDRAW  @D:L,0 

1310         RDRAW  @D:-L,0 

1340         RMOVE  (30:0,-5 

1350         PRINT  @0:"H_H_H_"; 500-1 

1360  NEXT  I 

1370  MOVE  @D:0,0 

1380  DRAW  @D: 0,500 

1390  DRAW  (50:0,0 

1400  Z7=l 

1410  FOR  1=50  TO  L  STEP  50 
1420         MOVE  @0:I,0 
1430         IF  Z7=3  THEN  1480 
1440         RDRAW  @D : 0 , 6 
1450         RDRAW  @D:0,-6 
1460  Z7=Z7+1 
1470         GO  TO  1530 
1480         RDRAW  @D:0,12 
1490         RDRAW  @D : 0 , -12 
1500         RDRAW  @D:-2.0 
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1510         PRINT  @D:"J_H_":I 

1520  Z7=l 

1530  NEXT  I 

15A0  MOVE  @D:L/2,0 

1550  PRINT  @D: "J_J_H_H_H_H_H_H_H_H_H_H_H_H_H_H_H_H_RANGE  AND  PATTERN  LENGTH  (FEET)' 

1560  IF  D=l  THEN  1580 

1570  GO  TO  1640 

1580  MOVE  (11:0,250 

1590  PRINT  "H_H_H_H_" 

1600  SET  DEGREES 

1610  PRINT  @1, 25:90 

1620  PRINT  (30 : "K_K_H_H_H_H_H_H_H_H_H_H_H_H_H_H_H_DROP  HEIGHT  (FEET  ABOVE  GROUND)" 
1630  PRINT  @1,25:0 
1640  MOVE  @D:20,500 

1650  PRINT  QD:"K_       ":T$;"   "; L2 ;" -GALLON  DROP" 

1659  PRINT  "G_G_G_G_" 

1660  END 

1700  VIEWPORT  G8,G9,20,60 

1710  WINDOW  0,90,0,500 

1720  AXIS  @D:30,500 

1725  PRINT  @1, 17:1. 125, 2. 576 

1730  MOVE  @D:0,0 

1740  PRINT  QD: "J_0" 

1750  FOR  1=30  TO  60  STEP  30 

1760         MOVE  @D:I,0 

1770         RDRAW  @D:0,10 

1780        RDRAW  @D: 0,-10 

1790         PRINT  @D:"J_H_";I 

1800  NEXT  I 

1810  MOVE  @D:0,0 

1820  PRINT  (50: "J_J_HALF-WIDTHJ_H_H_H_H_H_H_H_H_(FEET)" 

1830  FOR  1=100  TO  500  STEP  100 

1840  MOVE  @D:0,I 

1850         RORAW  @D:90,0 

1860         RDRAW  @D:-90,0 

1870  NEXT  I 

1880  PRINT  "G_" 

1890  END 

1900  WINDOW  0,1000,0,500 

1910  VIEWPORT  12,124,20,60 

1920  MOVE  @0: 0,500 

1930  K9=l 

1940  07=1 

1950  GO  TO  2500 

1960  END 

2000  PAGE 

2010  DIM  P(20) 

2020  C9=0 

2030  PRINT  "USER  DEFINABLE  KEYS:" 
2040  PRINT 


2050  PRINT 
2060  PRINT 
2070  PRINT 
2080  PRINT 
2090  PRINT 
2100  PRINT 
2110  PRINT 


1  ENTER  DATA   (WIDTHS,    TENSION,   FILE  #'S,  etc.)" 

2  ENTER  PARAMETERS  FOR  AXES  AND  DRAW" 

3  DRAW  AXES  FROM  PREVIOUSLY  ENTERED  PARAMETERS" 

4  DRAW  FOOTPRINTS" 

5  DRAW  WIDTHS" 

6  PRINT  LABELS  ON  PLOTTER" 

7  INSERT  VALUE  INTO  ANY  GPC  FOOTPRINT  " 
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2120  PRINT  "  8      PRINT  DATA  ON  PRINTER" 

2130  PRINT  "  9      DRAW  WIDTH  PORTION  OF  AXES" 

2140  PRINT  "10       PRINT  ID  ON  LOWER  LEFT  OF  GRAPH" 

2150  PRINT  "11       PRINT  LABELS  ON  WIDTH  CHART  WITH  LINES" 

2160  PRINT  "12      PRINT       STATEMENT  ON  GRAPH" 

2163  PRINT  "13       PRINT  X'S  ON  FOOTPRINT" 

2170  PRINT  "14       PRINT  X'S  ON  HALFWIDTH" 

2180  PRINT  "15       SAVE  ENTERED  DATA  ON  SISC  " 

2181  PRINT  "16       RECALL  ENTERED  DATA  FROM  DISC  " 

2189  PRINT  "20       CHANGE  FOOT  PRINT  ARRAY" 

2190  END 
2300  PAGE 

2310  PRINT  "PRESS  RETURN  WITH  NO  ENTRY  TO  EXIT  THIS  ROUTINE" 

2320  PRINT  "POSITION  PEN  THEN  INPUT  LABEL  AND  HIT  RETURN.  " 

2321  PRINT  @1, 17:1. 125, 2. 576 
2330  PRINT 

2340  PRINT  "INPUT  LABEL:"; 

2350  INPUT  A$ 

2360  IF  A$=""  THEN  2390 

2370  PRINT  @D:A$i "H_" ;A$i "H_" ; 

2380  GO  TO  2340 

2390  END 

2400  P(12)=l 

2500  DELETE  X,Y,C1,C2 

2540  07=1 

2550  IF  K9=l  THEN  8260 
2560  MOVE  (30:0,500 
2570  N=6 

2580  DIM  X(6) ,Y(6) ,C1(6) ,C2(6) 
2590  P(12)=l 

2600  GO  TO  07  OF  4330,4390,4430,4470,4510,4550 
2610  END 

2620  GO  TO  07  OF  8280,8370,8460,8560,8650,8740 

2630  DELETE  PS 

2635  IF  N=l  THEN  4550 

2636  Adot=SUM(X) 

2637  IF  Adot=0  THEN  4550 
2640  DIM  P8(N) 

2650  P3=X(2)-X(1) 

2660  P5=Y(2)-Y(1) 

2670  Q1=SQR(P3''P3+P5"P5) 

2680  Q3=P3/Q1 

2690  Q5=P5/Q1 

2700  IF  N>2  THEN  2750 

2710  C1=0 

2720  C2=0 

2730  GO  TO  3560 

27  4  0  REM     ''"VVrFIND  SLOPES  AT  ENDPOINTS''"'"' 

2750  P2=SQR( (X( 3 ) -X( 2 ) ) ^2+ ( Y ( 3 ) - Y ( 2 ) ) ^2 ) 

2760  Q4=Q1+P2 

2770  P3=-(Q4+Q1)  /  (Q4'''Q1) 

2780  P4=Q4/  (Q1'''P2) 

2790  P5=-Q1/ (Q4'''P2) 

2800  P1=P3'''X(1)+P4'''X(2)+P5'''X(3) 

2810  P2=P3'''Y(1)+P4'''YC2)+P5-'Y(3) 

2820  P3=SQR(P1'''P1+P2"P2) 

2830  C1(1)=Q3-P1/P3 
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2840  C2(1)=Q5-P2/P3 

2850  P2=SQR(  (X ( N-2 ) -X ( N- 1)  ) '2+ ( Y  (N-2 )  - Y  ( N- 1) ) 

2860  P6=SQR( (X(N) -X(N-l) ) ^2+( Y(N) -Y(N-1 ) ) ^2 ) 

2870  Q4=P2+P6 

2880  P3=(Q4+P6)  /  (Q4'''P6) 

2890  P4=-Q4/ (P6'-'P2) 

2900  P5=P6/ (Q4''P2) 

2910  Cl(N)=P3''X(N)+P4'-X(N-l)+P5''X(N-2) 
2920  C2(N)=P3"Y(N)+P4'''Y(N-l)+P5'''Y(N-2) 
2930  REM     '■"•"»SET  UP  TRI-DIAGONAL  SYSTEM*''"^ 
2940  P8(1)=Q1 
2950  Q6=Q1 

2960  IF  N<3  THEN  3120 
2970  P2=2 

2980  FOR  Pl=3  TO  N 

2990  P3=X(P1)-X(P2) 

3000  P5=Y(P1)-Y(P2) 

3010  Q2=SQR(P3--''P3+P5'-'P5) 

3020  P4=P3/Q2 

3030  P6=P5/Q2 

3040  C1(P2)=P4-Q3 

3050  C2(P2)=P6-Q5 

3060  P8(P2)=Q2 

3070  Q3=P4 

3080  Q5=P6 

3090  Q6=Q6+Q2 

3100  P2=P1 

3110  NEXT  PI 

3120  P(11)=Q6 

3130  P3=SQR(C1(N)*2+C2(N)^2) 
3140  C1(N)=C1(N) /P3-Q3 
3150  C2(N)=C2(N) /P3-Q5 

3160  REM     '"'"■'UN-NORMALIZE  TENSION  FACTOR**"' 

3170  P0=ABS(P(2)  )'-'CN-l) /Q6 

3180  REM     '"■"■'FORWARD  ELIMINATION'''"* 

3190  P3=P0"P8(1) 

3200  P4=EXP(P3) 

3210  P5=0 . 5*(P4-1/P4) 

3220  P6=l/ (P8(1)"P5) 

3230  Q1=P6"(P3"0 .5"(P4+1/P4)-P5) 

3240  Q5=1/Q1 

3250  C1(1)=Q5"C1(1) 

3260  C2(1)=Q5''C2(1) 

3270  Q4=P6''(P5-P3) 

3280  P8(1)=Q5''Q4 

3290  IF  N<3  THEN  3450 

3300  P2=l 

3310  FOR  Pl=2  TO  N-1 

3320  P3=P0"P8(P1) 

3330  P4=EXP(P3) 

3340         P5=0 .5''(P4-1/P4) 

3350         P6=l/ (P8(P1)''P5) 

3360  Q2=P6''(P3''0.5'MP4+1/P4)-P5) 

3370         Q5=l/ (Q1+Q2-Q4'*P8(P2) ) 

3380        C1(P1)=Q5''(C1(P1)-Q4"C1(P2) ) 

3390        C2(P1)=Q5'MC2(P1)-Q4''C2(P2) ) 

3400  Q4=P6''(P5-P3) 


3410  P8CP1)=Q5'''Q4 
3420  Q1=Q2 
3430  P2=P1 
3440  NEXT  PI 
3450  P2=N-1 

3460  Q5=l/ (Q1-Q4'-'P8(P2)  ) 

3470  C1(N)=Q5'-'(C1(N)-Q4'''C1(P2)  ) 

3480  C2(N)=Q5''(C2(N)-Q4'''C2(P2)  ) 

3490  REM     ^'"''^'"'f^'f'fBACK  SUBSTITUTION^'''''"'^'^*^'^'^'^^^"'''''^'^"''^'^'^'^' 
3500  ?2=N 

3510  FOR  P1=N-1  TO  1  STEP  -1 

3520  C1(P1)=C1(P1)-P8(P1)'''C1(P2) 

3530  C2(P1)=C2(P1)-P8(P1)"C2(P2) 

3540  P2=P1 

3550  NEXT  PI 

3560  P(13)=P9 

3570  DELETE  P8 

3580  GO  TO  3680 

3590  IF  P9019  THEN  3540 

3600  GOSUB  64 

3610  GO  TO  2600 

3620  END 

3630  GO  TO  8  OF  4560,4630,4700,4770,4840,39000 

3640  IF  P9>19  THEN  100 

3650  P9=P9+30 

3660  GOSUB  16 

3670  REM 

3680  REM 

3690  P5=0 

3700  GOSUB  3900 

3710  REM 

3720  DELETE  11,12 

3730  DIM  11(400) , 12(400) 

3740  11=0 

3750  12=0 

3760  FOR  Pl=l  TO  400 
3770  P5=P5+P0 
3780         GOSUB  3900 
3790  I1(P1)=P6 
3800  I2(P1)=P7 
3810  NEXT  PI 

3820  DRAW  @D:I1,I2 

3821  REM  '■'VtiV>v>v,v,vv.,wc.vvc/rVc*ABOVE  IS  DRAW  ROUTINE''''"^''"''^''"'"'"'"'"'' 
3840  MOVE  @D:0,0 

3850  DRAW  @D: 0,500 
3870  IF  K9=l  THEN  2620 
3880  GO  TO  2600 
3890  END 

3900  Q5=ABS(P(2)  )'''(N-P(12) )  /  P(ll) 
3910  Q7=ABS(P5'''P(11)  ) 
3920  IF  P5<0  THEN  4000 
3930  P3=2 
3940  P4=0 

3950  IF  P5>0  THEN  4000 
3960  IF  P(12)=0  THEN  4000 
3970  P6=X(1) 
3980  P7=Y(1) 
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3990  RETURN 

4000  IF  P3=l  THEN  4200 

AOlO  REM     '■"■"■'DETERMINE  WHICH  SEGMENT  P5  IS  MAPPED  INTO"** 

4020  Q8=P3-1 

4030  FOR  Q6=P3  TO  N 

4040  Q1=X(Q6)-X(Q8) 

4050  Q2=Y(Q6)-Y(Q8) 

4060  Q3=SQR(Q1"Q1+Q2"Q2) 

4070         IF  P4+Q3=>Q7  THEN  4180 

4080  P4=P4+Q3 

4090  Q8=Q6 

4100  NEXT  Q6 

4110  IF  P(12)=0  THEN  4150 

4120  P6=X(N) 

4130  P7=Y(N) 

4140  RETURN 

4150  Q6=l 

4160  Q8=N 

4170  Q3=P(11)-P4 

4180  P3=Q6 

4190  REM     '■"■"■'PERFORM  INTERPOLATION''"'* 

4200  Q1=Q7-P4 

4210  Q2=Q3-Q1 

4220  Q4=EXP(Q5'-'Q1) 

4230  Q6=0.5'''(Q4-1/Q4) 

4240  Q5=EXP(Q5*Q2) 

4250  Q7=0  .5'MQ5-1/Q5) 

4260  Q5=Q4'''Q5 

4270  Q5=0.5'''(Q5-1/Q5) 

4280  P6=(C1(P3)'^'Q6+C1(Q8)*Q7)  /Q5 

4290  P6=P6+( (X(P3)-C1(P3))'^Q1+(X(Q8)-C1(Q8) )*Q2) /Q3 
4300  P7  =  (C2(P3)'^'Q6+C2(Q8)>-'Q7)  /Q5 

4  310  P7=P7+(  (Y(P3)-C2(P3)  )''Q1  +  (Y(Q8)-C2(Q8)  )'''Q2)  /Q3 

4320  RETURN 

4330  DELETE  X 

4340  DIM  X(6) 

4  350  X=X1 

4360  Y=Y1 

4370  07=2 

4380  GO  TO  2630 

4  3  90  X=X2 

4400  Y=Y2 

4410  07=3 

4420  GO  TO  2630 

4  4  30  X=X3 

4440  Y=Y3 

4450  07=4 

4460  GO  TO  2630 

4470  X=X4 

4480  Y=Y4 

4490  07=5 

4500  GO  TO  2630 

4510  X=H1 

4520  Y=H2 

4530  07=6 

4540  GO  TO  2630 

4550  END 
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4560  DELETE  X,Y 

4570  DIM  X(N5) , Y(N5) 

4580  N=N5 

4590  X=X5 

4600  Y=Y5 

4610  08=2 

4620  GO  TO  11410 

4630  DELETE  X,Y 

4640  DIM  X(N6) , Y(N6) 

4650  N=N6 

4660  X=X6 

4670  Y=Y6 

4680  08=3 

4690  GO  TO  11410 

4700  DELETE  X,Y 

4710  DIM  X(N7) , Y(N7) 

4720  N=N7 

4730  X=X7 

4740  Y=Y7 

4750  08=4 

4760  GO  TO  11410 

4770  DELETE  X,Y 

4780  DIM  X(N8) , Y(N8) 

4790  N=N8 

4800  X=X8 

4810  Y=Y8 

4820  08=5 

4830  GO  TO  11410 

4840  END 

4850  N=N5 

4860  GO  TO  11220 

4870  N=N6 

4880  GO  TO  11220 

4890  N=N7 

4900  GO  TO  11220 

4910  N=N8 

4920  GO  TO  11220 

4930  END 

4940  GO  TO  1106 

4950  VIEWPORT  G8,G9,20,60 

4960  PRINT  13931 

4970  GO  TO  1350 

4980  VIEWPORT  12,124,20,60 

4990  PRINT  13951 

5000  GO  TO  19 

5010  VIEWPORT  G8,G9,20,60 

5020  PRINT  13971 

5030  GO  TO  67 

5040  VIEWPORT  12,124,20.60 

5050  PRINT  13991 

5060  GO  TO  8200 

5070  VIEWPORT  G8,G9,20,60 

5080  PRINT  14011 

5090  WINDOW  0,90,0,500 

5100  GO  TO  2500 

5110  REM  WHAT  IS  THIS  FOR?'"'  FIND  9 

5120  REM  INPUT  @33:AS 
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5130  REM  PRINT  A$ 

5140  REM  GO  TO  5120 

5150  REM  '"'"''Subroutine  to  determine  smallest  x  vlue   in  array'"'''' 

5160  M=500 

5170  FOR  1=1  TO  C 

5180        M=M  MIN  A2(I,2) 

5190  NEXT  I 

5200  RETURN 

5500  PAGE 

5501  PRINT   "RET  TYPE?  "; 
5510   INPUT  R$ 

5512  PRINT  "J_DEVICE?  ": 
5514   INPUT  Dl 

5530  PRINT  001 : "L_J_J_J_K_K_" : "     ";T$;"   "; L2 :" -GALLON  DROP  ":R$ 
5540  PRINT  @D1:"J_J_  WIDTHSJ_" 
5550  FOR  1=2  TO  6 

5560  PRINT  001;  USING  5980 : XI ( I) "2 , X2 ( I ) "2 , X3 ( I) "2 , X4 ( I ) "2 , HI ( I ) "2 
5570  PRINT  @D1:  USING  6010 : 500-Yl ( I) , 500-Y2 ( I) , 500-Y3 ( I) , 500-Y4 ( I) 
5580         PRINT  dDl:    USING  5990 : 500-H2 ( I) 

5590         PRINT  @D1:"  " 

5600  NEXT  I 
5610  DELETE  J$ 
5620  DIM  J$(132) 

5630  J$="         I  I  I  I  I  I  I" 

5640  J$=J$&"  I  I  I  I  I  I  " 

5650  J$=J$St"  I" 

5660  PRINT  @D1:"J_J_      DOWNRANGES  &  HEIGHTSJ_" 
5670  FOR  1=2  TO  N5 

5680         PRINT  dDl:   USING  6000 : X5 ( I , 1) , " | " ' 500-X5 ( I , 2 ) 
5690  NEXT  I 
5700  PRINT  dDl: 
5710  PRINT  (3D1:J$ 
5720  FOR  1=2  TO  N6 

5730         PRINT  ODl;    USING  6000 : X6 ( I , 1 ) , " | " ■ 500-X6 ( I , 2 ) 
5740  NEXT  I 
5750  PRINT  (501: 
5760  PRINT  @01:J$ 
5770  FOR  1=2  TO  N7 

5780        PRINT  @D1:   USING  6000 : X7 ( I , 1) , " | " , 500-X7 ( I , 2) 
57  90  NEXT  I 
5800  PRINT  @D1: 
5810  PRINT  @D1:J$ 
5820  FOR  1=2  TO  N8 

5830         PRINT  @D1:   USING  6000 : X8 ( I , 1) , " | " , 500-X8 ( I , 2 ) 
5840  NEXT  I 
5850  PRINT  (301  : 
5860  PRINT  @D1:J$ 
5870  FOR  1=2  TO  N9 

5880         PRINT  @D1:   USING  6000 : X9 ( I , 1) , " | " , 500-X9 ( I , 2 ) 

5890  NEXT  I 

5900  PRINT  @D1: 

5910  PRINT  (a01:J$ 

5920  FOR  1=1  TO  132 

5930        PRINT  (501  :  "_"  ; 

5940  NEXT  I 

5950  PRINT  @D1: 

5970  END 
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5980   IMAGE  5(5D) 
5990   IMAGE  5D 
6000  IMAGE  4D,A,3D,X,S 
6010  IMAGE  A(5D) ,S 
6020  END 

6030   INPUT  (133:  AS 

6200  REM  DRAW  HORIZ  LINES  AND  PRINT  GPC  ON  WIDTH 
6205  Ard=0 

6210  PRINT  (11,17:1.125,2.576 

6220  PRINT  "POSITION  PEN  WHERE  CHARACTER  WILL  BE  WRITTEN  THEN  RETURN" 
6230  INPUT  QS 
6240  GIN  @1:H8,H9 

6250  PRINT  "POSITION  PEN  ON  GRAPH  &  ENTER  CHARACTER  TO  BE  PRINTED" 

6260  PRINT  "RETURN  WITH  NO  CHARACTER  TO  EXIT" 

6270  INPUT  Q$ 

6290  IF  Q$=""  THEN  6360 

6295  GIN  @1:H6,H7 

6300  DRAW  @1:H8,H9 

6330  RMOVE  (11:5,-5 

6340  PRINT  @1:Q$: 

6345  MOVE  @1:H6,H7 

6346  IF  Ard=l  THEN  6200 
63  4  7  Ard=l 

6348  GO  TO  6300 
6350  GO  TO  6200 
6360  END 

6370  IF  C9=l  THEN  6440 

6380  OPEN  "C123SIDE";1,"R",ZS 

6390  READ  #1 : S6 

6410  DELETE  Sidel,Side2 

6420  DIM  Sidel(S6) ,Side2(S6) 

6425  FOR  1=1  TO  S6 

6430         READ  #1 : Sidel(I) ,Side2(I) 

6435  NEXT  I 

6436  CLOSE 

6437  Side2=Side2+14 
6440  VIEWPORT  0,24,60,84 
6445  V7=24/112*1000/2 
6450  WINDOW  -V7  ,  V7  ,  0  ,  V7'-'2 
6460  MOVE  (lD:Sidel(l)  ,Side2(l) 
6480  DRAW  @D : Sidel , Side2 

6485  HATCH  ROTATE  0 

6486  HATCH  SPACE  0.1 
6491  HATCH  @D : Sidel , Side2 
6520  IF  C9=l  THEN  6600 

6530  OPEN  "C123END";2,"R",Z$ 
6540  READ  #2:S8 

6570  DIM  Uel(S8) ,Ue2(S8) 

6571  Uel=0 

6572  Ue2=0 

6575  FOR  1=1  TO  S8 

6580  READ  #2:Uel(I) ,Ue2(I) 

6581  NEXT  I 

6582  CLOSE 
6584  Ue2=Ue2+14 
6590  C9=l 

6600  VIEWPORT  G8-12 , 08+12 , 60 , 84 
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6610  REM  WINDOW  -78.26,78.26,0,139. 

6611  WINDOW  -V7,V7,0,V7*2 
6620  MOVE  @D:Uel(l) ,Ue2(l) 
6640  DRAW  @D:Uel,Ue2 

6650  HATCH  @D:Uel,Ue2 

6660  RETURN 

6670  DIM  X5(K,2) 

6690  FOR  1=1  TO  K 

6700  X5(I,1)=A2(I,1) 

6710  X5(I,2)=A2(I,2) 

6730  NEXT  I 

6740  C=K 

6750  N5=K 

6760  GOSUB  5160 

6770  IF  M<0  THEN  6790 

6780  Y1(INT( (500-M) /100+2) )=M 

6790  RETURN 

6800  DIM  X6(K,2) 

6810  FOR  1=1  TO  K 

6820  X6(I,1)=A2(I,1) 

6830  X6(I,2)=A2(I,2) 

6840  NEXT  I 

6850  N6=K 

6860  C=K 

6870  GOSUB  5160 

6880  IF  M<0  THEN  6900 

6890  Y2(INT( (500-M) /100+2) )=M 

6900  RETURN 

6910  DIM  X7(K,2) ,Y7(K,2) 

6920  FOR  1=1  TO  K 

6930  X7(I,1)=A2(I,1) 

6940  X7(I,2)=A2(I,2) 

6950  NEXT  I 

6960  N7=K 

6970  C=K 

6980  GOSUB  5160 

6990  IF  M<0  THEN  7010 

7000  Y3(INT( (500-M)/100+2) )=M 

7  010  RETURN 

7020  DIM  X8(K,2) ,Y8(K,2) 

7030  FOR  1=1  TO  K 

7040  X8(I,1)=A2(I,1) 

7050  X8(I,2)=A2(I,2) 

7060  NEXT  I 

7070  N8=K 

7080  C=K 

7090  GOSUB  5160 

7100  IF  M<0  THEN  7120 

7110  Y4(INT( (500-M) /100+2) )=M 

7120  RETURN 

7130  DIM  X9(K,2) ,Y9(K,2) 

7140  FOR  1=1  TO  K 
7150  X9(I,1)=A2(I,1) 
7160  X9(I,2)=A2(I,2) 

7170  NEXT  I 


7180  N9=K 

7190  C=K 

7200  GOSUB  5160 

7210  IF  M<0  THEN  7230 

7220  H2(INT( (500-M) /100+2) )=M 

7230  RETURN 

7300  PRINT  "J_J_AFTER  WHAT  POSITION?  "; 

7501  INPUT  02 

7502  DELETE  01 

7503  DIM  01(N5+1,2) 

7504  FOR  1=1  TO  02 

7505  01(I,1)=X5(I,1) 

7506  01(I,2)=X5(I,2) 

7507  NEXT  I 

7508  PRINT  "J_ENTER  DOWNRANGE  THEN  HEIGHT  TO  BE  INSERTED 

7509  INPUT  01(02+1,1) ,01(02+1,2) 

7510  FOR  1=02+2  TO  N5+1 

7511  01(I,1)=X5(I-1,1) 

7512  01(I,2)=X5(I-1,2) 

7513  NEXT  I 

7514  DELETE  X5 

7515  DIM  X5(N5+1,2) 

7516  X5=01 

7517  N5=N5+1 

7518  END 

7519  FOR  1=1  TO  5 

7520  01(I,1)=X9(I,1) 

7521  01(I,2)=X9(I,2) 

7522  NEXT  I 

7523  END 

7600  PRINT  "J_J_AFTER  WHAT  POSITION?  "; 

7601  INPUT  02 

7602  DELETE  01 

7603  DIM  01(N6+1,2) 

7604  FOR  1=1  TO  02 

7605  01(I,1)=X6(I,1) 

7606  01(I,2)=X6(I,2) 
760  7  NEXT  I 

7608  PRINT  "J_ENTER  DOWNRANGE  THEN  HEIGHT  TO  BE  INSERTED 

7609  INPUT  01(02+1,1) ,01(02+1,2) 

7610  FOR  1=02+2  TO  N6+1 

7611  01(I,1)=X6(I-1,1) 

7612  01(I,2)=X6(I-1,2) 

7613  NEXT  I 

7614  DELETE  X6 

7615  DIM  X6(N6+1 ,2) 

7616  X6=01 

7617  N6=N6+1 

7618  END 

7700  PRINT  "J_J_AFTER  WHAT  POSITION?  "; 

7701  INPUT  02 

7702  DELETE  01 

7703  DIM  01(N7+1,2) 

7704  FOR  1=1  TO  02 

7705  01(I,1)=X7(I,1) 

7706  01(I,2)=X7(I,2) 

7707  NEXT  I 
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7708  PRINT  "J_ENTER  DOWNRANGE  THEN  HEIGHT  TO  BE  INSERTED?  ' 

7709  INPUT  01(02+1,1) ,01(02+1,2) 

7710  FOR  1=02+2  TO  N7+1 

7711  01(I,1)=X7(I-1,1) 

7712  01(I,2)=X7(I-1,2) 

7713  NEXT  I 

7714  DELETE  X7 

7715  DIM  X7(N7+1,2) 

7716  X7=01 

7717  N7=N7+1 

7718  END 

7800  PRINT  "J_J_AFTER  WHAT  POSITION?  "; 

7801  INPUT  02 

7802  DELETE  01 

7803  DIM  01(N8+1,2) 

7804  FOR  1=1  TO  02 

7805  01(I,1)=X8(I,1) 

7806  01(I,2)=X8(I,2) 

7807  NEXT  I 

7808  PRINT  "J_ENTER  DOWNRANGE  THEN  HEIGHT  TO  BE  INSERTED?  " 

7809  INPUT  01(02+1,1) ,01(02+1,2) 

7810  FOR  1=02+2  TO  N8+1 

7811  01(I,1)=X8(I-1,1) 

7812  01(I,2)=X8(I-1,2) 

7813  NEXT  I 

7814  DELETE  X8 

7815  DIM  X8(N8+1,2) 

7816  X8=01 

7817  N8=N8+1 

7818  END 

7900  PRINT  "J_J_AFTER  WHAT  POSITION?  "; 

7901  INPUT  02 

7902  DELETE  01 

7903  DIM  01(N9+1,2) 

7904  FOR  1=1  TO  02 

7905  01(I,1)=X9(I,1) 

7906  01(I,2)=X9(I,2) 

7907  NEXT  I 

7908  PRINT  "J_ENTER  DOWNRANGE  THEN  HEIGHT  TO  BE  INSERTED?  " 

7909  INPUT  01(02+1,1) ,01(02+1,2) 

7910  FOR  1=02+2  TO  N9+1 

7911  01(I,1)=X9(I-1,1) 

7912  01(I,2)=X9(I-1,2) 

7913  NEXT  I 

7914  DELETE  X9 

7915  DIM  X9(N9+1,2) 

7916  X9=01 

7917  N9=N9+1 

7918  END 

8000  WINDOW  0,1000,0,500 
8010  VIEWPORT  12,124,20,60 
8020  MOVE  @D:0,0 
8030  PRINT  (51,17:0.9,2.1 

8040  PRINT  @D: "J_J_J_J_J_J_'*  INDICATES  COVERAGE  LEVELS  ABOVE 
8050  PRINT  @D:"5  GAL/100  SQ.  FT." 
8060  PRINT  @1 , 17 : 1 . 125,2. 576 
8070  END 
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8200  PAGE 

8210  PRINT   "ENTER  LABEL:  "; 

8220  INPUT  N$ 

8230  MOVE  @l:-720,-410 

8240  PRINT  (ai:N$ 

8250  END 

6.;ou  P(12)  =  1 

8270  GO  TO  2620 

8280  DELETE  X,Y,C1,C2 

8290  DIM  X(N5) ,Y(N5) ,C1(N5) ,C2(N5) 

8300  N=N5 

8310  07=2 

8320  FOR  1=1  TO  N 

8330  X(I)=X5(I,1) 

8340  Y(I)=X5(I,2) 

8350  NEXT  I 

8360  GO  TO  2630 

8370  DELETE  X,Y,C1,C2 

8380  DIM  X(N6) ,Y(N6) ,C1(N6) ,C2(N6) 

8390  N=N6 

8400  07=3 

8410  FOR  1=1  TO  N 

8420  X(I)=X6(I,1) 

8430  Y(I)=X6(I,2) 

8440  NEXT  I 

8450  GO  TO  2630 

8460  DELETE  X,Y,C1,C2 

8470  DIM  X(N7) ,Y(N7) ,C1(N7) ,C2(N7) 

8480  N=N7 

8490  07=4 

8500  FOR  1=1  TO  N 

8510  X(I)=X7(I,1) 

8520  Y(I)=X7(I,2) 

8530  07=4 

8540  NEXT  I 

8550  GO  TO  2630 

8560  DELETE  X,Y,C1,C2 

8570  DIM  X(N8) ,Y(N8) ,C1(N8) ,C2(N8) 

8580  N=N8 

8590  07=5 

8600  FOR  1=1  TO  N 

8610  X(I)=X8(I,1) 

8620  Y(I)=X8(I,2) 

8630  NEXT  I 

8640  GO  TO  2630 

8650  DELETE  X,Y,C1,C2 

8660  DIM  X(N9) , Y(N9) ,C1(N9) ,C2(N9) 

8670  N=N9 

8680  07=6 

8690  FOR  1=1  TO  N 

8700  X(I)=X9(I,1) 

8710  Y(I)=X9(I,2) 

8720  NEXT  I 

8730  GO  TO  2630 

8740  END 

50000  MOVE  X5(l,l) ,X5(1,2) 
50010  FOR  1=1  TO  N5 
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50020        DRAW  X5(I,1) ,X5(I,2) 
50030  NEXT  I 

500A0  MOVE  X6(l,l) ,X6(1,2) 
50050  FOR  1=1  TO  N6 
50060         DRAW  X6(I,1) ,X6(I,2) 
50070  NEXT  I 

50080  MOVE  X7(l,l) ,X7(1,2) 
50090  FOR  1=1  TO  N7 
50100         DRAW  X7(I,1) ,X7(I,2) 
50110  NEXT  I 

50120  MOVE  X8C1,1) ,X8(1,2) 
50130  FOR  1=1  TO  N8 
50140         DRAW  X8(I,1) ,X8(I,2) 
50150  NEXT  I 

50160  MOVE  X9(1,1),X9(1,2) 

50170  FOR  1=1  TO  N9 

50180         DRAW  X9(I,1) ,X9(I,2) 

50190  NEXT  I 

50200  END 

51000  REM  CHANGE  X5 

51010  PAGE 

51020  FOR  1=1  TO  N5 

51030  PRINT  USING  "3D , 2( 6D) " : I ,X5(I , 1 ) ,X5( I , 2) 
51040  NEXT  I 

51050  PRINT  "WHICH  ELEMENT  ENTER  ROW  THE  COL"; 

51055  PRINT   "     ENTER  0  TO  EXIT  "; 

51060  INPUT  Row, Col 

51065  IF  Row=0  THEN  51200 

51070  PRINT  "ENTER  CHANGE"; 

51080   INPUT  Chg 

51090  X5(Row,Col)=X5(Row,Col)+Chg 
51100  GO  TO  51050 
51200  END 

52000  REM  CHANGE  X6 

52010  PAGE 

52020  FOR  1=1  TO  N6 

52030  PRINT  USING  " 3D , 2 ( 6D ) " : I , X6 ( I , 1 ) , X6 ( I , 2 ) 
52040  NEXT  I 

52050  PRINT   "WHICH  ELEMENT  ENTER  ROW  THE  COL"; 

52055  PRINT  "     ENTER  0  TO  EXIT  " : 

52060  INPUT  Row, Col 

52065  IF  Row=0  THEN  51200 

52070  PRINT  "ENTER  CHANGE"; 

52080  INPUT  Chg 

52090  X6 (Row, Col)=X6(Row, Col )+Chg 

52100  GO  TO  52050 

53000  REM  CHANGE  X7 

53010  PAGE 

53020  FOR  1=1  TO  N7 

53030  PRINT  USING  " 3D , 2 ( 6D) " : I , X7 ( I , 1 ) , X7 ( I , 2) 
53040  NEXT  I 

53050  PRINT  "WHICH  ELEMENT  ENTER  ROW  THE  COL"; 

53055  PRINT  "     ENTER  0  TO  EXIT  "; 

53060  INPUT  Row, Col 

53065  IF  Row=0  THEN  51200 

53070  PRINT  "ENTER  CHANGE"; 

53080  INPUT  Chg 


53090  X7 (Row , Col ) =X7 (Row, Col)+Chg 

53100  GO  TO  53050 

54000  REM  CHANGE  X8 

54010  PAGE 

54020  FOR  1=1  TO  N8 

54030  PRINT  USING  " 3D , 2 ( 6D) " : I ,X8 ( I , 1) ,X8 ( I , 2) 
54040  NEXT  I 

54050  PRINT  "WHICH  ELEMENT  ENTER  ROW  THE  COL"; 

54055  PRINT  "     ENTER  0  TO  EXIT  "; 

54060  INPUT  Row, Col 

54065  IF  Row=0  THEN  51200 

54070  PRINT  "ENTER  CHANGE"; 

54080  INPUT  Chg 

54090  X8(Row,Col)=X8(Row,Col)+Chg 

54100  GO  TO  54050 

55000  REM  CHANGE  X9 

55010  PAGE 

55020  FOR  1=1  TO  N9 

55030  PRINT  USING  " 3D , 2 ( 6D ) " : I ,X9 ( I , 1) , X9 ( I , 2 ) 
550  40  NEXT  I 

55050  PRINT  "WHICH  ELEMENT  ENTER  ROW  THE  COL": 

55055  PRINT  "     ENTER  0  TO  EXIT  " ; 

55060  INPUT  Row, Col 

55065  IF  Row=0  THEN  51200 

55070  PRINT  "ENTER  CHANGE" ; 

55080  INPUT  Chg 

55090  X9 ( Row , Co 1) =X9 ( Row , Co 1) +Chg 

55100  GO  TO  55050 

60000  REM  PRINT  X  ON  FOOT 

60010  FOR  1=1  TO  N5 

60020         MOVE  @D:X5(I,1) ,X5(I,2) 

60030         GOSUB  62000 

60040  NEXT  I 

60050  FOR  1=1  TO  N6 

60060        MOVE  @D:X6(I,1) ,X6(I,2) 

60070        GOSUB  62000 

60080  NEXT  I 

60090  FOR  1=1  TO  N7 

60100         MOVE  @D:X7(I,1) ,X7(I,2) 

60110         GOSUB  62000 

60120  NEXT  I 

60130  FOR  1=1  TO  N8 

60140         MOVE  @D:X8(I,1) ,X8(I,2) 

60150         GOSUB  62000 

60160  NEXT  I 

60170  FOR  1=1  TO  N9 

60180         MOVE  @D:X9(I,1) ,X9(I,2) 

60190         GOSUB  62000 

60200  NEXT  I 

60210  END 

62000  RDRAW  @D : 3 , 5 
62010  RDRAW  @D:-6,-10 
62020  RMOVE  @D ; 3  ,  5 
62030  RDRAW  @D:-3,5 
62040  RDRAW  @D:6,-10 
62050  RETURN 
63000  FOR  1=1  TO  6 
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63010         MOVE  @D:X1(I),Y1(I) 

63015         IF  X1(I)=0  AND  Y1(I)=0  THEN  63030 

63020         GOSUB  63700 

63030  NEXT  I 

63100  FOR  1=1  TO  6 

63110  MOVE  @D:X2(I),Y2(I) 

63115         IF  X2(I)=0  AND  Y2(I)=0  THEN  63130 

63120        GOSUB  63700 

63130  NEXT  I 

63200  FOR  1=1  TO  6 

63210         MOVE  @D : X3 ( I) , Y3  C I ) 

63215         IF  X3(I)=0  AND  Y3(I)=0  THEN  63230 

63220        GOSUB  63700 

63230  NEXT  I 

63300  FOR  1=1  TO  6 

63310         MOVE  (3D:X4(I)  ,Y4(I) 

63315         IF  X4(I)=0  AND  Y4(I)=0  THEN  63330 

63320         GOSUB  63700 

63330  NEXT  I 

63400  FOR  1=1  TO  6 

63410         MOVE  @D:H1(I) ,H2(I) 

63415         IF  H1(I)=0  AND  H2(I)=0  THEN  63430 

63420        GOSUB  63700 

63430  NEXT  I 

63500  END 

63700  RDRAW  @D : 1 ,  5 
63710  RDRAW  (aD:-2,-10 
63720  RMOVE  @D : 1  .  5 
63730  RDRAW  (30: -1,5 
63740  RDRAW  @D:2, -10 
63750  RETURN 
63760  END 

64000  PRINT  "WHAT  FILE  FOR  RETRIEVAL?  "; 
64005  INPUT  B$ 

64008  DELETE  XI , X2 , X3 , X4 , X5 , X6 , X7 , X8 , X9 , HI , H2 , Yl , Y2 , Y3 , Y4 

64009  DIM  Xl(6) ,X2(6) ,X3(6) ,X4(6) ,H1(6) ,H2(6) ,Y1(6) ,Y2(6) , Y3(6) ,Y4(6) 

64010  OPEN  B$: 1 , "R" ,CS 

64020  READ  #1 : N5 , N6 , N7 , N8 , N9 , XI , X2 , X3 , X4 , HI , Yl , Y2 , Y3 , Y4 , H2 
64030  DIM  X5(N5,2) ,X6(N6,2) ,X7(N7,2) ,X8(N8,2) ,X9(N9,2) 
64040  READ  #1 : X5 , X6 , X7 , X8 , X9 
64050  PRINT  "TENSION?  "; 
64060  INPUT  P(2) 
64070  END 

65000  PRINT  "WHAT  FILE  FOR  STORAGE?  "; 

65005  INPUT  B$ 

65006  CALL  "FILE" , 0 , B$ , C$ 

65007  PRINT  LEN(C$) 

65008  IF  LEN(C$)<>0  THEN  65020 
65010  CREATE  B$;1000,0 

65020  OPEN  B$;l, "F" ,C$ 

65030  WRITE  #1 : N5 , N6 , N7 , N8 , N9 ,X1 , X2 ,X3 , X4 , HI , Yl , Y2 , Y3 , Y4 , H2 
65040  WRITE  #1 : X5 , X6 , X7 , X8 , X9 
65060  END 
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optimum  retardant  distribution  from  a  specific  airplane  and  tanking  system.  Text  and 
drawings  are  assembled  into  a  sample  guide  for  the  Evergreen  P2V-5  airplane. 
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The  Inlermountain  Research  Station  provides  scientific  knowledge  and  technology  to  im- 
prove management,  protection,  and  use  of  the  forests  and  rangelands  of  the  Intermountain 
West.  Research  is  designed  to  meet  the  needs  of  National  Forest  managers,  Federal  and 
State  agencies,  industry,  academic  institutions,  public  and  private  organizations,  and  individu- 
als. Results  of  research  are  made  available  through  publications,  symposia,  workshops, 
training  sessions,  and  personal  contacts. 

The  Intermountain  Research  Station  territory  includes  Montana,  Idaho,  Utah,  Nevada,  and 
western  Wyoming.  Eighty-five  percent  of  the  lands  in  the  Station  area,  about  231  million 
acres,  are  classified  as  forest  or  rangeland.  They  include  grasslands,  deserts,  shrublands, 
alpine  areas,  and  forests.  They  provide  fiber  for  forest  industries,  minerals  and  fossil  fuels  for 
energy  and  industrial  development,  water  for  domestic  and  industrial  consumption,  forage  for 
livestock  and  wildlife,  and  recreation  opportunities  for  millions  of  visitors. 

Several  Station  units  conduct  research  in  additional  western  States,  or  have  missions  that 
are  national  or  international  in  scope. 

Station  laboratories  are  located  in: 

Boise,  Idaho 

Bozeman,  Montana  (in  cooperation  with  Montana  State  University) 
Logan,  Utah  (in  cooperation  with  Utah  State  University) 
Missoula,  Montana  (in  cooperation  with  the  University  of  Montana) 
Moscow,  Idaho  (in  cooperation  with  the  University  of  Idaho) 
Ogden,  Utah 

Provo,  Utah  (in  cooperation  with  Brigham  Young  University) 
Reno,  Nevada  (in  cooperation  with  the  University  of  Nevada) 

USDA  policy  prohibits  discrimination  because  of  race,  color,  national  origin,  sex,  age,  reli- 
gion, or  handicapping  condition.  Any  person  who  believes  he  or  she  has  been  discriminated 
against  in  any  USDA-related  activity  should  immediately  contact  the  Secretary  of  Agriculture, 
Washington,  DC  20250. 


